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CHAPTER  I 


INTRODUCTION 


Transmission  of  infrared  radiation  through  the  atmosphere  is 
adversely  affected  by  turbulence,  scattering  by  aerosols,  and  mole- 
cular absorption.  This  study  is  restricted  to  certain  aspects  of 
the  last  problem. 

An  investigation  of  a low  resolution  solar  spectrum  from  1 
to  15  ym, such  as  that  shown  in  Fig.  1[1],  reveals  that  there  are 
certain  regions  of  the  spectrum  where  there  is  relatively  little 
absorption  by  atmospheric  gases.  One  of  these  regions  from  3.5  to 
4.0  microns  corresponds  to  the  wavelength  range  of  the  deuterium  fluor- 
ide (OF)  laser.  This  correspondence  of  DF  laser  wavelengths  with 
an  atmospheric  window,  along  with  the  high  power  capability  of  this 
laser  makes  it  very  attractive  for  certain  applications  such  as  com- 
munications or  laser  radar. 

Although  the  molecular  absorption  in  the  pF  laser  region  is 
small,  it  does  vary  rapidly  with  frequency.  This  may  be  seen  more 
clearly  in  the  medium  resolution  spectrum  obtained  by  Streete[2] 
over  a 24.9  km  sea-level  path  (Fig.  2).  Straight  lines  correspond- 
ing to  selected  DF  laser  lines  have  been  superimposed  on  the  spec- 
trum in  Fig.  2. 

In  order  for  intelligent  decisions  to  be  made  regarding  the 
use  of  DF  lasers  in  systems  involving  transmission  of  radiation 
through  the  atmosphere,  a precise  knowledge  of  the  atmospheric 
transmittance  at  the  laser  frequencies  is  required.  As  may  be  seen 
in  Fig.  2 this  precise  knowledge  is  very  difficult  to  obtain  using 
conventional  spectroscopy  techniques  since  the  laser  lines  frequent- 
ly occur  at  frequencies  where  the  transmittance  is  changing  very 
rapidly  with  frequency.  Even  the  highest  resolution  conventional 
spectrometers  tend  to  smooth  this  rapid  variation  because  of  their 
finite  resolution.  Thus,  if  precise  knowledge  of  the  atmospheric 
transmittance  at  a laser  frequency  is  needed,  it  it  necessary  to 
measure  it  directly  using  the  laser  as  the  source.  This  is  the 
approach  taken  in  this  study. 

It  can  be  seen  from  Fig.  1 that  there  are  contributions  to 
molecular  absorption  in  the  DF  laser  region  from  HDD,  methane, 
nitrous  oxide,  and  carbon  dioxide.  In  addition  there  is  a contri- 
bution from  water  continuum  absorption  presumably  caused  by 


Fig.  1.  Low  resolution  solar  spectrum  compared  with 
laboratory  spectra  of  atmospheric  gasses 
(Shaw) , 


the  wings  of  the  strong  water  vapor  bands  at  2.6  and  6.3  microns, 
and  a small  pressure-induced  nitrogen  continuum. 

The  purpose  of  this  study  was  to  measure  accurately  the  ab- 
sorption of  selected  DF  laser  lines  by  HDO,  CH^,  N2O,  CO2  and  the 
water  continuum. 

The  absorption  was  measured  for  each  constituent  separately 
using  a long  path  laboratory  absorption  cell  and  a small,  pulsed, 
single-frequency  laser.  DF  laser  lines  selected  for  study  were 
the  2-1  P(6),  P(7),  P(8),  P(10),  and  P(ll)  lines  and  the  3-2  P(6), 
P(7),  and  P(8)  lines.  These  lines  are  reported  to  be  relatively 
intense  in  high-power  lasers.  Ideally  one  would  wish  to  study  the 
absorption  on  all  lines.  However,  the  lasers  used  in  this  study 
were  capable  of  operation  on  about  30  lines,  and  it  would  have  been 
prohibitively  costly  in  terms  of  time  and  money  to  measure  the  ab- 
sorption on  all  lines. 

McClatchey,  et  al.[3  ] have  compiled  line  data  on  the  mole- 
cular absorption  for  the  most  important  absorbers  in  the  infrared 
region!  and  have  made  these  data  available  to  other  workers  on  mag- 
netic tape.  The  data  are  based  on  conventional  spectra  from  many 
different  sources.  Using  this  tape  it  is  possible  to  compute  the 
atmospheric  absorption  at  desired  frequencies.  The  accuracy  of 
such  a computation  is  of  course  limited  by  the  accuracy  of  the  data 
on  which  it  is  based.  However,  even  if  the  computation  is  off  by  a 
factor  of  two  or  more,  it  is  still  very  useful  as  an  aid  in  planning 
the  necessary  absorber  concentration  and  path  required  to  make 
accurate  measurements  of  the  absorption  at  the  laser  frequencies  of 
interest.  This  was  the  approach  taken  in  this  study. 

In  Chapter  II  the  theoretical  basis  of  molecular  absorption 
calculations  is  discussed,  computer  programs  which  calculate  the 
absorption  are  described  and  calculations  are  presented  for  the 
DF  laser  lines.  Also  discus  ed  in  Chapter  II  are  measurements  made 
by  Burch  [4]  of  the  water  vapor  and  nitrogen  continue  using  a 
conventional  spectrometer  and  a long  path  absorption  cell. 

In  Chapter  III  the  experimental  equipment  used  to  make  the 
absorption  measurements  is  described. 

In  Chapter  IV  the  absorption  measurements  are  described  and 
the  results  are  discussed  and  compared  with  the  experimental  results 
of  other  workers  and  with  the  computer  predictions-  Also  described 
is  a laser  spectroscopy  technique  which  can  be  used  to  measure  the 
frequency  separation  of  an  absorption  line  and  a nearly  coincident 
laser  line. 
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CHAPTER  II 


MOLECULAR  ABSORPTION  CALCULATIONS 


For  many  applications  it  is  useful  to  compute  both  single 
frequency  absorption  and  absorption  spectra,  for  various  situations 
of  interest,  from  spectral  data.  The  most  complete  set  of  standard 
spectral  data  currently  available  is  the  AFCRL  Line  Compilation[3] . 
This  compilation  lists  data  for  over  100,000  absorption  lines  from 
0,307  cm'l  to  15,000  cm*i , and  is  available  on  magnetic  tape,  In 
this  chapter  the  theoretical  background  for  the  absorption  calcu- 
lations will  be  presented,  two  computer  programs  which  have  been 
written  to  compute  fixed  frequency  absorption  and  absorption  spec- 
tra will  be  discussed,  and  results  of  calculations  made  for  the 
OF  laser  region  of  the  spectrum  will  be  presented. 


A.  Theoretical  Background 

Assuming  the  photon  density  is  low  enough  that  nonlinear  ex- 
tinction processes  do  not  occur,  the  change  in  intensity  dl  of  a 
parallel  beam  of  electromagnetic  radiation  at  frequency  v travers- 
ing an  infinitesimal  path  dt  is  given  by  Lambert's  law 


(1)  dl^  = - k(v)I^  da 

where  da  is  the  amount  of  absorber  in  a volume  of  unit  cross- 
sectional  area  and  length  d£.  The  factor  k(v)  is  called  the  ab- 
sorption coefficient. 

The  transmittance  T through  an  absorber  amount  a is  found 
by  integrating  Eq.  (1). 

(2)  T = — = g"k(v)a 

V jO 


Here  I is  the^intensity  of  the  beam  after  passing  through  the 
absorber  and  I is  the  intensity  of  the  beam  before  it  enters  the 
absorber, 
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The  absorber  amount  a and  absorption  coefficient  k(v)  may  in 
general  be  expressed  in  any  convenient  units.  The  requirement  is 
that  the  product  k(v)a  be  dimensionless.  In  making  calculations 
using  the  AFCRL  Line  Compilation  it  is  convenient  to  express  the 
absorber  amount  in  mol  ;m-2.  The  absorption  coefficient  is  there- 
fore in  (mol-cm"^)-1  . In  the  experimental  work  presented  later  it 
is  conveni'^nt  to  characterize  the  absorber  amount  by  the  length  of 
the  path  cjntaining  the  absorber.  The  absorber  amount  is  then 
given  in  km  and  the  absorption  coefficient  in  km"'.  A more  com-_ 
plete  discussion  of  the  units  and  conversion  factors  between  various 
units  is  given  in  Appendix  A. 

The  absorption  of  electromagnetic  radiation  by  an  absorbing 
medium  is  caused  by  the  interaction  of  the  electromag”*etic  fields 
with  the  energy  states  of  the  absorber.  The  energy  cf  a molecule 
consists  of  its  electronic,  vibrational,  and  rotational  energies 
and  the  interaction  energies  between  these  three  energy  modes. 

Energy  levels  are  quantized  so  that  the  molecule  can  have  only  certain 
allowed  values  of  energy.  The  absorption  of  one  photon  of  electro- 
magnetic radiation  then  involves  a change  of  energy  of  one  molecule 
from  one  energy  state  to  a higher  energy  state, ignoring  Raman  and 
multiphoton  processes.  The  difference  in  the  energy  states  must 
be  equal  to  the  energy  of  the  photon  h’.'.  In  addition  for  a given 
type  of  molecule  only  transitions  between  certain  energy  states 
are  allowed,  and  therefore  only  photons  having  energy  equal  to  the 
difference  in  energy  betweei  those  states  will  be  absorbed.  Thus, 
for  a given  absorber  type,  electromagnetic  radiation  will  be  ab-_ 
sorbed  at  only  chose  frequencies  where  the  photon  energy  is  equal 
to  the  energy  difference  of  an  allowed  transition.  The  above  dis- 
cussion applies  equallv  „ell  to  emission  of  radiation. 

Allowed  electronic  energy  states  are  normally  widely  separated, 
and  the  absorption  or  emission  of  radiation  due  to  changes  in 
electronic  energy  usually  occurs  in  the  visible  or  ultra-violet 
portion  of  the  spectrum.  Absorption  or  emission  in  tiie  infrared 
portion  of  the  spectrum  of  interest  in  this  study  is  caused  by 
changes  in  both  the  vibrational  and  rotational  energies  of  the  mole- 
cule while  the  electronic  energy  remains  unchanged.  Changes  in 
rotational  energy  only  result  in  emission  or  absorption  in  the  far- 
infrared  to  microwave  region  of  the  spectrum. 

Transitions  occur  from  a higher  energy  state  u to  a lower 
energy  state  x,  at  a rate 


(3)  * ''u*’ 

where  B is  the  Einstein  coefficient  For  stimulated  emission, 
p(v)  is^lhe  radiation  density  per  unit  frequency  interval,  and  A^^^ 
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is  the  Einstein  coefficient  for  spontaneous  emission.  The  transi- 
tion rate  from  lower  level  to  upper  level  is  given  by 


(4) 


“ju  ' 


-1 


where  is  the  Einstein  coefficient  for  absorption.  In  consider 
ing  Planck's  radiation  law  Einstein  showed 


(5) 


B = B — 
£u  uii  g 


where  g and  g are  the  degeneracies  of  the  upper  and  lower  states 
respectively,  and 


(6) 


-1 


(7) 


The  transition  rates  may  now  be  written  in  terms  of  A : 

U ■Jv 


U = f t l"!  A 


sec 


-1 


(8)  W = -y-  A sec 

h 8.hv’ 


-1 


The  transition  rates  are  related  to  the  transition  moment  by 


(9) 


fia  4 3 

. _ 0471  V 

uz  3hg„ 


-1 

sec 


where  m is  the  electric  dipole  transition  momen;t  <E^-|m|E^>.  It 
is  necessary  to  know  |E^>  and  |E^^  to  calculate  |m^£|.  Sihce  the 
states  have  a finite  lifetime,  the  uncertainty  principle  says  it 
is  not  possible  to  know  the  energy  of  the  energy  levels  pre- 
cisely. In  addition,  for  a system  of  colliding  molecules,  the 
collision  interaction  energy  must  be  included  for  the  whole  system. 
This  is  a many  body  problem  which  cannot  be  solved.  The  effect  of 
these  considerations  is  that  the  energy  levels  are  smeared  and 
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there  is  a finite  transition  probability  for  interaction  with  the 
radiation  field  over  a range  of  frequencies. 

This  effect  is  taken  into  account  by  introducing  a frequency 
dependent  transition  probability.  The  transition  rates  may  then 
be  written. 


(10) 


us, 


\ 8tt  hv'^  / 


and 


(11) 


s,u 


^S,  8tr  hv 


where 


and  g(v-v^)  is  a normalized  shape  factor. 

If  q(v-v  ) is  small  except  for  a narrow  range  of  frequencies 

I 0 

near  v^,  and 

p(v) 

3 

V 

is  essentially  constant  over  this  range,  then  the  transition  rates 
are 


(12) 

and 


p(vq)c^ 

Sirhv^ 

0 


(13) 


g p(v„)c^ 
_y_  " 0 
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where  the  spontaneous  emission  term  in  W has  been  ignored  (since  its 
contribution  to  a collimated  beam  would  Ve  negligible)  and 


(14) 


-00 

g(v-Vy)  dv  = 1 

*00 


When  p(v)  is  monochromatic. 


(15)  p(v')  = Pq  6 (v-v')  , 


and  Eqs.  (11)  and  (10)  become 
9..  Pg  ^ 


(16) 

and 

(17) 


W„  = -2 ^ g(v-v J 

9j,  Stt  hv^  ° 


P-  c 

8, 


respectively. 


The  number  of  molecules  per  unit  absorber  in  the  upper  energy 
level  is  related  to  the  number  of  molecules  per  unit  absorber  in 
the  lower  energy  level  by  the  Boltzmann  distribution 


(18) 


g (E  -E^) 
N = N ^ e kT 
u £ g. 


A monochromatic  beam  passing  through  an  infinitesimal  amount  of 
absorber  da  will  have  an  intensity  change 


(19)  dl^  = (N^  hv  da 


or,  from  Eqs.  (13)  and  (18) 
( '"’o 


(20) 


dl  = 

V 


kT 


V 


- 1 I Nj  “w  *"> 
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From  Eqs.  (1 ) and  (20) 


(21)  dl  = - k(v)I  , da  = e 


N^W^„h„da  . 


Equation  (21)  may  be  solved  for  the  absorption  coefficient 

(22) 


k(v)  = 1 ~ e 


hv  \ 
_JO  \ 
kT 


/ 


N.  hv 

I £U 


I 


Using  I = c p and  W from  Eq.  (16) 

\)  vJ  X>  LI 


9m  a c 

(23)  k(v)  = N / 

^ 8ir  V 


I ""o\ 


1 - e 

V 


kT 


g(v-VQ) 


/ 


The  line  strength  S is  defined  as  the  integral  over  all  v of  k(v)dv 

2 


(24) 


S = 


g A c 

k(v)d  v=  N,  -M.  II  - e 


kT 


i g„  8it 


r 9(v-v„) 


dv 


-00  V 


Again  assuming  g(v-v  ) is  small  except  for  a narrow  range  of  fre- 
quencies about  vqi  tne  integral  in  Eq,  (24)  becomes 


(25) 


r . 1 

dv  ■=  : 


2 


The  line  strength  is  then 


g A „ c‘ 

(26)  S = N„  - 


^ g 


Z 8ir  V, 


1 - e 

V 


kT 


y 
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and  the  extinction  coefficient  is  simply 


(27)  k(v)  = S g(v-v^) 


If  vibration-rotation  interaction  is  neglected,  N is  related 
to  the  total  number  of  molecules  of  absorber  by  the  rotational  and 
vibrational  partition  functions  and  Qy.  The  relationship  of 
to  Ny  is  given  by 


(28) 


Qv  Qr 


q e 


kT 


The  line  strength  then  becomes 


The  line  strength  per  molecule  per  square  centimeter  S°  at 
standard  temperature  Tq  is  found  by  evaluating  Eq.  (29)  at  tempera- 
ture Tq  and  dividing  by  the  number  density  Nj. 


(30) 


S° 


vv 


3it  V ^ 
0 


The  line  strength  at  any  other  temperature  is  then 


The  AFCRL  line  compilatic  [3]  gives  for  each  line  the  line 
strength  at  296°K,  Equation  (3  ) is  sed  to  correct  the  line 
strength  for  use  in  making  calculations  at  other  temperatures. 

The  vibrational  partition  function  Qy  can  be  determined  as 
fol 1 ows . 
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For  the  Maxwel  1 -Boltzmann  distribution,  the  number  of  mole- 
cules in  a vibrational  state  v,  N^,  is  prooortional  to  the  Boltz- 
mann factor 

-Ey/kT 

e 

For  the  vibrational  states  the  term  values  G(v)  are  given  b>[5] 

(32)  G(v)  = ii)g(v  + %)  - + %)^  + 

where 

E T 

(33)  G(v)  = cm  ' 

Now  let 

Gq(v)  = G(v)  - G(0)  % o)g  V , 

assuming  that  the  simple  harmonic  oscillator  function  approximation 
is  adequate.  Then 

- G (v)  hc/kT 

(34)  ® 

where  C,  is  a constant  of  proportionality.  Again  assuming  the 
validity  of  the  simple  harmonic  oscillator  approximation 


(35) 

- Gn(l)  hc/kT 

(36) 

N,  = 

u 

e 

- G.(2)  hc/kT 

(37) 

"2  = 

'^1 

e ° 

etc. 

From  Eq. 

(34) 
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1 

-0)  hc/kT 


Now  (uhc  is  approximately  the  difference  between  the  vibrational 
energy  levels  of  the  absorbing  transition,  and  since  the  difference 
in  rotational  energy  levels  is  small  it  is  also  approximately  the 
energy  of  the  absorbed  photon.  Therefore 


McClatchey,  et  al[3]  have  given  a table  of  Q for  the  various 

, . . . T • I 4-U_Vn  r-f.ni  1 ; rr.m_ 


molecules  whose  absorption  lines  are  listed  in  the  AFCRL  Line  Com 
pilation  which  is  repeated  here  as  Table  1. 


TABLE  1 

VIBRATIONAL  PARTITION  FUNCTION  VALUES  FROM  McCLATCHEY[3] 


Molecule  175 

200 

Temperature 
225  250 

°K 

275 

296 

325 

H2O 

1 .000 

1 .000 

1.000 

1.000 

1 .000 

1 .000 

1 .001 

CO, 

1 .0095 

1 .0192 

1 .0327 

1.0502 

1 .0719 

1 .0931 

1 .1269 

“3 

1 .004 

1 .007 

1 .013 

1 .022 

1 .033 

1 .046 

1 .066 

NpO 

1 .017 

1 .030 

1 .048 

1 .072 

1 .100 

1 .127 

1 .170 

CO 

1 .000 

1 .000 

1 .000 

1 .000 

1 .000 

1 .000 

1 .000 

CH4 

1 .000 

1.000 

1 .001 

1 .002 

1 .004 

1 .007 

1 .011 

02 

1.000 

1 .000 

1 .000 

1 .000 

1 .000 

1 .000 

1 .001 

From  Table  1 it  appears  that  may  be  ignored  for  H^O,  CO, 
and  CH4.  For  the  other  molecules  it  appears  that  Eq.  (44)  is  not 
a very  good  approximation.  The  table  indicates  no  dependence  on 
frequency  although  the  above  discussion  indicates  that  there  is  a 
definite  frequency  dependence.  Therefore  it  would  seem  that  cal- 
culations made  at  frequencies  lower  tnan  about  800  cm"'  would  have 
reduced  accuracy. 


The  rotational  partition  function  is  defined  by 

-ER(J)/kT 

(45)  g,  = I 9j  e 


-E^(J)/kT 


The  degeneracy  g is  different  for  different  types  of  molecules. 
For  a linear  molecule  (N2O,  C02>  CO)  Herzberg[5]  gives  as  an 


approximation 
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(46) 


Qr(T)  = 


kT 

hcB 


so  that 


For  asymmetric  top  molecules  {H2O,  0^)  Herzberg[5]  gives  as  an 
approximation 


For  spherical  top  molecules  (CH^)  Herzberg[5]  gives 


The  line  strength  correction  may  now  be  written 


where  BX  is  either  1.0  or  1.5  depending  on  the  molecule. 
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There  are  three  effects  w'lich  contribute  to  the  line  shape 
factor  g(v-vo).  The  first  of  these  is  called  natural  broadening 
and  results  from  the  smearing  of  energy  levels  due  to  the  Heisen- 
berg uncertainty  principle  since  the  energy  states  have  a finite 
lifetime.  At  temperatures  and  pressures  normally  occuring  in  the 
atmosphere  this  effect  is  negligible  compared  to  other  mechanisms 
and  will  not  be  considered  further. 


The  second  effect  is  called  Doppler  broadening  and  is  caused 
by  the  rapid  movement  of  the  molecules  while  they  are  emitting  or 
Absorbing.  The  line  shape  factor  due  to  Doppler  broadening  is 
given  by[6] 


(53) 


g(v-v^) 


-[x,n2(v-VQ)^/oiQ^] 

e 


where  ap,  is  the  Doppler  half-width  defined  as  half  the  width  of 
the  g(v-vo)  versus  v curve  at  half  the  maximum  amplitude,  is 
given  by 


(54) 


1/2 


where  vq  is  the  center  frequency  of  the  line,  k is  Boltzmann's 
const^.nt,  T is  temperature  in  degrees  Kelvin,  m is  the  molecular 
mass,  and  c is  the  speed  of  light.  For  convenience  of  calculation 
this  may  be  rewritten 


-7  /tA 

(55)  =3.5812  X 10  ' 


where  M is  the  molecular  weight. 

The  third  effect  is  caused  by  collisions  between  rapidly  moving 
molecules  in  the  absorbing  medium  and  is  called  collision-broadening. 
This  is  the  predominant  effect  at  pressures  higher  than  50_to  100 
torr.  For  frequencies  near  Vq  the  Lorentz  line  shape[7]  gives  a 
reasonable  representation  of  the  line  shape  factor  g(v-v^). 


(56)  g(v-v^j)  = - 


( ^2  . 2 
(v-Vp)  + a|_ 
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Here  a,  is  the  Lorentz  half-width  defined  as  half  the  width  of  the 
glv-Vg)  versus  v curve  at  one  half  maximum  amplitude. 

The  AFCRL  line  compilation[3]  gives  for  each  line,  the  Lorentz 
half-width  at  296°K  and  1 atmosphere  total  pressure.  For  calcula- 
tions at  other  temperatures  and  pressures  the  half-width  must  be 
modified.  The  form  of  the  modification  may  be  determined  fron.  a 
consideration  of  the  Lorentz  half-width  given  by  kinetic  theory[8] 


where  F is  the  collision  frequency,  N^-  is  the  number  of  molecules 
of  the  ith  type  per  unit  volume,  ^ is  the  sum  of  the  optical 
collision  diameters  of  the  absorbini^i  molecule  and  a molecule  of  the 
ith  type,  k is  Boltzmann's  constant,  T is  the  temperature  in  degrees 
K,  m,  is  the  mass  of  the  absorbing  molecule,  and  m,-  is  the  mass  of 
a moTeculo  of  the  ith  type. 

For  a binary  mixture  of  an  absorbing  gas  a and  a broadening 
gas  b,  Eq.  (57)  reduces  to 


(58) 


or 


(58a) 


'L  ■ Att 


+ N.  (D  . ) 
b a,b' 


“L  ■ i 87t  J 

V f 


-- 

•Na<“a,a>' 

a. 

1/2 


/ m +m,  \ 

1/2" 

2i.kT  “ 1 

1—  _J 

> 

) 

Na(‘>a.a)'(rT'^V“a,b) 


,2  /V\ 

^Vb 


I/2I 


Now  define  and  B as  follows 
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(67) 


From  the  gas  law  the  partial  pressure  of  gas  of  the  ith  type  is 
(62)  P^.  = N.  kT 

Using  Eqs.  (59),  (60),  and  (62)  Eq.  (58a)  becomes 

/ 1 V/2 

(«)  “L  = t'’a  * Pb  <^abJ  ' 

Substituting  B from  Eq.  (61)  into  Eq.  (63)  yields 


( 1 V/2 

“L  = (^j)  <=ab  8 * V 
The  total  pressure  P is  + P|^  so  Eq.  (64)  may  be  rewritten 

- i^V'" 

(65)  “L  ■ \^8»kTj  C^lj  [P  + P3(b-1)] 

or 

/ 1 

(66)  “L  " Utt  klj  ^ab  ^e 


where  the  equivalent  pressure  P^  may  be  defined 


(67)  Pg=P+Pg(B-l) 

The  modification  to  the  Lorentz  half-width  a|_Q  given  in  th 
AFCRL  line  compilation  for  different  temperature  and  pressure  is 
then 


(68) 


» 
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where  P is  one  atmosphere  and  T is  296°K.  The  exponent  on  T /T 
is  giveR  ns  CX  rather  than  1/2  since  it  has  been  determined  experi- 
Fr«ntal1y  that  CX  is  actually  about  .62  for  H2O  [9]  and  .58  for 
C02-  The  value  of  .58  for  CO2  is  quoted  by  Deutschman  and  Cal  fee 
[lO]  as  a private  communication  from  Benedict. 

B is  called  the  self-broadening  coefficient  and  represents 
the  ratio  of  the  effect  on  the  half-width  of  collisions  of  absorb- 
ing molecules  with  each  other  to  the  effect  or  collisions  between 
absorbing  and  non-absorbing  molecules.  B is  usually  defined  with 
respect  to  nitrogen  as  the  broadening  gas  since  nitrogen  is  the 
major  constituent  of  the  atmosphere,  and  is  not  generally  an  infra- 
red absorber. 

Another  situation  of  interest  in  this  study  is  that  of  triple 
mixtures  of  an  absorbing  gas  a and  two  non-absorbing  gases  b and 
nitrogen.  For  this  situation  Eq.  (57)  beco.  s 


Using  Eqs.  (70),  (71),  (72)  and  (62),  Eq.  (69)  becomes 
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(73) 


1/2 

[PC  + P.  C . + P,,  C 
a aa  b ab  N2  aN2 

Now  B will  be  defined  as  before  as  the  self-broadening  coefficient 
of  a with  respect  to  N2* 

C 

(74)  B = 

^aN2 

For  a broadening  gas  other  than  nitrogen  a foreign  broadening 
coefficient  is  defined 


(75) 


Using  Eqs.  (74)  and  (75),  Eq.  (73)  becomes 


/ 1 

" (sTTrj  ^aN2  ^'^a  ^ '^b  ^b  '^N2^ 


where  here  P is 
e 

(78)  P^  " P B + P.  F.  + P^, 
' ' e a b b N 


Note  that  in  the  above  discussion  Da^-j,  the  sums  of  the 
optical  collision  diameters  was  assumed  to’be  independent  of  tem- 
perature, pressure  and  frequency.  This  leads  to  the  expectation 
that  the  half-width  is  the  same  for  all  lines  in  a vibration-rota- 
tion  band  and  that  the  half-width  depends  on  the  temperature  as 
(1/T)'/^.  In  fact  the  half-width  varies  from  transition  to  transi- 
tion as  do  B and  F-j  and  tomparature  dependence  is  not  always  (1/T)^ 
RicCll]  gives  a good  discussion  of  the  complexity  of  the  determi- 
nation of  the  Lorentz  half-width  and  a discussion  of  other 
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derivations  of  the  Lorentz  half-width.  For  the  purposes  of  the 
calcu'uHiions  used  in  this  investigation  the  above  theory  will  be 
assumed  adequate. 

At  frequencies  removed  from  line  center,  the  Lorentz  line 
shape  does  not  in  general  give  an  accurate  representation  of  the 
true  line  shape.  Long,  et  al[12l  found  that  water  vapor  window 
absorption  between  5 and  6 microns  could  be  more  accurately  pre- 
dicted by  a line  shape  having  more  absorption  in  the  wings  than 
the  Lorentz  profile.  Also  Benedict,  et  al[13]  found  that  for  the 
fundamental  band  of  CO  the  absorption  in  the  wings  of  lines  was 
less  than  predicted  by  the  Lorentz  shape. 

At  pressures  below  about  0.5  torr  the  line  shape  is  dominated 
by  Doppler  broadening,  and  at  pressures  greater  than  about  75  torr 
it  is  dominated  by  collision  broadening.  Between  0.5  torr  and  75 
torr  both  collisior  broadening  and  Doppler  broadening  are  import- 
ant. In  this  pressure  range  a line  shape  called  the  Voigt  profile 
is  used  which  reduces  to  the  Doppler  profile  at  low  pressures  and 
the  Lorentz  profile  at  high  pressures.  The  line  shape  factor  for 
the  Voigt  profile  is  given  by[14] 


(79) 


g(v-v^)  = 


2 


“l 

r e-^ 

dt 

3/2 

IT  ^ 

-00  2 

“d  ^ 

T~ 

“L  + 

where  oq  is  the  Doppler  half-width  as  in  Eq.  (53)  and  a,  is  the 
Lorentz  half-width  as  in  Eq.  (56).  A Fortran  program  for  the 
evaluation  of  Eq.  (79)  has  been  made  available  by  Charles  Young 
[15] 


B.  Continuum  Absorption 

Continuum  absorption  in  the  DF  laser  region  from  3.6  to  4 
microns  arises  from  two  sources.  One  is  water  vapor  absorption 
which  is  thought  to  be  caused  by  the  far  wings  of  the  strong  water 
bands  at  2.7  and  6.3  microns.  The  other  is  pressure  induced  nitro- 
gen absorption.  Burch[4]  has  investigated  both  these  effects  ex- 
perimentally using  a black-body  source  and  a conventional  spectro- 
meter. Burch's  results  have  been  used  in  this  study  to  help  pre- 
dict absorption  at  each  of  the  DF  laser  lines. 

In  Burch's  work  the  absorption  coefficient  has  units  mole- 
cules"*-cm2  and  will  be  called  k'  in  the  following  discussion. 

The  absorber  amount  is  then  expressed  in  molecules-cm2  and  will 
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be  called  u here.  For  a given  sample  u may  be  determined  from  the 
partial  pressures  of  the  absorbing  gas,  the  length  of  the  sample 
and  the  temperature  as  follows 

(80)  u(molecules/cm^)  = 2.69  x 10^^  p(atm)  L(cm)  (273/T) 


where  p is  the  pressure  of  the  absorbing  gas  in  atmospheres,  L is 
the  path  length  in  cm  and  T is  the  temperature  in  degrees  Kelvin. 
The  transmittance  is  then 


(81) 


Assuming  for  the  moment  that  the  Lorentz  line  shape  is  valid 
in  the  far  wings  of  lines,  k'  for  a single  line  is  given  by  Eqs. 
(27)  and  (56) 

Ct 

(82)  k'  2 2 (mol-cm~^)~\ 


For  v-v  >>  a,  this  becomes 
0 L 


(83)  k'  = i!. 

where  is  from  Eq.  (63) 

(8^)  \ = (b^Tkt)  ["a  'aa  0,^] 


Since  a|_  is  proport  i;-“'al  to  pressure,  k'  is  also  proportional  to 
pressure.  It  followb  then  that  for  continuum  absorption  due  to 
far  wings  k'  would  be  expected  to  have  the  form 

(85)  k'  = C^p^  + Pf  (mol-cm’^)'l  . 

In  order  to  obtain  enough  absorption  to  measure  the  water 
continuum  accurately  Burch  was  forced  to  use  higher  temperatures 
so  that  the  water  vapor  pressure  could  be  increased  substantially 
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over  that  normally  occurring  in  the  atmosphere.  He  measured  the 
water  continuum  between  2400  cm"'  and  2800  cm"'  by  tuning  his  spectro- 
meter to  points  in  the  spectrum  which  appeared  to  be  free  of  local 
water  lines  and  measuring  the  transmittance.  He  found  that  the  ab- 
sorption coefficient  k'  was  proportional  to  pressure  as  indicated  by 
Eq.  (85).  This  incidentally  does  not  indicate  that  the  Lorentz  line 
shape  is  valid  in  the  far  wings,  but  rather  that  the  valid  line 
shape  has  the  same  pressure  dependence  as  the  Lorentz  shape. 

Burch  measured  the  continuum  absorption  of  pure  water  samples 
at  three  temperatures  and  found  that  Cg  is  proportional  to  A exp 
(B/T)  where  A and  B vary  with  frequency.  Using  this  relationship 
Burch  deduced  a curve  for  the  water  vapor  continuum  at  296K  from 
the  data  at  higher  temperatures.  His  data  are  reproduced  in  Fig.  3. 

He  measured  Cf^2  higher  temperature  and  found  that 


(86)  ~ = 0.12  + 0.03. 

Unfortunately  0^2  cannot  be  measured  at  normal  atmospheric  pressure 
and  temperature  since  the  absorption  is  too  small  so  it  is  assumed 
that  0^2  ratio  at  296K  as  at  the  higher  tem- 

perature. 

It  is  useful  to  represent  the  extinction  coefficient  in  units 
km"l . From  Eqs.  (80)  and  (85) 


1 Q 

(87)  k'u  = (Cs  Ps(atm)  + p^(atm))  2.69  x 10  Pg(atm)L(cm) 

X (273/T) 

-1 

Rearranging  gives  k in  units  cm 

(88)  k = ^ = 7.34  X 10^^  ^ (C3  P3  + p^)  cm"^ 

The  extinction  coefficient  k in  km  ^ is  then 

(89)  k = 7.34  X (C3P3  + p^)  km"^ 

For  the  water  continuum  data  at  296K,  k can  be  written 
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Spectral  pi 
(Burch) . 


(90) 


k = 2.48  X 10^^  Ps^^S  '^S  * 


or 

(91)  k = ^''*V33°~  PsS  '*•”  f'S  ^ Pn>  '^'"■'  • 

The  total  pressure  Pj  is  plus  P,^  so  Eq.  (91)  may  be  rewritten 

(92)  k=  2.98  X 10^^  C^p^  (7.33  Pg  + Pj)km"‘. 

Figure  4 was  derived  by  evaluating  Eq.  (92)  at  14.26  torr  H2O  and 
one  atmosphere  total  pressure  with  Cg  taken  from  Fig.  3. 

The  nitrogen  continuum  absorption  arises  from  a pressure- 
induced  band  near  2200  cm"'.  A number  of  workers  have,  investigated 
this  absorption  between  2400  cm"'  and  2640  cm"'*  Burch[4]  has 
determined  from  his  measurements  and  those  of  others  that  (-InT) 
is  proportional  to  p^L  just  as  for  the  pure  water  vapor  continuum. 
Figure  5 is  adapted  from  Burch  and  shows  the  nitrogen  continuum 
between  2400  cn:“'  and  2650  cm-1  at  296K. 


C.  Description  of  Calculation  Programs 

The  basic  quantity  calculated  by  the  programs  is 


(93)  (-J-n  T)  = k'u 

2 

where  k‘  is  in  units  (mol  - cm"^)"^  and  u ’s  in  units  mol  - cm  . 
There  are  two  different  equations  used  to  calculate  k'  depending  on 
the  total  pressure.  For  pressures  above  about  75  torr  the  Lorentz 
profile  is  used: 


(94) 


k'  = - 


S “L 


“2  -1 
(mol -cm  ) , 


For  lower  pressures  the  Voigt  profile  is  used. 

.2 


(95) 
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Nitrogen  absorption  from  Burch,  T=296K, 
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In  the  above  equations  is  given  by 

(96)  a.  = 3.5812  X 10'^  v 

u Mo 

and  is  given  by 

“L'law)"'  CaN/e 

where  is  the  effective  broadening  pressure  given  in  general  by 


(98)  P,  = B.ljP,  F, 


Here  P is  the  partial  pressure  of  the  absorbing  gas,  B is  the  self- 
broadening  coefficient,  P-j  is  the  partial  pressure  of  broadening  gas 
1 and  F-j  is  the  foreign  broadening  coefficient  for  broaden* nq  qas  i. 
Note  that  F for  nitrogen  is  defined  to  be  1 . 

For  each  absorption  line  the  AFCRL  line  compilation  gives  the 
line  frequency  in  cm-1 , the  line  strength  in  (mol-cm-2)-l-cm-l  at 
296K,  the  Lorentz  half-width  in  cm~*-atm“l  at  296K,  the  energy  of 
the  lower  energy  state  of  the  transition  in  cm“  , the  energy  levels 
involved  in  the  transition,  and  the  isotope  and  type  of  molecule. 

For  use  at  temperatures  other  than  296K  and  pressures  different 
from  one  atmosphere,  the  line  strength  and  half-width  must  be  cor- 
rected using  the  following  equations  derived  earlier. 


The  programs  described  here  have  been  used  exclusively  for 

temperatures  between  195K 
and  400K.  For  those  conditions  the  term  in  brackets  in  Eq.  (99)  is 
very  nearly  one  and  may  be  safely  ignored.  Also  the  ratio  of 
vibrational  partition  functions  Qy°/Qv  has  been  taken  to  be  nearly 
one  and  therefore  ignored  although  it  could  be  included  if  uncer- 
tainties in  line-shape  and  half-width  were  reduced  to  the  point 
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where  errors  in  Qy°/Qv  were  significant.  Assuming  296K  for  T , 
1 atm  for  P and  using  k = .6951  cm"V°K  the  strength  and  half- 
width  correStions  become 


The  assumed  values  for  the  self-broadening  coefficient  B, 
BX,  and  CX  a'^e  given  for  each  molecule  in  Table  2. 


TABLE  2 

VALUES  USED  FOR  B,  BX,  AND  CX 
IN  ABSORPTION  CALCULATIONS 


B 

BX 

CX 

N20 

1.24 

1.0 

.5 

CH4 

1.3 

1 .5 

.5 

H20 

5.0 

1 .5 

.62 

COp 

1.3 

1 .0 

.58 

°3 

1.0 

1 .5 

.5 

CO 

1 .02 

1 .0 

.5 

Equation  (94)  gives  the  absorption  coefficient  k'  for  one 
absorption  line.  At  any  frequency  there  might  be  several  lines 
contributing  to  the  absorption,  k'  would  then  be  given  by  the 
summation  of  11  the  individual  line  contributions. 

S . “l  . 

I — (mol-cm*^)’^  . 

i (V-V^)  + 

The  first  program  to  be  described  calculates  and  plots  k in 
km"l  versus  wavenumber.  Note  that  k in  km"^  is  k'u  where  u is 
evaluated  using  Eq.  (80)  with  the  path  length  L being  one  kilometer 
Input  parameters  are  the  number  of  absorbers. to  be  considered;  the 
number  of  plots  to  be  made;  the  beginning  wavenumber  of  the  first 


(103)  k'  = - 

TT 
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plot;  the  number  of  wavenumbers  per  plot;  the  total  pressure  in  torr, 
the  temperature;  the  I.D.  number,  desired  isotope,  partial  pressure 
in  torr,  and  broadening  coefficient  for  each  absorber  type;  and 
descriptive  information  which  is  written  on  the  plots.  The  user  may 
specify  a continuum  of  the  form  k (v)  = + A,v  + Ap  v which  is 

added  to  the  calculated  k(v)  at  each  point  on  the  plot.  The  user 
may  also  specify  a set  of  frequencies  corresponding  to  laser  lines 
or  just  frequencies  of  special  interest,  and  the  program  will  draw 
a vertical  line  on  the  plot  at  each  specified  frequency. 

Figure  6 shows  an  example  of  the  interactive  (teletype)  print- 
out which  occurs  during  execution  of  tne  program  with  the  informa- 
tion provided  by  the  operator  underlined.  The  program  then  calcu- 
lates u for  each  absorber  from  the  partial  pressure  and  temperature 
using  Eq.  (78)  and  Pg  for  each  absorber  from  the  partial  pressure, 
total  pressure  and  self-broadening  coefficient.  Data  is  read  from 
the  AFCRL  tape  as  required  and  the  strength  and  half-width  are  cor- 
rected for  temperature  and  pressure  using  Eqs.  (101)  and  (102). 

There  is  a parameter  SLOW  in  the  program  which  may  be  set  to  ignore 
all  lines  with  strength  less  than  SLOW  and  thus  speed  the  calculation. 

The  heart  of  the  program  is  the  subroutine  which  calculates 
k'u.  It  is  adapted  from  one  written  by  Deutschman  and  Calfee[10J. 

At  any  point  v the  subroutine  considers  contributions  from  absorp- 
tion lines  within  BOUND  of  v in  either  direction.  Recognizing  the 
uncertainty  in  the  Lorentz  line  shaoe  at  frequencies  far  from  line 
center,  BOUND  is  usually  set  at  20  cm-'  or  25  cm~l . For  each  separate 
absorption  line  in  the  range  v ± BOUND  the  subroutine  calculates 
the  contribution  using  either  Eq.  (94)  (Lorentz  profile)  or  (95) 

(Voigt  profile)  and  the  appropriate  u for  the  absorber  being  con- 
sidered. Whether  the  Lorentz  or  Voigt  profile  is  used  is  deter- 
mined by  the  ratio  ai/ag  for  each  particular  absorption  line.  If 
ai/aD  is  greater  than  five  the  Lorentz  profile  is  used,  otherwise 
the  Voigt  profile  is  used.  The  contributions  from  each  absorption 
line  are  then  added  together  to  get  the  total  k'u  for  local  line 
absorption.  If  the  user  has  specified  a continuum,  it  is  evalu- 
ated at  this  time  at  the  frequency  v and  added  to  k'u.  This  pro- 
cedure is  repeated  for  each  point  on  the  plot. 

The  program  plots  the  absorption  coefficient  in  km"i  on  a 
logarithmic  scale  with  the  scale  limits  calculated  automatically 
for  the  most  effective  presentation. 

Another  version  of  this  program  plots  the  transmittance  rather 
than  the  absorption  coefficient.  It  would  seem  however  that  this  is 
less  desirable  since  by  improper  choice  of  path  length  and  absorber 
amount  it  would  be  possible  to  calculate  a spectrum  showing  total 
absorption  or  total  transmittance,  whereas  the  absorption  coeffici- 
ent plot  will  always  show  the  spectral  structure  no  matter  what 
absorber  amount  is  specified. 


¥■ 
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Copy  avdlabJe  to  DDC  do5s  not 
permit  fully  legible  reproduction 


CATr.  ?8/ 17/75  TIME  19!5Aie6 

ENTER  THE  NUMBER  OF  ABSORBERS  ANC  THE  NUMBER  OF  PLOTSl,  1 

ENTER  THE  BEGINNING  VAVHJU’MEE? 

ANL  THE  VAVENIMBER/PLOT.  S65?,  \? 

ENTER  THE  TOTAL  PRESSURE  IN  TORP.  7^ 

ENTER  THE  TEMPERATURE  IN  DEGREES  F.69.8 

ENTER  ID  NIMEEP,  PARTIAL  PRESSURE,  AND  BROADENING 
GOEFFICIENT  FOR  EACH  ABSORBER. 

ENTER  ID  Nl'MBER  OF  ABSORBER  NO.  1 
:»H20  2=G02  3=03  A=N20  5=  CO  6=CHA  7=02 


ENTER  NIMEER  OF  DESIRED  ISOTOPE. 

IF  ALL  ISOTOPES  APE  DESIRED  ENTER  P.  1 62 

ENTER  THE  ABSORBER  AMOUNT  IN  TORR.  1A.26 

ENTER  THE  SELF- BROADENING  COEFFICIENT.^ 

DO  YOU  NISH  A CONTINUU'M?  NO 

ENTER  AESOPEER.  DESCRIPTION  INFORMATION  ON  T'.'O  :MES 
OF  UP  TO  A2  CHARACTERS  EACH.  THE  FIRST  LINE  'TLL 
APPEAR  IMMEDIATELY  TO  THE  RIGHT  OF  ABSORBERS:  , 

THE  SECOND  LINE  VILL  APREAP  TO  THE  RIGHT  OF 
AMOIINTC  TORR) ; AND  IMMEDIATELY  BELOV  THE  FIRST  LINE. 

HDO 

1A.26 

VHERE  IS  THE  ABSORPTION  LINE  DATA? 

ENTER.  .MTe  OR  .MTl  .MTl 

EvlTER  TAPE  START  POINT:  FILE  NIMBER  AND  RECORD  NIMBERA,  8 50 

VHICH  TAPE  l^NIT  VILL  PLOT  DATA  BE  ’(Til  TTEN  ON? 

ENTER  .MT0  OF  .MTl  .MT0 

ENTER  NU1EER  OF  -FIRST  FILE  TO  BE  '.T1ITTEN.216 

VHERE  IS  THE  LASER  LINE  DATA? 

ENTER  FILE  AND  USER  NAME  ON  TVO  LINES. 

NDFLN 
3271 A 

DO  YOU  VI SH  TO  RUN  THIS  PROGRAM  IN  EACHGROU’ND?  YES 


Fig,  6.  Example  run  of  spectra  plotting  program. 


31 


The  program  which  calculates  single  frequency  absorption  is 
very  similar  with  a few  exceptions.  The  absorption  coefficient  is 
calculated  for  only  one  absorber  at  a time,  however  calculations 
may  be  made  for  a number  of  absorber  pressures  and  frequencies  at 
the  same  time. 


The  other  main  difference  is  that  the  Lorentz  line  shape  is 
not  used,  but  rather  a modification  of  the  Lorentz  line  shape  sug- 
gested by  Trusty[12].  The  alternate  shaoe  has  the  form 


(104) 


k = 


C(v,n)S  \ 
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C(»m.n)  = 
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k dv 
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is  added  so  that  the  line  strength  has  the  usual  definition  of  the 
integrated  absorption  coefficient  over  the  entire  line.  The  inte- 
grations in  Eq.  (105)  may  be  performed  directly  to  give  for  C(v  ,n) 


(106) 


2 [Tan‘^(v^)  + Vn,/((^n,^+l ) (n-D)] 


The  modification  location  v and  modifier  power  n are  specified  by 
the  user.  For  ^ 30<  (X|_  and  n = 2 the  modified  shape  reduces  to 
the  Lorentz  shape.  For  r\  < 2 the  alternate  line  shape  has  more 
wing  absorption  than  the  Lorentz  shape  and  for  n > 2 the  wing  ab- 
sorbtion  is  less. 

Figure  7 shows  an  example  of  the  interactive  (teletype)  print- 
out which  occurs  during  execution  of  the  program  with  the  informa- 
tion provided  by  the  operator  underlined. 

Appendix  B gives  listings  of  the  programs  and  a somewhat  more 
detailed  discussion  of  the  mechanics  of  the  programs, 
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prtss^  AES^  TLWP,  MALFA,  LTA>  LIlMtS,  EPOAE^  DATA,  E'<E^  llltj  Sic,  TTOPE, 
'ZAP. 


LATE  08/17/7F  THIE  19:A8:22 

H4TE?  THE  TOTAL  P^LSSI'PE  I’J  TOPP.  7*g 

EMTEP  THE  TEMP  I'l  LEG  F.  6'?. 8 

EMTEP  IL  M3.  OF  AEEOPBEP 
l = H20  2=C02  3=02  A=JEO  £=  CO  6=CHA  7=02 


E.MTEP  MlMBEP  OF  EESIPEE  IE0T3FE. 

IF  ALL  ISOTOPES  APE  LEEir'EL  LMTEP  0.  1 62 

HOV  MA:JY  ABSOPEEF  PFESSVPES  VILL  EE  VEEL?J_ 

EMTEP.  THEM  OM  OME  LIM£>  FPEE  FORMAT. 

1 A.  26 

E.MTEP  THE  SEL  F- EPO  ALEMIMG  COEFFI  Cl  EMT.  ^ 

EMTEP  MO  LI  FI  CATION  LOCATION  IM  HALFVIDTHS.  37 
EMTEP.  MOCIFIEP  POVEP.  £ 

VILL  CALCINATION  FPEPVEMCIES  DE  PEAL  FPOM  A ■’FILE"  OP  THE  "TTY"?  riLc 

EMTEP  FILE  AML  USEP  NAME  ON  T''0  LIMES. 

MLFLN 

3271A 

. '■•HEPE  IS  THF  LINE  LATA? 

EMTEP  .MT7  .MTl 
.MT  1 

I EMTEP  TAPE  CTAPT  POINT!  FILE  Ml?4BEP.  AML  PECOPD  Me.  A,  7S3 

CO  YOl'  VI SH  TO  PPM  OFF-LINE?  YES 


Fig.  7.  Example  run  of  fixed-frequency  program. 
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D. 


Calculations  in  the  DF  Laser  Region 


The  programs  described  in  the  last  section  and  the  AFCRL  line 
compilation  were  used  to  calculate  both  single  frequency  absorption 
and  synthetic  spectra  for  each  laser  line  available  from  the  probe 
laser  used  in  this  study. 

Calculations  were  made  for  the  AFCRL  mid-l<j  !:itude,  sea-level, 
summer  model [icl-  The  molecules  considered  in  the  calculation  are 
H2O,  N2O,  CH4,  and  CO2.  The  contributions  from  pressure-induced  N2 
absorption  and  water  continuum  absorption  are  obtained  from  Fig.  5 
and  Fig.  4 respectively.  Parameters  for  the  mid-latitude  summer 
model  are  given  in  Table  3. 

TABLE  3 

AFCRL  MID-LATITUDE  SUMMER  MODEL 


T 

= 294‘’K 

P 

= 760  torr 

P(H20) 

= 14.26 

torr 

P(N20) 

= 2.13  ) 

< 10“ 

P(CH4) 

= 1.216 

X 10 

P(C02) 

= .251 

torr 

torr  (.28  ppm) 

^ torr  (1 .6  ppm) 


The  HDO  line  strengths  on  the  AFCRL  tape  are  scaled  assuming 
an  isotopic  abundance  of  HDO  relative  to  total  water  of  0.03%. 

SLOW  was  set  at  10"27  for  these/Calculations  so  that  essentially 
all  the  absorption  lines  are  included.  BOUND  was  20  cm"'.  B,  BX, 
and  CX  used  for  each  molecule  are  those  given  in  Table  2.  The  DF 
laser  frequencies  used  are  those  measured  by  Heath,  et  al[17]  and 
are  presumed  accurate  to  ± .003  cm”  • 

The  calculations  at  the  individual  laser  frequencies  are 
listed  in  Table  4.  The  calculations  for  each  molecule  are  listed 
separately  and  then  added  together  for  the  total  predicted  absorp- 
tion. For  the  water  contribution,  the  local  HDO,  local  H2O  (ex- 
cluding HDO),  and  H2O  continuum  are  tabulated  separately  so  that 
the  relative  importance  of  each  factor  can  be  noted. 

The  synthetic  spectra  are  presented  in  Figs.  8 through  35. 

They  are  vertical  lines  drawn  on  the  plots  to  indicate  laser  line 
positions.  The  laser  lines  are  identified  by  transition,  and  the 
major  absorption  lines  near  each  laser  line  are  identified.  The 
water  continuum  and  nitrogen  absorption  have  not  been  included  since 
it  is  slowly  changing  with  respect  to  frequency  and  the  effect  would 
be  to  obscure  the  spectral  structure. 
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Fig.  10,  Calculated  spectrum  near  the  3-2  P(ll) 
DF  laser  line. 


laser  line. 


Fig.  15.  Calculated  spectrum  near  the  3-2  P(7) 
DF  laser  line. 
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WAVENUMBER 

Fig.  20.  Calculated  spectrum  near  the  2-1  P(8) 
DF  laser  line. 


Calculated  spectrum  near  the  1-0  P(7) 
and  3-2  P(3)  DF  laser  lines, 
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HRVENUMBER 

Fig,  25,  Calculated  spectrum  near  the  1-0  P(9) 
OF  laser  line. 
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Calculated  spectrum  near  tne  2-1  P(4) 
DF  laser  line. 
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WAVENUMBER 

Fig.  29.  Calculated  spectrum  near  the  1-0  P(7) 
DF  laser  line. 
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Calculated  spectrum  near  the  2-1  P(3) 
DF  laser  line. 
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CHAPTER  III 

EXPERIMENTAL  EQUIPMENT 


A very  important  part  of  this  study  was  developing  the  experi- 
mental system  necessary  for  making  measurements  of  the  desired  high 
accuracy.  The  experimental  apparatus  making  up  this  system  can  be 
divided  into  four  groups: 

1 . Single  line  DF  laser 

2.  Auxiliary  optics 

3.  Detectors  and  electronics 

4.  Multiple  traversal  absorption  cell 

The  equipment  in  these  groups  will  be  described  in  this 
chapter. 


A.  Design  and  Construction 

of  the  DF  Laser 

1 .  Introduction 

A probe  laser  for  atmospheric  absorption  measurements  should 
have  the  following  characteristics: 

1.  It  should  operate  on  only  a single  line  and  be  easily 
tunable  from  line  to  line. 

2.  It  should  be  fairly  easy  to  use  , that  is  it  should  not 
require  a full-time  technician  to  operate. 

3.  It  should  be  economical  to  operate. 

4.  Amplitude  stability  should  be  as  good  as  possible.  ^ 

At  the  time  this  study  was  begun  and  indeed  even  now  such  a 
DF  laser  is  not  commercially  available.  It  was  therefore  necessary 
to  design  and  construct  the  laser. 

The  DF  laser  is  a chemical  laser,  that  is  the  upper  laser 
levels  are  populated  by  the  release  of  energy  from  a chemical  re- 
action. The  lasers  used  in  this  study  are  similar  but  not  identical 
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to  one  described  by  UUee[l8].  A mixture  of  helium,  sulfur  hexa- 
fluoride (SFc),  deuterium,  and  oxygen  gases  is  introduced  into  the 
discharge  tube.  A 12  to  15  kv  electrical  pulse  of  less  than  one 
microsecond  duration  dissociates  the  SFs  p-oducing  free  fluorine 
which  reacts  with  deuterium  to  form  excited  DF  which  then  emits  a 
laser  pulse  about  one  microsecond  long.  Helium  is  added  to  serve 
as  a heat  sink  and  promote  discharge  stability.  A small  amount  of 
oxygen  is  adied  to  keep  sulfur  deposits  from  forming  on  the  walls 
of  the  discharge  tube.  Note  that  all  the  input  gases  are  non-cor- 
rosive and  non-poisonous  and  require  no  speciail  handlirig, 


2.  Original  DF  laser 

a.  Gases 

The  gases  used  were  helium,  SF5,  deuterium,  and  oxygen  as 
described  in  part  1.  A conventional  15  cfm  mechanical  pump  was 
used  to  maintain  a pressure  of  3 to  12  terr  in  the  discharge  tube 
and  a flow  rate  sufficient  to  change  the  gas  mixture  between  pul  res 
It  was  found  that  the  optimum  pressure  and  gas  mix  varied  according 
to  the  particular  line  which  was  oscillating, 

Estimated  mass  flow  rates  for  operation  at  6 torr  are:  helium 
21  gm/hour;  SF6  .49  gm/hour;  deuterium,  1.4  gm/hour.  Estimated 
operating  cost  with  commercially  supplied  deuterium  is  $4. 00/hour. 


b.  Laser  tube  geometry 

The  laser  tube  geometry  is  shown  schematically  in  Fig.  36. 

The  tube  was  constructed  of  10  mm  I.D.  pyrex  tubing  with  two  inlet 
ports  and  three  output  ports.  The  electrodes  were  Kovar  glass  to 
metal  seal.  Swagelok  fittings  were  used  to  make  connections  at  the 
inlet  and  outlet  ports. 

The  windows  were  calcium  ^luoride  with  the  Brewster  angle  in 
the  proper  orientation  for  operation  with  a grating. 


c.  Optical  cavity 

The  optical  cavity  consisted  of  a 300  line/nm  grating  blazed 
at  3 pm  and  a 20  m radius  of  curvature  germanium  mirror  coated  for 
greater  than  80%  reflectivity  between  3 pm  and  4 pm.  The  grating 
and  the  mirror  were  separated  by  121  cm. 

The  grating  mount  was  designed  by  Professor  Damon  at  Ohio 
State  so  that  the  grating  could  be  easily  aligned,  and  so  that  the 
grating  could  be  easily  and  reproducibly  tuned  to  a specific  laser 
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Fig,  36,  Schematic  diagram  of  first  DF  laser  tube. 


line.  A photograph  of  the  grating  mount  is  shown  in  Fig.  37.  This 
mount  proved  to  be  quite  satisfactory  and  was  used  in  all  versions 
of  the  laser. 

The  output  mirror  was  mounted  in  a piezoelectric  drive  : 
which  was  part  of  a commercial  qimbal  mount  which  was  mounted 
on  a horizontal  translation  stage.  The  mirror  could  then  be  moved 
to  adjust  the  longitudinal  cavity  mode  for  optimum  operation  on  a 
single  line.  It  was  found  however  that  longitudinal  adjustment  had 
no  effect  on  the  laser  output,  probably  because  of  the  extremely 
high  gain  for  pulsed  operation. 

The  grating  mount,  laser  tube  and  mirror  mount  were  all  attach- 
ed to  a limBstone  slab  five  feet  long,  one  foot  wide,  arid  three  inches 
thick.  This  arrangement  proved  to  be  quite  stable.  Once  the  optical 
cavity  was  aligned  it  maintained  alignment  indefinitely  with  only 
minor  adjustments  required  when  the  laser  tube  was  changed. 
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d.  Power  supply 

The  laser  reaction  was  initiated  bv  disrharninn  a n no  re 
capacitor  through,^  thyratron.  The  capacitor  wa  ^cKLoeSlo  2 to 
16  ky  ujing  a commercial  high  voltage  power  supply.  Thfpulser 
supplied  by  Aerospace  Corporation  and  was  capable  of  singfe  0015^ 
operation  or  repetitive  pulsing  up  to  about  100  pu^Lrjer  second. 

necessary  to  repackage  the  Aerospace  pulser  to  reduce 
electrical  noise  problems.  This  is  discussed  in  mSre  ^tail  ^ 


-•  Laser  aliqhment 

n +-n  alignment  of  the  laser  was  relatively  simple  althouah  tne 
ouput  mirror  was  opaque  to  visible  light.  A helium-neon  (HeNe) 

coll  inear  with  the  laser  tube  with  the 

ti.hP  ^ ^ mounted  at  either  end  of  the  laser 

tube.  Then  a flat  mirror  was  attached  to  the  grating  table  and  the 
grating  mount  was  adjusted  so  the  grating  rotation  axis  wa.s  exactlv 

reflected^the  HeN^^i  occurred  when  the  flat  mirror 

retiected  the  HeNe  laser  beam  exactly  back  on  it<;pif  thic  f'miivi 

determined  because  the  reflected  beam  fed  back  into  the  HeNe 

amplitude  to  fluctuate  noticably.  Then 
the  flat  mirror  was  removed  and  the  grating  was  installed  The 
grooves  of  the  grating  were  then  aligned  paraifel  wi?h  the  al?s  rf 

orders  of  the  HtNe  laser  beam  were  reflected  st^'aiaht  hart  tr.  the 
laser  as  the  grating  angle  was  tuned.  Thirinsured  tha^L  gra?i^ 
was  aligned  for  the  wavelength  rang''  from  2.5'  ym  to  4.5  urn.  The  ^ 
output  mirror  and  mirror  mount  were  now  installed  so  that  the  outnut 
mirror  was  centered  on  the  HeNe  laser  beam.  Then  the  output  4r™f 

hark*^f“th®'^  '■'’I  cf  the  HeNe  laser  beam  from  the 

back  of  the  mirror  went  straight  back  to  the  laser  The  nii+nn+  miv 

ror  could  be  adjusted  closely  enough  this  way  so  that  with  the  lasei- 

0 c^'^?tio^'  anglelh^^e  wouto  bl 

oscillation.  Fine  adjustment  of  outpuL  mirror  orientation  could 

then  be  made  to  optimize  the  laser  output.  The  gracing  mJcromeSer 

tuning  calibration  could  be  determined  from  the  micrometer  nosit  on«; 

fleeted.  If  a least  square  fit  of  a quadratic  is  then  nade  to  these 
four  points,  the  micrometer  position  for  any  given  lase-  line  mav 
thcd  be  determined  from  the  fitted  quadratic.^  This  calibration  me- 
thod was  tested  by  observing  the  laser  output  with  a 1 meter  Czerny- 
Turner  spectrometer  and  found  to  be  valid.  This  -’s  hardly  suroris'ir.a 

Si-S^slUs!'''^  - ^^--lent  to  a 1.21  Lter  sJeaV^L^^S 
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f . Laser  performance 

The  laser  was  operated  and  found  to  oscillate  on  25  lines  with 
average  power  at  20  pulses  per  second  of  from  .04  to  1 milliwatts 
per  line.  A drawing  of  an  oscilloscope  trace  of  the  pulse  is  shown 
in  Fig.  38.  The  laser  could  be  operated  up  to  about  35  pulses  per 
second  before  the  laser  output  per  pulse  started  to  deteriorate. 

The  pulse  rate  was  probably  limited  by  the  gas  flow  rate.  There  | 

were  however,  two  serious  problems.  | 

One  was  the  great  amount  of  RF  electrical  noise  generated  by  i 

the  high-voltage  pulser  interfered  with  the  operation  of  the  elec-  | 

tronics  associated  with  the  absorption  measurement  experiment.  An 
attempt  was  made  to  reduce  the  noise  by  carefully  repackaging  the  ’ 

pulser.  An  improvement  was  observed,  but  the  noise  was  still  in-  » 

tolerable.  The  entire  laser  including  the  power  supply  was  then 
placed  in  a shielded  room  with  a double  thickness  of  copper  screen 
on  all  sides.  The  power  to  the  shielded  room  was  connected  through 
line  filters,  and  the  room  was  separately  grounded  to  an  eight  foot 
ground  rod  driven  through  the  floor.  The  gas  line  and  pumping  lines 
came  through  brass  or  copper  pipes.  A one  inch  hole  was  cut  in  the 
wall  to  let  the  laser  beam  out.  The  pulse  trigger  signal  was  coupled 
through  the  wall  using  an  infrared  light  emitting  diode  and  a photo- 
transistor. These  rather  drastic  measures  completely  eliminated  the 
RF  noise  problem. 
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The  second  problem  which  occured  was  sulfur  depositing  on 
Brewster  windows  and  on  the  inside  of  the  laser  tube.  The  deposit 
rate  was  such  that  the  windows  had  to  be  cleaned  after  about  one 
hour  of  operation.  An  attempt  was  made  to  reduce  the  sulfur  deposit 
rate  by  adding  more  oxygen  to  the  gas  mixture.  However  this  '-ad 
the  effect  of  drastically  reducing  the  laser  power. 

3.  Second  DF  laser 

To  solve  the  problem  of  window  contamination  a new  tube  was 
designed  similar  to  the  first  but  with  ports  added  near  each  Brew- 
ster window.  A schematic  of  the  new  tube  is  shown  in  Fig.  39.  The 
idea  is  to  admit  helium  through  these  additional  ports  to  keep  re- 
action products  out  of  the  area  near  the  windows. 

Other  modifications  were  also  made  in  the  tube  at  this  time. 
The  tube  I.D.  was  reduced  to  about  8.5  mm  to  improve  discharge  sta- 
bility. The  electrode  configuration  was  changed  from  Kovar  glass 
to  m^tal  seals  to  1/8  inch  copper  tubing  inserted  as  close  to  the 
active  region  as  possible  without  obstructing  the  optical  path. 

This  was  done  to  reduce  the  voltage  and  energy  required  to  initiate 
the  chemical  reaction. 
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Fig.  39.  Schematic  diagram  of  second  DF  laser  tube, 
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The  new  tube  worked  quite  satisfactori'iy.  The  added  helium 
flow  at  the  windows  not  only  had  the  expected  effect  of  eliminating 
window  contamination,  but  also  the  unexpected  effect  of  reducina 
sulfur  deposits  in  the  tube.  Thus  the  oxygen  flow  »,ould  be  safely 
reduced  to  a point  where  the  laser  output  was  no  longer  affected. 
The  helium  flow  rate  past  the  windows  was  such  that  with  the  main 
gas  mixture  shut  off  the  tube  pressure  was  about  5 torr. 

Both  the  output  power  and  the  pulse  to  pulse  stability  were 
somewhat  improved  with  the  new  tube.  The  laser  also  operated  on 
30  lines  instead  of  25  with  the  added  lines  being  primarily  the 
lower  rotational  transitions  of  the  1-0  band'.  This  may  be  because 
the  helium  flow  at  the  windows  removed  unexcited  DF  from  the  opti- 
cal path. 

This  second  version  of  the  laser  worked  reasonably  well  and 
most  of  the  measurements  described  in  the  next  chapter  were  made 
using  this  laser. 


4.  The  portable  DF  laser 

After  operating  the  laser  described  above  for  a considerable 
time  it  was  decided  that  a smaller  more  portable  laser  was  needed. 

In  particular  a laser  was  needed  which  had  very  little  RF  noise 
generation,  since  the  large  shielded  room  was  taking  up  much  needed 
laboratory  space.  It  was  also  felt  that  improved  optical  stability 
could  be  obtained  if  the  laser  could  be  mounted  on  the  same  table 
v.ith  the  external  optics. 

The  pulsed  laser  built  by  UlteePS]  used  a different  power 
supply  design  than  that  used  in  the  Aerospace  pulser.  IHtee's  power 
supply  charged  a 1 microfarad  capacitor  to  about  600  voits  and  dis- 
charged it  through  a pulse  transformer  to  get  the  high  voltage  re- 
quired by  the  laser  tube.  This  design  has  the  advantage  that  the 
high-voltage  pulses  are  restricted  to  the  output  of^ the  pulse  trans- 
former and  the  laser  tube.  These  components  can  be  fairly  easily 
shielded  to  eliminate  RF  noise  emission.  Using  circuit  diagrams 
provided  by  Ulteep9],  a new  pulser  was  constructed. 

Also  at  this  time  a new  tube  was  designed  and  built  which 
had  an  active  length  of  20  centimeters  and  inside  diameter  of  5.5 
millimeters.  This  compares  to  an  active  length  of  about  30  centi- 
meters and  inside  diameter  8.5  millimeters  for  the  previous  tube. 

The  new  tube  had  one  inlet  at  the  center  of  the  tube  for  the  main 
gas  mixture,  and  an  exhaust  port  to  either  end  as  well  as  the  helium 
inlet  ports  at  the  Brewster  windows. 

It  had  been  observed  on  the  previous  laser  that  the  copper 
electrodes  became  fouled  after  a period  of  time,  particularly  the 
one  at  the  high  voltage  end.  This  could  possibly  have  been  copper 
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sulfate.  A couple  of  different  electrode  configurations  were  tried 
to  eliminate  the  fouling.  The  most  satisfactory  solution  seems  to 
be  a nickel  rod  at  the  low  voltage  end  and  a tungsten  rod  at  the 
high  voltage  end  insulated  with  Teflon  tape  except  aL  the  tip. 

The  optical  cavity  for  the  new  laser  consisted  of  a 300  line/ 
millimeter  grating  blazed  at  3.5  ym  and  the  same  20  meter  radius  of 
curvature  germanium  mirror  used  in  the  previous  laser.  The  grating 
and  end  mirror  were  separated  by  80  cm.  The  grating  mount  was  of 
t‘he  same  design  as  that  used  on  the  previous  laser.  The  output 
mirror  was  mounted  in  a commercial  gimbal  mount  as  before,  but  . 
there  was  no  provision  made  for  tuning  the  cavity  length  since  it 
had  been  found  to  have  little  effect  in  the  earlier  laser. 

The  laser  was  constructed  on  a 3/4-inch  aluminum  plate,  ten 
inches  wide  and  thirty-nine  inches  long.  Along  one  side  of  the 
laser  was  a 1/4-inch  aluminum  plate.  Attached  to  the  plate  were  a 
box  containing  the  oscillator,  feedthroughs  for  the  gases  and  pump- 
ing port,  and  a vacuum  gauge  which  reads  \.he  pressure  at  the  exhaust 
ports  of  the  laser  tube.  The  mirror  mount,  grating  mount,  laser 
tube,  storage  capacitor,  thyratron,  and  pulse  transformer  were  all 
mounted  on  the  base  plate.  There  was  an  aluminum  box  attached  to 
the  base  plate  and  side  plate  which  eiiclosed  everything  but  the 
mirror  mount  and  grating  mount.  There  were  small  holes  in  the  ends 
of  the  box  to  let  the  laser  beam  out.  Figure  40  shows  the  laser 
with  the  shield  box  removed. 

Since  the  laser  tube  was  not  exactly  the  same  as  'ITtee's,  it 
, was  necessary  to  use  a different  pulse  transformer  and  storage  capa- 
' citor.  The  DC  voltage  to  charge  the  capacitor  is  supplied  by  Power 
Designs  Pacific,  Inc.,  Model  HV-1547  photomultiplier  tube  supply. 

The  capacitor  is  charged  to  about  1200  to  1300  volts.  Figure  41 
shows  the  electrical  schematic  of  the  Ultee  power  supply  as  modi- 
fied. Ultee  used  an  t.G/iG.  Model  TS-185  pulse  transformer  rather 
than  an  E.G&G.  Model  TS-146A  and  a one  microfrad  capacitor  instead 
of  a .333  yfd  capacitor.  When  the  laser  is  operating  at  50  pulses 
per  second  the  photomultiplier  supply  provides  about  1600  volts 
at  about  20'ma.  The  capacitor  is  charged  to  about  1250  volts  since 
the  charging  time  between  pulses  is  1.5  time  constants. 

The  new  laser  completely  eliminated  the  RF  noise  emission 
experienced  with  the  earlier  laser,  and  in  addition  the  average 
power  in  the  new  laser  was  increased  by  a factor  of  four  or  five. 

Part  of  the.  increase  comes  from  a faster  pulse  rate  (50  pps  rather 
than  35  pps)  made  possible  because  the  tube  was  smaller  and  the 
vacuum  pump  was  the  same.  Figure  42  shows  a photograph  of  an  oscil- 
loscope trace  of  the  pulse. 

A laser  very  similar  to  that  just  described  was  built  and 
loaned  to  Professor  Rao  of  the  Ohio  State  University  Physics  depart- 
ment so  the  laser  line  frequencies  could  be  accurately  determined 
using  a very  accurate  grating  spectrometer.  Using  the  laser  emission 
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Fig.  41.  Circuit  diagram  of  the  power  supply  for  the 
portable,  pulsed  DF  laser. 


spectra  in  conjunction  with  DF  absorption  spectra,  the  values  for 
the  laser  frequencies  were  determined  to  an  accuracy  of  ± .003  cm-1. 
These  frequency  values  are  given  in  Table  5 for  each  of  the  lines 
observed  in  the  pulsed  lasers  described  here. 


B.  Optics 

The  next  important  aspect  of  the  experimental  system  is  the  op- 
tics necessary  to  focus  the  laser  beam  into  the  absorption  cell  and 
onto  the  detectors.  There  were  two  optical  setups  used,  one  with 
the  DF  laser  mounted  separately  in  a shielded  room,  and  the  other 
with  the  DF  laser  mounted  on  the  main  optical  table  with  the  White 
cell  entrance  optics. 

The  first  optical  layout  to  be  described  is  the  one  used  with  the 
DF  laser  mounted  separately  in  the  screen  room.  This  is  the  arrange- 
ment which  was  used  in  making  most  of  the  measurements  described  in 
the  next  chapter.  The  layout  is  shown  schematically  in  Fig.  43. 

The  DF  laser  beam  comes  from  the  shielded  room  and  is  directed 
onto  the  main  optical  table  by  mirror  Ml.  Mirror  M2  is  used  to  align 
a visible  HeNe  laser  beam  col linear  with  the  DF  laser  beam  and  is  re- 
moved during  measurements.  The  N2O  cell  and  mirror  M7  are  used  in 
detector  calibration  (described  later)  and  are  also  removed  during 
a measurement.  The  apertures  are  used  as  an  aid  in  aligning  the 
visible  HeNe  laser  beam.  Mirror  M5  and  M6  are  used  to  direct  the 
beam  into  the  White  cell.  M6  is  a spherical  mirror  with  0.76  m 
focal  length  which  is  positioned  to  focus  the  DF  laser  beam  at  the 
plane  of  mirror  MIO  in  the  absorption  cell  without  overfilling  mirror 
Mil.  The  beamsplitter  is  an  uncoated  BaFo  wedge  which  reflects  a 
small  percentage  of  the  incident  beam  to  the  reference  detector  which 
is  positioned  near  the  focal  point  of  the  beam.  Mirrors  M8  and  M9 
are  used  to  direct  the  laser  beam  onto  the  signal  detector  after  it 
has  passed  through  the  White  cell,  with  M9  being  used  to  focus  the 
beam  on  the  detector.  All  mirrors  used  in  the  system  are  aluminum 
coated  first  surface  reflectors. 

The  proper  position  of  the  focusing  mirrors  M6  and  M9  and  the 
proper  detector  positions  were  determined  with  the  aid  of  a computer 
program  based  on  the  optical  resonator  chart  developed  by  Collins 
[20].  The  program  is  a generalization  of  one  written  by  Trusty[21  ], 
and  calculates  the  spot  size  and  distance  from  the  focal  point  at 
any  desired  location  in  an  optical  system  consisting  of  a laser  and 
up  to  ten  focusing  elements  (mirrors  o>^  lenses).  A listing  of  the 
program  is  given  in  Appendix  C. 

The  program  assumes  a laser  resonator  consisting  of  a flat  mirror 
on  one  end  and  a spherical  mirror  at  the  output  end.  This  is  equiva- 
lent to  a resonator  with  a flat  grating  on  one  end  and  a spherical 
output  mirror  such  as  the  one  used  with  the  DF  lasers  described 
earlier.  The  program  accounts  for  the  fact  that  the  output  mirror 
is  also  a diverging  lens  v.-'-fan  making  the  calculations. 
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TABLE  5 

FREQUENCIES  OF  OBSERVED  DF  LASER  LINES 


Identification 
Band  Transition 

1-0  P(2) 

1-0  P(3) 

1-0  P(4) 

1-0  P(5) 

1-0  P(6) 

1-0  P(7) 

1-0  P(8) 

1-0  P(9) 

1-0  P(10) 

1-0  P(ll) 

1-0  P(12) 

1-0  P(13) 


Frequency 
(cm-l)  (a) 

2862.653 

2839.791 

2816.380 

2792.434 

2767.968 

2742.998 

2717.539 

2691 .607 

2665.219 

2638.392 

2611.142 

2583.486 


2772.340 

2750.094 

2727.309 

2703.999 

2680.179 

2655.863 


Identification 
Band  Transition 

2-1  P(8) 

2-1  P(9) 

2-1  P(10) 

2-1  P(ll) 

2-1  P(12) 

2- 1  P(13) 

3- 2  P{2) 

3-2  P(3) 

3-2  P{4) 

3-2  P(5) 

3-2  P{6) 

3-2  P(7) 

3-2  P(8) 

3-2  P(9) 

3-2  P(10) 

3-2  P(ll) 

3-2  P(12) 

3-2  P(13) 


Frequency 
(ctir‘)  (a) 

2631 .068 
2605.807 
2580.097 
2553.953 
2527.391 
2500.428 

2683.890 
2662.246 
2640.074 
2617.386 
2594.198 
2570.522 
2546.375 
2521 .769 
2496.721 
2471.245 
2445.356 
2419.070 


(a)  Frequencies  determined  to  ±0.003  cm  ^ by  Heath,  et.al[18] 


Fig.  43.  Block  diagram  of  the  optics  used  with  the 
OF  laser  located  in  the  shielded  room. 


Aligning  the  HeNe  laser  beam  col  linear  with  the  DF  laser  beam  was 
rather  difficult  since  the  DF  laser  beam  was  invisible  and  the  aver- 
age power  was  too  low  to  permit  using  the  fluorescent  screens  which 
are  useful  in  observing  higher  power  infrared  laser  beams. 

The  first  alignitient  method  made  use  of  an  uncoated  calcium  fluo- 
ride flat  in  place  of  mirror  M2  in  Fig.  43.  A small  hand  held  lead 
selenide  detector  was  used  to  find  the  approximate  location  of  the 
DF  laser  beam  on  the  calcium  fluoride  flat.  Mirror  M3  was  then  used 
to  direct  the  visible  laser  beam  to  that  spot.  Then  the  small  de- 
tector was  used  to  find  the  location  of  the  DF  laser  beam  near  mirror 
M5.  The  calcium  fluoride  flat  was  adjusted  so  that  the  HeNe  laser 
beam  struck  the  small  detector.  The  coincidence  of  the  HeNe  and  DF 
laser  beams  was  then  checked  again  near  the  calcium  fluoride  flat 
and  mirror  M3  used  again  to  refine  the  alignment  of  the  visible  beam. 
If  adjustment  was  required,  the  alignment  was  again  checked  near 
mirror  M5  and  adjusted  if  necessary.  The  last  two  steps  in  the  a- 
lignment  were  repeated  until  the  alignment  of  the  two  laser  beams 
was  as  close  as  possible. 

This  method  proved  to  be  somewhat  unsatisfactory  since  the  DF 
laser  spot  was  fairly  large  and  the  center  of  the  spot  was  difficult 
to  determine. 

The  second  method  and  the  one  which  proved  to  be  most  satis- 
factory made  use  of  the  adjustable  apertures  shown  in  Fig.  43.  The 
idea  is  to  exactly  center  the  apertures  on  the  DF  laser  beam  with 
mirror  M2  removed  and  then  insert  mirror  M2  and  adjust  M2  and  M3  so 
that  the  helium  neon  laser  beam  goes  exactly  through  the  apertures. 

The  difficult  part  of  this  method  is  getting  the  apertures 
centered  properly  on  the  DF  laser  beam. 

The  DF  laser  beam  was  directed  to  the  reference  detector  using 
the  small  handheld  detector  to  make  sure  the  beam  fell  entirely  on 
Ml,  M5,  M6,  and  the  barium  fluoride  (BaF^)  beamsplitter.  The  beam- 
splitter was  then  adjusted  for  maximum  signal  on  the  reference  de- 
tector. The  reference  detector  was  located  near  the  focus  of  the  DF 
laser  beam  and  was  large  enough  to  collect  the  entire  beam. 

The  location  and  approximate  size  of  the  DF  laser  beam  at  the 
apertures  was  determined  by  observing  the  reference  detector  signal 
as  the  beam  was  obstructed  by  slowly  moving  a card  into  the  beam 
from  the  top  and  bottom  and  from  either  side.  Using  this  technique, 
the  apertures  could  be  accurately  positioned. 

With  the  HeNe  laser  beam  adjusted  col  linear  with  the  DF  laser 
beam  and  the  apertures  opened  up,  the  rest  of  the  optical  align- 
ment could  be  performed  rather  easily  using  the  HeNe  laser  beam 
as  the  visible  reference. 


Mirrors  M5  and  M6  were  used  to  direct  the  beam  into  the  White 
cell  properly.  The  White  cell  was  adjusted  to  the  desired  path 
length  using  the  HeNe  laser  as  a reference.  The  BaF?  beamsplitter 
was  then  adjusted  to  direct  the  HeNe  laser  beam  to  trie  reference 
detector,  and  mirrors  M8  and  M9  were  used  to  direct  the  laser  beam 
to  the  signal  detector  after  it  ,iad  passed  through  the  White  cell. 
Fine  adjustment  of  the  detector  alignment  was  mcide  with  mirror  M2 
removed  and  the  DF  laser  operating. 

The  second  optical  layout  (Fig.  44)  is  the  one  used  when  the 
DF  laser  was  mounted  on  the  main  optical  tabl a with  the  White  cell 
entrance  optics. 

The  apertures  were  used  to  align  the  visible  HeNe  laser  beam 
col  linear  with  the  DF  laser  as  before. 

Since  the  DF  laser  u^ed  in  this  configuration  had  more  average 
power,  Eppley  thermopiles  w°rr  used  as  the  reference  and  signal  de- 
tectors. Because  the  tppley  thermopiles  are  relatively  slow  they 
could  not  be  easily  used  to  adjust  the  DF  laser  for  optimum  operation 
on  a desired  line.  Therefore  provision  was  made  for  directing  the 
reference  beam  to  a leadselmide  detector  for  laser  adjustment  and 
to  the  Eppley  thermopile  for  absorption  measurements.  This  was  done 
by  mounting  mirror  M7  on  a kinematic  mount  which  could  be  removed 
and  replaced  without  affecting  the  optical  alignment. 


C.  Detectors 


Early  attempts  to  measure  absorption  in  the  White  cell  with 
the  first  pulsed  laser  made  use  of  Eppley  thermopiles  as  the  refer- 
ence and  signal  detectors.  The  average  laser  power  was  too  low 
however  and  air  currents  and  temperature  changes  in  the  room  caused 
detector  responses  which  were  a noticable  fraction  of  the  laser 
signal.  This  resulted  in  a low  signal  to  noise  ratio  and  very  poor 
repeatability  in  the  measurements.  It  was  necessary  therefore  to 
use  fast  detectors  which  could  respond  to  the  peak  laser  power  in 
each  pulse  rather  than  the  average  power. 

The  first  fast  detectors  used  were  mercury  cadmium  telluride 
photoconducti ve  detectors  manufactured  by  Mullard  Inc.  The  detectors 
had  an  operating  temperature  range  of  273-300°K,  a time  constant  of 
0.15  ysec,  and  were  .25  mm  square.  The  detectors  had  adequate  sensi- 
tivity, however  their  small  size  made  it  necessary  that  the  optical 
alignment  be  very  closely  maintained.  This  is  almost  impossible, 
particularly  at  the  output  of  the  absorption  cell  after  the  beam 
has  traveled  three-fourths  of  a kilometer  or  more. 

The  next  detectors  were  lead  selenide  photoconductive  detec- 
tors manufactured  by  Santa  Barbara  Research  Center.  These  detectors 
operated  at  ambient  temperature  ('v  296°K)  with  a time  constant  of 
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1-3  usec  and  were  4 mm  square.  These  detectors  also  had  the  required 
sensitivity  and  their  relatively  large  size  made  alignment  less  criti- 
cal, although  great  care  was  still  needed. 

The  use  of  pulse  detectors  presented  three  additional  problems. 
The  pulse  to  pulse  instability  of  the  laser  made  it  necessary  to  use 
pulse  averaging  electronics,  and  the  propagation  delay  of  about  2.5 
microseconds  in  the  absorption  cell  made  it  necessary  to  look  at  the 
reference  and  signal  pulses  at  different  times  or  to  "save"  the  re- 
ferer.,-e  pulse.  The  solution  to  these  problems  is  discussed  in  the 
next  section. 

The  third  problem  was  the  possibility  of  detector  non-linear- 
ity. Although  the  average  laser  power  was  only  about  1 milliwatt, 
the  peak  pulse  power  was  about  25  or  30  watts.  Accounting  for  the 
fact  that  the  reference  detector  received  a small  percentage  of  this 
power  and  the  signal  detector  only  about  a fourth  (because  of  White 
cell  insertion  loss),  the  detectors  still  received  from  1 to  about 
8 watts  peak  power. 

The  possibility  of  detector  non-linearity  was  investigated  by 
intercomparing  the  reference  detector  with  an  Eppley  thermopile. 

While  the  thermopile  detectors  could  not  be  used  to  make  the  ab- 
sorption measurements  because  of  the  low  reflectivity  of  the  beam- 
splitter and  White  cell  insertion  loss,  it  was  possible  to  use  a 
thermopile  to  make  the  detector  comparison  by  placing  it  where  it 
received,  most  of  the  laser  power.  The  Eppley  thermopile  was  known 
to  be  linear  in  this  low  average  power  region,  so  if  the  lead  sele- 
nide  detector  was  also  linear, a plot  of  voltage  from  the  thermopile 
detector  versus  voltage  from  the  lead  selenide  detector  would  be  a 
straight  line.  In  fact  the  relationship  observed  was  not  a straight 
line,  so  the  lead  selenide  detector  was  indeed  non-linear  as  sus- 
pected. 


The  detectors  were  calibrated  against  the  Eppley  thermopile 
using  the  following  procedure. 

The  experimental  layout  is  that  shown  in  Fig.  43.  The  idea 
is  to  compare  the  reference  detector  PbSe  ^^1  to  the  Eppley  thermo- 
pile over  a wide  range  of  signal  levels  with  mirror  M7  in  place. 

Then  the  mirror  is  removed  and  the  reference  detector  and  signal 
detectorPbSe  #2  are  compared  with  the  White  cell  evacuated. 

At  first  ’t  was  thought  that  the  signal  level  to  the  detectors 
could  be  changed  by  changing  the  gas  mix  and  excitation  voltage  in 
the  laser.  This  procedure  gave  unsatisfactory  and  unrepeatable  re- 
sults. This  could  have  been  caused  by  pulse  shape  changes  as  gas 
mix  and  discharge  voltage  are  changed  since  the  electronics  associ- 
ated with  the  lead  selenide  detectors  responded  to  only  the  peak  of 
the  pulse  while  the  thermopile  responded  to  the  energy  in  the  pulse. 
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It  was  necessary  therefore  to  use  an  external  attenuator  and 
leave  the  laser  undisturbed  during  the  calibration  runs.  The  sig- 
nal to  the  detectors  was  attenuated  step  by  step  by  introducing  N?0 
into  the  21  cm  ^.ell  shown  in  Fig.  43.  The  first  part  of  the  experi- 
ment comparing  the  thermopile  and  the  reference  detector  PbSe  #] 
was  repeated  three  times  on  three  separate  days,  as  was  the  second 
part  of  the  experiment  comparing  the  reference  detector  and  the  sig- 
nal detector  PbSe  #2.  At  each  signal  level  the  signals  from  the 
two  detectors  were  read  and  recorded  on  the  typewriter  by  the  com- 
puter using  the  data  taking  program  which  will  be  described  later. 

A least  squares  fit  of  the  data  from  the  first  part  of  the 
experiment  was  made  to  a third-order  polynomial.  This  results  in 
an  expression  of  the  form 


(107)  E]  = ^ + A^X^^  + A3X^^ 


relating  the  reference  detector  voltage  X]  to  the  Eppley  thermopile 
voltage  E]  (Fig.  45).  The  voltages  are  those  read  by  the  computer 
after  arbitrary  fixed  amplifications. 

Similarly  the  data  from  the  second  part  of  the  experiment  was 
used  to  obtain  an  expression  of  the  form 


(108)  = Cq  + C^X^  + C^X^^  + 63X2^ 


where  X2  is  the  signal  detector  voltage  as  'eaJ  by  the  computer  and 
X]  is  the  reference  detector  voltage  (Fig.  16  . 

From  Eqs.  (107)  and  (108)  the  followng  expression  relating 
the  signal  detector  voltage  X2  to  an  equiviient  Eppley  voltage  Ep 
v/as  determined  (Fig.  47):  ^ 


(109)  E^  = Bg  + B^X^  + B^X^^  + 63X2^ 


The  linearity  correction  expressions  determined  by  the  above 
procedure  were  tested  by  incorporating  them  into  the  data-taking 
program  and  repeating  the  second  part  of  the  calibration  experiment. 
At  each  signal  level  the  computer  read  the  voltages  Xi  and  X„  from 
the  reference  and  signal  detectors  and  usee  Eqs.  (107)  and  (t09)  to 
determine  E]  and  E2  respectively.  If  the  linearity  corrections  are 
good,  a plot  of  E2  versus  E]  should  be  a straight  line.  Figure  48 
shows  that  this  was  indeed  the  case. 
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PbSe  #I(X,) 


Fig,  45.  Eppley  thermopile  calibration  signal  versus  lead 
selenide  detector  number  one,  Units  are  volts 
at  A/D  converter  after  arbitrary  amplifications, 


Fig,  48.  Results  of  experiment  to  check  calib*'ation 
linearity. 


The  first  attempt  to  determine  the  linearity  correction  used 
quadratic  rather  than  cubic  expressions.  It  was  found,  however, 
that  the  quadratic  expressions  were  not  sufficient  to  correct  the 
nonlinearity. 

Subsequent  calibrations  used  a somewhat  different  procedure. 
The  experimental  technique  used  was  the  same.  However,  before  the 
second  part  of  the  experiment  was  performed,  the  data  from  the  first 
part  of  the  experiment  was  used  to  obtain  the  expression  in  Eq. 

(107).  This  expression  was  incorporated  into  the  data  taking  pro- 
gram. The  data  from  the  second  part  of  the  experiment  was  then  in 
the  form  corresponding  to  Fig.  47  or  Eq.  (109)  rather  than  Fig.  46 
or  Eq.  008).  The  two  procedures  yield  the  same  results,  however 
the  second  requires  less  hand  manipulation  of  the  data. 

The  detectors  described  above  were  used  with  the  original 
laser  mounted  in  the  shielded  room.  The  small  portable  laser  which 
was  built  later  could  be  mounted  on  the  main  optical  table  and  had 
four  or  five  times  as  much  average  power.  This  made  possible  the 
use  of  thermopile  detectors.  This  was  desirable  since  the  thermo- 
pile detectors  respond  to  the  average  laser  power  and  they  are  known 
to  be  linear  in  the  low  average  power  region. 

Some  care  was  required  in  the  use  of  the  thermopile  a^tectors 
since  they  readily  respond  to  room  temperature  fluctuations  and  air 
currents.  To  help  eliminate  these  problems  the  thermopiles  were 
each  encased  in  a 4 inch  block  of  styrofoam  with  a small  hole  to 
admit  the  laser  beam.  Boxes  were  also  constructed  of  1 inch  thick 
styrofoam  to  completely  enclose  the  detectors,  with  a small  hole 
cut  in  the  boxes  to  admit  the  laser  beam. 

The  styrofoam  boxes  and  insulating  blocks  protected  the  de- 
tectors from  short  term  temperature  fluctuations  and  air  currents. 
However  ’:he  detectors  did  respond  to  long  term  temperature  changes. 
Therefore  if  the  room  thermostat  setting  was  changed  or  a heat  load 
was  suddenly  introduced  in  the  room  it  was  necessary  to  wait  until 
the  temperature  had  stabilized  before  attempting  to  make  measure- 
ments. This  usually  took  no  more  than  thirty  minutes  to  an  hour. 


D.  Electronics 


Another  important  part  Of  the  system  was  the  electronics  used 
to  amplify  and  record  the  detector  signals.  There  were  two  separate 
systems  used  with  the  two  different  types  of  detectors.  The  lead 
selenide  detectors  which  were  fast  enough  to  respond  to  each  laser 
pulse  required  the  use  of  commercial  box- car  integrators  and  some 
other  specially  built  electronics.  The  Eppley  thermopiles  were 
relatively  slow  and  therefore  responded  only  to  the  average  laser 
power.  DC  amplifiers  were  used  with  these  detectors.  The  output 
of  the  electronics  in  either  case  was  read  by  an  XDS  920  computer 
through  a 24  channel,  14  bit  A/D  converter. 
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The  boxcar  integrators  and  associated  electronics  will  be 
described  in  part  1,  the  DC  amplifiers  will  be  described  in  part  2, 
and  the  computer  and  A/D  converter  will  be  described  in  part  3 below. 

1 . Boxcar  integrators  and  associated 
electronics 

Boxcar  integrators  were  needed  with  the  fast  lead  selenide 
detectors  for  two  reasons:  Ij  to  average  out  variations  caused 

by  pulse  to  pulse  instability  and  2)  to  provide  a DC  voltage  pro- 
portional to  the  pulse  height  which  could  be  read  by  the  computer. 

A block  diagram  of  the  pulse  detector  electronics  is  shown  in 
Fig.  49.  The  pulse  amplifiers,  dual  gate  generator,  and  differential 
gated  integrators  were  manufactured  by  Molectron  Corporation.  The 
2-channel  peak  holding  circuit,  the  optical  trigger  coupler,  and  the 
detector  bias  circuits,  were  specially  built  for  this  particular 
experiment. 

The  key  components  of  the  system  are  the  differential  inte- 
grators (Molectron  Mdl . 112).  Typical  operation  of  one  of  the  dif- 
ferential gated  integrators  is  depicted  in  Fig.  50  and  Fig.  51. 

The  diode  bridges  are  back-biased  when  the  gates  are  off. 

When  a gate  is  on  (gate  input  greater  than  2 volts)  the  diode  bridge 
is  forward  biased  and  the  integrator  side  of  the  bridge  follows  the 
signal  side.  Thus  if  the  voltage  on  the  integrator  capacitor  is 
less  than  the  signal  voltage,  current  flows  from  the  bridge  to  charge 
the  capacitor,  and  if  the  signal  voltage  is  less  than  the  capacitor 
voltage,  the  capacitor  discharges.  In  the  back-biased  mode  the 
diodes  have  extremely  high  impedance,  and  t'"  buffer  amplifier  follow- 
ing the  integrator  has  extremely  high  impeoanc^.  These  factors  re- 
duce leakage  currents  from  the  integrating  capacitor  to  less  than 
0.5  pi  coamperes. 

The  integrating  time  constant  is  selected  by  changing  R and 
C in  the  integrator.  If  the  input  is  a continuous  string  of  pulses 
as  shown  in  Fig.  51,  and  the  two  gates  are  adjusted  as  shown,  the 
output  is  just  the  pulse  height  and  varies  with  a time  constant 
equal  to  RC(T/Xg). 

The  gate  pulses  and  the  input  signal  are  superimposed  at  the 
monitor  jack  for  convenience  in  adjusting  the  gate  timing.  The 
signal  at  the  monitor  jack  is  one  tenth  the  input  signal  and  the 
gate  pulses  are  fixed  at  25  millivolts. 

The  time  constants  of  the  integrators  in  both  channels  are 
matched  to  within  0.5% 
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The  output  of  the  gated  integrators  was  fed  directly  into 
the  A/D  converter  on  the  XDS  920  computer. 

The  dual  gate  generator  (Molectron  Model  122)  accfepted  the 
trigger  pulse  and  generated  the  gate  signals  required  by  the  two 
differential  integrators.  The  delay  time  from  the  trigger  pulse  to 
gate  1,  the  gate  separation,  and  the  gate  width  are  all  adjustable. 

The  response  of  the  lead  selenide  photoconductive  detectors 
to  an  increase  in  the  intensity  of  infrared  radiation  striking  them 
is  an  increase  in  electrical  conductivity.  The  primary  purpose  of 
the  bias  circuit  is  to  convert  this  increase  in  conductivity  to  a 
voltage  change.  A circuit  diagram  of  the  bias  circuits  used  in 
this  study  is  shown  in  Fig.  52.  R1  is  adjusted  to  maintain  the  bias 
voltage  at  160v  as  the  battery  ages.  Cl  is  a filter  capacitor,  R2  is 
the  load  resistor,  C2  insures  that  the  steady  state  voltage  at  the 
output  connector  is  zero,  R3  and  C3  integrate  the  pulse,  and  S2  is 
used  to  short  the  output  and  prevent  high  voltage  transients  when 
SI  is  opened  or  closed. 

The  signals  from  the  detector  bias  circuits  were  processed  by 
the  pulse  amplifiers  (Molectron  Model  131)  which  shaped  the  pulses 
and  amplified  them  up  to  five  volts. 
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Fig,  52,  Lead  selenide  detector  bias  circuit  diagram, 
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One  problem  which  had  to  be  solved  was  the  delay  between  the 
reference  pulse  and  the  signal  pulse  caused  by  the  transit  time  of 
the  White  cell.  The  solution  to  this  problem  was  a peak-holding 
circuit  which  held  the  peak  of  each  pulse  until  after  gate  1 was 
off.  This  circuit  is  shown  in  Fig.  53.  When  the  input  voltage 
exceeds  the  output  voltage  D1  conducts  and  charges  Cl  to  the  new 
input  voltage.  As  the  input  voltage  decreases,  D1  becomes  reversed 
biased  and  the  peak  voltage  is  held  on  Cl  (discharging  slowly 
through  R5).  R1  in  series  with  D1  limits  the  charging  current.  To 
discharge  Cl,  a positive  pulse  of  approximately  2 volts  (Gate  1)  is 
applied  to  the  reset  terminal.  Q1  inverts  the  reset  pulse  and  C2 
and  R4  differentiate  the  inverted  reset  pulse.  D2  suppresses  the 
negative  spike  from  the  differentiator  (caused  by  the  leading  edge 
of  the  reset  pulse).  The  positive  spike  (caused  by  the  trailing 
edge  of  the  reset  pulse)  causes  the  two  FET's  to  conduct  and  dis- 
charge the  storage  capacitors  Cl. 

Another  problem  was  caused  by  the  very  large  high  frequency 
noise  associated  with  the  trigger  pulse  from  the  laser.  An  attempt 
was  made  to  suppress  the  noise  with  an  RC  filter.  This  attempt  was 
unsuccessful.  The  trigger  pulse  was  then  coupled  optically  through 
the  screen  room  wall  using  an  LED  and  a photo  transistor  as  shown 
in  Fig.  54.  This  optically  coupled  trigger  worked  quite  well  and 
eliminated  the  noise  completely. 
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Fig,  53.  Peak  holding  circuit  diagram. 
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Fig.  54,  Optical  trigger  circuit  diagram, 


2.  DC  amplifiers 

When  the  thermopiles  were  used,  the  required  electronics  was 
much  simpler  since  the  only  requirement  was  that  the  thermopile 
signal  be  amplified  to  about  a plus  or  minus  ten  volt  range  without 
introducing  noise  so  that  full  advantage  could  be  taken  of  the  A/D 
converter  resolution. 

The  signals  were  first  fed  into  HP425A  microvoltmeters  whose 
output  was  a maximum  of  plus  or  minus  one  volt.  The  outputs  of  the 
microvoltmeters  were  then  amplified  to  plus  or  minus  ten  vclts  by 
Alinco  Model  518A  differential  amplifiers.  The  output  of  the  Alinco 
amplifiers  was  read  directly  by  the  A/D  converter  on  the  XDS  920 
computer. 


3.  A/D  converter  and  computer 

The  XDS  920  computer  used  in  this  experiment  has  4096  24  - 
bit  words  of  core  memory  with  a cycle  time  of  8 microseconds  and 
a fairly  powerful  instruction  set, and  input/output  capability  which 
makes  it  easily  adaptable  to  laboratory  experimental  applications. 
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The  cycle  time  of  8 microseconds  is  slow  compared  with  modern  com- 
puters, however  the  speed  is  quite  adequate  for  the  present  appli- 
cation . 

The  analog  to  digital  converter  was  manufactured  by  Epsco, 
Incorporated.  It  has  a 24  channel  multiplexer  with  an  input  imped- 
ance of  100,000  ohms  and  an  input  range  from  -10  to  10  volts.  The 
digital  output  word  is  14  bits  long  which  means  the  resolution  is 
about  1.24  millivolts  per  bit. 

Programs  for  the  computer  are  written  primarily  in  Fortran 
except  for  some  subroutines  which  control  special  devices  such  as 
the  analog  to  digital  converter  which  are  written  in  assembly 
language.  The  program  used  to  take  the  data  in  this  study  is  de- 
scribed in  the  next  chapter  along  with  the  experimental  procedure. 


E.  The  Absorption  Cell 

The  absorption  cell  used  in  this  study  was  originally  designed 
and  built  by  Long[22].  The  cell  is  0.61  meter  in  diameter  and  16.15 
meters  long  with  a sample  volume  of  4.72  cubic  meters  or  4720  liters. 
The  cell  walls  were  honed  smooth  during  construction  in  order  to 
reduce  water  vapor  adsorption  effects. 

Path  lengths  of  well  over  one  kilometer  are  easily  obtainable 
in  the  cell  using  a three  mirror  optical  system  of  the  type  first 
described  by  White[23].  The  White  optical  system  used  in  this  cell 
is  shown  in  Fig.  55.  The  mirrors  ail  have  a radius  of  curvature  of 
15.24  meters,  and  mirrors  A and  B are  separated  from  mirror  C by 
the  common  radius  of  curvature  15.24  meters.  Mirrors  A and  B were 
obtained  by  cutting  a 51  centimeter  diameter  mirror  into  two  halves. 
Mirror  C is  30.48  centimeters  in  diameter  with  notches  cut  out  as 
shown  for  the  entrance  and  exit  beams. 

The  number  of  traversals  and  nence  the  path  length  is  adjusted 
by  tilting  mirrors  A and  B as  shown  in  Fig.  55.  The  external  optics 
are  adjusted  so  that  the  incoming  laser  beam  is  focused  just  as  it 
passes  the  front  surface  of  mirror  C and  with  proper  optics  exter- 
nal to  the  cell  it  diverges  to  almost  fill  mirror  A.  The  White 
optical  system  then  has  the  property  that  the  beam  is  focused  at 
the  plane  of  mirror  C after  every  second  traversal  of  the  cell  if 
the  distance  from  mirror  C to  mirrors  A and  B is  exactly  the  radius 
of  curvature  of  mirrors  A and  B.  The  number  of  traversals  is  deter- 
mined by  counting  the  spots  on  mirror  C.  The  number  of  traversals 
is  just  two  more  than  twice  the  number  of  spots. 

The  entrance  and  exit  windows  of  the  cell  are  barium  fluoride 
half-degree  wedges.  Wedges  are  used  rather  than  flats  to  eliminate 
interference  effects. 


95 


Fig,  55,  Ray  diagram  of  White  cell  (8  traversals). 


For  most  of  the  measurements  reported  in  the  next  chapter,  the 
cell  was  set  for  48  traversals  or  731.7  meters.  For  some  of  the 
measurements  however  a longer  path  was  required.  The  cell  was  there- 
fore set  at  88  traversals  or  1.341  kilometers.  It  was  discovered 
that  at  the  longer  path  length  the  output  spot  was  so  large  that  the 
spots  on  mirror  C overlapped.  Using  the  spot  size  program  described 
earlier,  it  was  discovered  that  a small  error  in  mirror  separation 
would  cause  the  spot  on  mirror  C to  become  larger  as  the  number  of 
traversals  became  greater. 

The  mirror  separation  was  then  adjusted  using  a white  light 
source.  An  image  was  formed  at  point  0 in  the  plane  of  mirror  C 
which  then  diverged  to  fill  mirror  A.  The  image  was  observed  in 
the  plane  of  mirror  C at  the  output  with  the  cell  set  for  64  passes. 
The  mirror  separation  was  then  adjusted  to  make  that  image  as  sharp 
as  possible.  It  was  found  that  the  mirrors  had  been  about  1 cm  too 
far  apart.  Since  the  mirror  separation  is  15.24  meters,  this  is  not 
a large  error.  It  does  however  cause  problems  if  long  path  lengths 
are  desired. 

From  the  above  discussion  it  is  obvious  that  for  a White 
cell  which  is  designed  to  be  either  heated  or  cooled  some  provision 
must  be  made  in  the  design  for  keeping  the  mirrors  properly  separated 
as  the  cell  expands  and  contracts. 

The  cell  is  evacuated  using  a 100  cubic  foot  per  minute  mechani- 
cal pump  and  a six  inch  diffusion  pump  connected  near  the  center  of 
the  cell  through  pneumatically  operated  valves. 

Gases  are  admitted  to  the  cell  through  three  ports  located 
near  the  center  of  the  cell  and  at  each  end.  The  ports  are  con- 
nected to  a comnon  manifold  where  various  gases  can  be  admitteed. 

There  are  two  small  fans,  one  at  each  end,  inside  the  cell 
to  aid  in  mixing  the  gas  samples. 

There  are  several  different  types  of  vacuum  gauges  attached  to 
the  cell.  The  vacuum  pumps'  operation  is  monitored  by  a set  of 
thermocouple  gauges  and  an  Alphatron  530  gauge.  Gas  sample  pressure 
is  monitored  by  three  separate  mercury  gauges.  Sample  pressures  up 
to  two  torr  were  measured  with  a McLeod  gauge.  Sample  pressures 
from  1 torr  to  50  torr  are  measured  with  a Roger  Gilmont  Instruments 
Model  906  mercury  micrometric  manometer.  This  instrument  can  be 
read  to  .025  torr.  Prsssures  up  to  1000  torr  are  measured  with  a 
U-tube  manometer  which  can  be  read  to  1 torr. 

The  dew-point  of  the  water  vapor-air  or  water  vapor-nitrogen 
mixtures  was  measured  with  a dew-point  hygrometer  which  optically 
senses  the  formation  of  dew  or  frost  on  a thermo-electrical ly  cooled 
metal  mirror.  Two  different  hygrometers  were  used.  One  was  a Cam- 
bridge Systems  Model  880  which  used  a precision  aged  thermistor  to 


97 


measure  the  dew-poin^  ‘ mperature.  This  instrument  was  calibrated 
against  a mercury  mi  i ?tric  manometer  in  this  laboratory.  The 
other  hygrometer  usef.  ,,  a Cambridge  Systems  Model  992  which  used 
a precision  platinum  resistance  thermometer  to  sense  the  dew-point 
temperature.  The  calibration  on  this  instrument  is  traceable  to 
the  National  Bureau  of  Standards. 

The  cell  temperature  was  monitored  using  Stow  Laboratories 
platinum  resistance  thermometers  installed  near  the  center  of  the 
cell  and  at  either  end. 

Also  connected  to  the  cell  v>as  a mass  spectrometer  residual 
gas  analyzer  which  could  be  used  to  monitor  the  composition  of  a 
gas  sample  over  a period  of  time  to  determine  whether  selective 
adsorption  might  be  occurring. 


CHAPTER  IV 


DISCUSSION  OF  MEASUREMENTS 


Absorption  measurements  were  made  for  five  absorbers  and  eight 
different  DF  laser  lines.  The  absorbers  studied  were  N?0,  CH^,  C02» 
HDO,  and  HoO.  The  laser  lines  studied  were  the  2-1  P(6),  P(7],  P(8), 
P(10),  and  P(ll)  lines  and  the  3-2  P(6),  P(7),  and  P{8)  lines. 

Absorption  by  all  constituents  was  not  measured  for  each  line. 
For  the  2-1  P(6)  - P(8)  lines  and  the  3-2  P(6)  - P(8)  lines  the  goal 
was  to  measure  the  molecular  absorption  accurately.  Therefore  only 
those  constituents  which  caused  significant  absorption  for  these 
lines  were  measured.  This  determination  was  made  by  studying  the 
calculations  described  in  Chapter  II  and  preliminary  measurements 
made  by  Spencer,  et  al.[24].  In  addition  NpO  absorption  was  also 
measured  on  the  2-1  P(10)  and  P(ll)  lines  since  they  were  the  only 
other  lines  for  which  N2O  absorption  was  significant. 

Also  a laser  spectroscopy  technique  was  used  to  determine  the 
frequencies  of  the  3-2  P(7)  and  2-1  P(10)  lines  by  measuring  their 
separation  from  N2O  absorption  lines. 


A.  N2O  Absorption  Measurements, 

760  Torr,  23~^ 

Measurements  of  the  absorption  by  nitrous  oxide  samples  broad- 
ened to  760  torr  with  dry  nitrogen  were  made  for  the  3-2  P(6)  - P(8) 
lines  and  the  2-1  P(10)  and  P(ll)  lines.  For  these  measurements  the 
White  cell  was  set  for  a path  length  of  0.7317  kilometers. 

Assuming  the  calculations  in  Chapter  II  are  in  error  by  no  more 
than  ±100%,  it  is  not  possible  to  measure  the  N2O-N2  absorption  coeffi- 
cients with  acceptable  accuracy  in  a path  of  0.731/  kilometers  unless 
the  N2O  concentration  is  increased  above  the  normally  assumed  sea- 
level  concentration  of  0.28  ppm.  Therefore  the  N2O  concentrations  used 
were  from  10  to  1000  times  the  normal  concentration  (2.8  to  280  ppm). 
These  concentrations  are  still  quite  low  and  the  results  can  be 
linearly  extrapolated  to  normal  N2O  abundance  without  distortion 
which  might  arise  from  the  difference  between  the  sel f -broadened  and 
foreign-broadened  half-width.  This  can  be  seen  from  the  following. 
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Assuming  the  Lorentz  line  shape,  the  absorption  coefficient  in 
km"^  may  be  written  as  follows 

-IS  u 

(110)  k(km^)=-- -7- ^ 

(v-v  ) + a, 


where  u is  the  absorber  concentration  in  molecules-cm 


Now  uis  directly  proportional  to  the  absorber  partial  pressure?- 
and  a(_  is  ^iveh  for  1%K  from  Eqs.  (67)  and  (68)  by  aLo 
Therefore  the  absorption  coefficient  can  be  written 


(111) 


k(km  ) = — 


“Lo'i 

+ “LO^  (Pt+ 


where  C is  a proportionality  constant  relating  u to  P . If  the 
absorber  pressure  P,  is  very  small  compared  to  the  total  pressure, 
p_  + (B-l)P  is  very  nearly  Pj  and  the  absorption  coefficient  k is 
directly  proportional  to  the  absorber  pressure  P^.  Thus  if  the 
absorption  coefficient  is  known  at  some  absorber  pressure  P^,  the 
absorption  coefficient  at  some  lower  pressure  P^  is  just  PL/Pa  times 
the  absorption  coefficient  at  P^. 

For  each  experimental  run,  the  absorption  was  measured  on  a 
single  DF  laser  line  for  a number  of  NgO  concentrations.  First  the 
cell  was  evacuated  and  the  signal  detector  and  reference  detector 
voltages  were  read  by  the  computer  and  their  ratio  calculated  and 
recorded  on  the  computer  typewriter.  Any  detector  linearity  cor- 
rection determined  from  calibration  experiments  was  made  automatic- 
ally at  this  time.  Next  a predetermined  amount  of  a mixture  of  1% 

NpO  in  Nitrogen  was  admitted  to  the  cell  and  the  cell  filled  to 
7b0  torr  with  dry  nitrogen.  The  1%  mixture  was  :,ised  to  admit  the 
sample  rather  than  pure  N2O  so  the  absorber  amount  could  be  accurately 
measured.  Even  at  the  highest  N2O  concentration  used  (250  ppm)  the 
N2O  partial  pressure  was  still  only  0.2  torr.  With  the  mercury 
pressure  gauges  which  are  used  on  the  absorption  cell,  pressures 
between  1 and  50  torr  can  be  measured  with  greater  accuracy  than 
pressures  less  than  1 torr.  Therefore  a mixtrue  of  1%  N2O  in  nitro- 
gen was  used  as  the  absorbing  gas. 

The  cell  now  contained  the  absorbing  mixture  at  the  highest 
N2O  concentration  to  be  measured  in  the  current  run.  This  sample 
was  permitted  to  mix  for  about  one  hour.  Then  the  ratio  of  signal 
voltage  to  reference  voltage  was  read  and  recorded  by  the  computer 
as  before.  The  proper  mixing  time  was  determined  by  monitoring  the 
transmittance  of  the  cell  continuously  until  it  was  observed  to  be 
stable.  Next  the  mixture  was  partially  pumped  out,  the  pressure 
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measured,  and  the  cell  refilled  with  nitrogen  to  760  torr.  This 
procedure  reduced  the  N2O  concentration  by  a known  amount.  This 
new  sample  was  allowed  to  mix  for  about  one  hour  and  the  signal  to 
reference  ratio  was  read  and  recorded  as  before.  This  procedure 
was  repeated  to  obtain  data  at  several  N2O  concentrations.  At  the 
end  of  the  day,  the  cell  was  evacuated,  and  the  signal  to  reference 
ratio  recorded. 

The  above  procedure  was  repeated  several  times  on  different 
days  for  each  laser  line.  All  measurements  were  made  at  room  tem- 
perature 23°C) . 

Although  the  laser  was  tunable  from  line  to  line,  it  was  found 
the  desired  repeatability  in  transmittance  could  not  be  obtained  un- 
less the  laser  was  left  tuned  to  one  line  for  the  duration  of  a 
measurement  run.  This  is  presumably  due  to  optical  alignment  diffi- 
culties associated  with  moving  the  grating. 

For  these  measurements  the  optical  arrangement  shown  in  Fig.  43^ 
with  the  DF  laser  in  the  shielded  room  and  the  lead  selenide  detectors 
was  used.  The  procedure  used  for  making  each  measurement  took  ad- 
vantage of  the  averaging  capability  offered  by  the  computer.  Before 
each  measurement,  the  laser  beam  was  blocked.  On  command  from  the 
typewriter  the  computer  then  read  the  signal  channel  and  reference 
channel  300  times  in  about  15  seconds,  computed  the  averages  and 
recorded  them  on  the  typewriter.  These  zero  readings  represented 
the  offset  introduced  by  the  electronics  and  were  subtracted  from 
all  subsequent  readings.  Next  the  beam  was  unblocked  and  the  signal 
and  reference  readings  were  allowed  to  stabilize.  On  command  from 
the  typewriter  the  computer  then  made  a series  of  five  measurements. 
For  each  measurement,  the  signal  and  reference  voltage  were  read 
100  times  in  a period  of  about  5 seconds  and  averaged,  the  zero  off- 
sets measured  before  were  subtracted,  and  the  readings  were  corrected 
for  linearity  using  the  cubic  equation  for  each  detector  derived  from 
the  calibration  experiments.  The  computer  then  printed  the  corrected 
signal  voltage  and  reference  voltage,  the  standard  deviations  of  each 
measurement,  and  the  ratio  of  signal  to  reference  voltage.  After 
this  had  been  repeated  five  times,  the  average  signal  voltage, 
reference  voltage  and  ratio  were  recorded  on  the  typewriter.  A 
sample  data  run  is  shown  in  Fig.  56.  Five  sets  of  100  samples  were 
used  rather  than  one  set  of  500  samples  so  that  short  term  trends 
in  the  readings  such  as  zero  drifts,  alignment  instabilities  or  laser 
power  changes  could  be  observed.  All  of  the  parameters  of  the  mea- 
surement procedure  could  be  modified  by  the  user  such  as  number  of 
measurements,  number  of  samples  and  delay  between  sample  readings. 

In  addition  a threshold  value  of  standard  deviation  could  be  set 
which  if  exceeded  caused  the  measurement  to  be  ignored  in  the  final 
average. 

Figures  5/  - 61  show  the  measured  absorption  coefficients  versus 
N2O  concentration  for  the  five  lines  studied.  In  each  case  a least 
squares  fit  of  the  data  to  a straight  line  through  the  origin  was 
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9. 

nATA  TAKING  PR0GRAI1.  RCVISEO  OCT.  25,1975 


TODAYS  DATE  IS: 
ll/lC/73 


VARIABLE  VALUES  ARE: 

ADI  3 

AD2  18 

DLY  50.000 

SAM  100 

ASTDV  1.0000 

NO  A 5 


-OPTION 

17. 

ENTER  DETECTOR  CALIBRATION  COEFFICIENTS  FOR  THE  CUBIC  Y»A0AAl*XaA2*X**2RA3*X**'' 
ENTER  AO,..., A3  FDR  SIGNAL  DETECTOR. 

.2 51 1,1. 71 4 7,  1.8 98  7, -.  52 312, 

ENTER  AO,..., A3  FOR  REFERENCE  DETECTOR. 

.11762,1.3011, .729 59, -.015012, 

-DPTIDN 

5. 

ZERO  ADJUST.  BLOCK  BEAM,  CARRIAGE  RETURN  TO  CONTINUE. 


ZEROl  IS 
ZERD2  IS 


-.0418 

.0005 


STANDARD  DEVIATION  IS 
STANDARD  DEVIATION  IS 


-OPTION 

5. 

ZERO  ADJUST.  BLOCK  BEAM,  CARRIAGE  RETURN  TD  CONTINUE. 


ZERDl  IS  -.0425 
ZER02  IS  .0455 


STANDARD  DEVIATION  IS 
STANDARD  DEVIATION  IS 


.0268 
. 0425 


.0266 

.0140 


-DPTIDN 

5. 

ZERO  ADJUST.  BLOCK  BEAM,  CARRIAGE  KciUKN  TO  CONTINUE. 


ZERDl  IS 
ZERD2  IS 


-.0424 

.0188 


STANDARD  DEVIATION  IS 
STANDARD  DEVIATION  IS 


. 0265 
.0286 


-OPTION 


CHANNEL 

3 

CHANHFL 

18 

SIGNAL 

REFERENCE 

VOLTAGE 

STD  DEV 

VOLTAGE 

STD  DFV 

SIG/REF 

1 

5.4207 

.0271 

6.9426 

.0269 

.78078 

2 

5.4173 

.0266 

6.9555 

.0267 

. '7910 

3 

5.4770 

.0267 

6.9186 

.0270 

. 79163 

4 

5.4484 

.0265 

6.9521 

.0266 

. 78370 

5 

5.4292 

.0265 

6.9077 

.0264 

, 78595 

AVERAGE: 

5.4385 

6.9549 

. 78423 

TAKEN  5 TIMES 


Fig.  56,  Example  run  of  data-taking  program, 
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Fig.  57.  Measured  N2O  absorption  for  the  3-2  P(6) 
line  at  2594,198  cnri  (23°C), 
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Fig.  59.  Measured  N2O  absorption  for  the  3-2  P(8) 
line  at  2546,522  cm"l  (23°C). 


made.  The  fact  that  the  data  points  fall  almost  on  a straiqht  line 
is  the  best  indication  of  the  validity  of  the  detector  linearity 
correction  procedure  described  in  Chapter  III.  Table  6 gives  the 
measured  absorption  coefficient  per  part  per  million  along  with  the 
calculated  value  from  Table  IV.  The  absorption  coefficients  extra- 
polated to  the  normally  assumed  atmospheric  abundance  of  0.28  ppm 
are  also  given  alongwith  thei^esults  obtained  by  Spencer(24  ] and 
Deaton[  25] . 

The  measured  absorption  coefficients  and  calculated  absorption 
coefficients  agree  within  1 or  2 percent  except  for  the  2-1  P(ll) 
line  where  the  difference  is  about  10%. 


The  agreement  with  Spencer's  measurements  is  not  as  good. 
Spencer's  measurements  were  made  with  pure  N2O  in  a 1 or  10  cm  cell . 
The  values  given  in  Table  6 are  extrapolated  from  Spencer's  meas- 
urements at  760  torr  assuming  a self-broadening  coefficient  for 
N2O  of  1 . If  the  self-broadening  coefficient  in  not  1,  this 
extrapolation  could  introduce  a substantial  error, 


Deaton's  measurements  were  made  using  a differential  spectro- 
phone  with  flat  windows  and  N2O  concentrations  from  1 to  1000  times 
the  normally  assumed  atmospheric  concentration.  The  spectrophone 
was  calibrated  using  the  N2O  data  from  this  study  for  the  3-2  P{6), 

P(7)  and  P(8)  lines.  The  calibration  factors  determined  by  Deaton 
were  different  for  the  three  lines.  He  attributed  this  difference 
to  experimental  errors  and  took  the  average  of  the  three  as  the 
correct  value.  The  fairly  good  agreement  for  the  3-2  P{6),  P{7), 
and  P(8)  lines  then  is  to  be  expected  since  Deaton  used  the  data 
from  this  study  for  those  three  lines  as  his  calibration.  The  data 
for  the  2-1  P(10)  line  is  not  in  good  agreement.  There  are  a couple 
of  things  which  could  account  for  this.  The  data  in  this  study  is 
valid  for  760  torr  total  pressure.  Deaton's  measurements  and  cali- 
brations were  made  at  630  torr.  There  is  no  certainty  that  the  760  torr 
data  could  be  used  to  calibrate  at  630  torr  and  calibration  errors 
could  be  introduced  which  were  different  for  different  lines.  Also 
Peterson[26]  has  found  that  the  spectrophone  calibration  can  be 
frequency  dependent  since  the  flat  spectrophone  windows  are  also 
low-finesse  Fabry-Perot  etalons. 

For  two  of  the  laser  lines,  3-2  P{7)  and  2-1  P(10),  the  assumed 
laser  frequency  differs  from  the  assumed  frequency  of  an  N2O  absorp- 
tion line,  In  this  situation,  it  is  possible  to  determine  para- 
meters of  the  absorption  line  from  a plot  of  absorption  coefficient 
versus  total  pressure  with  the  absorber  pressure  constant. 

This  can  be  shown  by  the  following[27].  For  an  isolated  Lorentz 
line  the  following  expression  may  be  written 


(112) 


g 

(Av)^  + 


2 
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Measured  N2O  Absorption  Coefficients  for  Five  DF  Laser  Lines 


E 

s-  <= 

0-0  — ' 
00 

rtJ  ' 

-D  • fl)  E 


oo 

00 

CO 

Oi 

LD 

Oi 

0 

0 

CO 

r— 

0 

0 

0 

0 

• 

• 

• 

• 

0 

0 

0 

0 

CO 

0 

LO 

CO 

CO 

CO 

CO 

LD 

CO 

0 

CO 

CO 

r— 

0 

0 

0 

0 

0 

• 

• 

• 

• 

■ 

0 

0 

0 

0 

0 

B 

CL^ 

f^ 

CM 

CO 

CO 

<T\ 

LD 

CO 

r— 

uo 

ID  0. 

001 

0 

CO 

CM 

0 

00  X 

CM  E 

0 

0 

0 

0 

0 

0 OJ 

• ^ 

• 

• 

• 

• 

• 

O'— ' 

0 

0 

0 

0 

0 

(J  ^ 

0 

CO 

C^. 

00 

VD 

(O  E 1 

0 

CO 

LD 

CO 

U Q-  E 

0 

r— • 

0 

r— 

0 

«C  O- 

• 

• 

• 

• 

• 

0 

0 

0 

0 

0 

. 1 

r-“ 

LO 

1 — • 

00 

CO 

LD 

CO 

E • 

0 

CO 

LD 

CO 

Q-  E 

0 

r— 

0 

r“ 

0 

• 

• 

• 

• 

0 

0 

0 

0 

0 KJ- 

m 

00 

CM 

LIO 

CO 

1 ^ 

CM 

CTi 

LIO 

— V 

UO 

CO 

0 

CTi 

r— 

• 

• 

• 

• 

• 

p 1 

0 

LO 

0 

CO 

B 

cr» 

CO 

LIO 

u 

uo 

uo 

LIO 

LIO 

LIO 

CM 

CO 

CO 

CO 

CO 

N 

0 

r— 

VO 

CO 

r— 

Q- 

D_ 

CL 

CL 

CL 

CM 

CO 

CO 

T 

CO 

CO 

CO 

CO 

CO 

. *|<  I ilW  ■*  t 


(b)  From  Deaton,  et  al . 
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where  k is  the  absorption  coefficient. 


is  the  half-width  of  the  absorption  line,  Pq  = 1 atm,  Av  is  the 
difference  in  the  frequencies  of  the  laser  line  and  the  absorption 
line,  P is  the  total  pressure,  and  S is  che  line  strength.  For  a 
constant  temperature  aiid  P in  atmospheres,  a becomes 


(114)  a = a P 
0 


On  a curve  of  k versus  total  pressure  P,  the  condition  for 
zero  slope  with  Av  constant  is 


(115) 


or. 


^k 

da 


S 

TT 


(Av)^  + - 2g^ 

((Av)^  + g^)^ 


= 0 


(116)  g = Av 


Therefore,  as  a first  approximation  (i.e.,  assuming  a single 
isolated  Lorentz  line)  when  the  slope  of  the  k versus  P curve  is 
zero 


(117)  g^P  = Av  . 


Thus,  if  the  frequency  difference  between  the  laser  line  and 
the  absorption  line  is  known,  the  half-width  can  be  determined. 
Similarly,  if  the  half-width  is  known  the  frequency  difference  can 
be  determined. 


Now  if  g P = Av  is  substituted  into  Eq.  (112) 
is  obtained 


the  following 


(118) 


k 


Av 

* 2fsv)^ 


or 

(119)  S = 2TTk(Av) 
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Thus  the  line  strength  my  also  be  determined. 

When  this  experiment  was  performed,  the  laser  frequencies  had 
not  yet  been  accurately  measured.  Therefore  the  line  widths  and 
frequencies  listed  for  the  N2O  absorption  lines  on  the  AFCRL  tape[3] 
were  assumed  to  be  correct  and  this  technique  was  used  to  better 
determine  the  laser  line  frequencies. 

It  should  be  mentioned  that  Eq.  (117)  is  strictly  true  only 
if  the  pressure  is  high  enough  that  the  Voigt  profile  does  not  have 
to  be  used  to  describe  the  absorption  line  shape, ithere  is  no  pressure 
shift  in  the  frequency  of  the  absorption  line,  and  there  are  no  other 
absorption  lines  close  enough  to  contribute  to  the  absorption  at  the 
laser  frequency.  In  the  present  study,  the  first  assumption  is  valid, 
the  second  one  is  probably  valid,  but  the  third  assumption  is  not 
valid. 


A first  approximation  of  the  effect  of  additional  absorption 
lines  on  Eq,  (117)  car  be  obtained  by  including  the  nearest  absorp- 
tion line  on  either  side  of  the  nearly  coincident  line.  In  this 
case  Eq.  (112)  becomes 


(120) 
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where  the  non-subscripted  terms  refer  to  the  nearly  coincident  ab- 
sorption line  and  the  subscripts  1 and  2 refer  to  the  lines  on  either 
side  of  the  coincident  line. 


Assume  the  half-widths  and  line  strengths  are  all  about  the 
same  and  Av]  and  Av2  are  about  10  times  Av.  This  is  the  case  for 
the  N2O  band  under  consideration  here.  Equation  (120)  then  becomes 


(121) 


k = 5 
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Differentiating  Eq.  (118)  with  respect  to  a gives 


(122)  ^=1 

d a IT 


A 2 2 

Av  - a 


A 2 . 2 

Av.|  + a 


-A  2 2 

Av£  “ cx 


((Av)^  + a^)^  ((Ay^)^  + a^)^  ( (Av2)^  + ct^)^ 


Now  assume  that  the  two  additional  lines  cause  only  a small 
change  in  the  single  line  solution  which  is  a = Av.  Substitute  a = 

Av  in  the  denominator  of  the  first  term  of  Eq.  (122)  and  drop  a in 
the  second  and  third  terms  of  Eq.  (122)  since  (Av])2  and  (Av^)^  are 
both  much  greater  than  a2,  and  substitute  Av]  = Av2  - 100  Av.^Equation 
(122)  then  becomes 
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(123) 
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or 
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Solving  for  a‘~ 


= 0 


(126)  = ] .08  Av^ 


Substituting  Eq.  (126)  into  Eq.(122)  and 

(127)  1.0865  Av^ 


solving  for  a again  yields 


Repeating  the  process  one  more  time  gives 

(128)  a^=  1.087  Av^ 


or 


(129)  agP=1.0-f3Av 


coeff1c?ent“ve)sSs'to"?^^*;!ess'l!;e*I;a^'lls^^^  ^bsorptioo 

Of  the  laser  frequency.  Curves  nf  ^ first  estimate 

total  pressure  were  calrniatoH  fn  orption  coefficient  versus 

to  the  estimated  frequency  using  the^cnmnut^*^®'"  ^^^R^^i^cies  close 
perla^ota,  data  „as 
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Figure  62  ?hows  the  experimental  data  for  the  3-2  P(7)  line 
with  N2O  pressure  constant  at  6.51  x 10“^  torr  and  total  pressure 
varying  from  50  - 850  torr  along  with  calculated  curves  for  the 
same  N2O  pressure  and  laser  frequencies  2570.505  cm"l , 2570.510 
cm"’,  and  2570.515  cm"l.  The  laser  frequency  determined  from  Fig. 

62  is  2570.51  + .01  cm“ ' . This  compares  with  2570.522  + .003 
cm“l  determined  by  Heath,  et  al.[17]. 

Figures  63  and  64  show  the  experimental  data  for  the  2-1  P(10) 
line  with  N2O  pressures  6.82  x 10"^  torr  and  1.73  x 10"^  torr  re- 
spectively and  total  pressures  from  50  torr  to  760  torr.  Also  shown 
are  calculated  curves  for  the  same  N?0  pressures  and  laser  fre- 
quencies 2580.100  cm"l . 2580.105  cm~l,  and  2580.110  cm"'.  The  laser 
frequency  determined  from  the  plots  is  2580.10  + .005  cm"^ . This 
compares  with  2580.097  + .003  cm~l  determined  by  Heath,  et  al,[17]. 

For  the  2-1  P(10)  line  the  frequency  determined  results  in 
oqP  = 1.043  Av  at  the  zero  slope  point.  This  is  in  exact  agree- 
ment with  Eq.  (129).  This  is  probably  coincidental  considering  the 
assumptions  which  were  made  in  the  derivation.  For  the  3-2  P{7) 
line  at  the  zero  slope  point  oqP  = 1.068  Av  which  is  a little  higher 
than  predicted. 

The  technique  described  here  is  certainly  useful 
provided  proper  precautions  are  taken  in  interpreting  the  results. 


B.  CH^  Absorption  Measurements, 
760  Torr,  23^ 


Measurements  of  the  absorption  by  methane  samples  broadened  to 
one  atmosphere  total  pressure  by  dry  air  were  made  for  the  2-1  P(6), 
P(7)  and  P(8)  OF  laser  lines.  The  White  cell  was  set  at  0.7317  km 
for  these  measurements. 

As  with  the  nitrous  oxide  measurements  the  methane  concentra- 
tion had  to  be  increased  substantially  over  the  normal  atmospheric 
abundance  in  order  to  have  sufficient  absorption  so  that  it  could 
be  measured  accurately  in  a 0.7317  km  path.  The  highest  partial 
pressure  of  methane  used  was  1.72  torr  which  corresponds  to  2263 
parts  per  million.  The  normally  assumed  atmospheric  concentration 
is  1-6  parts  per  million.  Even  though  the  concentration  is  greatly 
enhanced  over  the  normal  concentration  the  methane  partial  pressure 
is  still  low  enough  that  the  results  can  be  linearly  extrapolated 
without  distortion  caused  by  self-broadening. 
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Fig.  63,  Absorption  coefficient  vs  total  pressure  for 
2-1  P(10)  line  at  6,82  x 10“3  torr  N2O, 
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Assuming  the  calculations  presented  in  Chapter  II  are  in  error 
by  no  more  than  'lOOi^j,  the  absorption  due  to  methane  for  the  six 
high-power  lines  in  the  2-i  and  3-2  bands  is  negligible  compared  with 
the  absorption  due  to  HDO  and  the  water  continuum.  One  of  the  pur- 
poses of  these  measurements  then  was  to  verify  that  the  methane 
absorption  was  indeed  small.  The  second  purpose  of  these  measure- 
ments was  to  provide  calibration  data  for  a spectrophone  experiment 
f/hich  has  been  reported  elsewhere[28]. 

The  procedures  and  experimental  setup  used  for  making  the 
methane  measurements  were  the  same  as  those  used  in  making  the  nitrous 
oxide  measurements  except  that  all  three  laser  lines  were  studied  in 
one  experimental  run. 

Figures  65,  66,  and  67  show  the  measured  absorption  coeffi- 
cients versus  methane  concentration  for  the  2-1  P(6),  P(7),  and 
P(8)  lines  respectively.  For  each  laser  line  a least  squares  fit 
of  the  absorption  data  to  a straight  line  through  the  origin  was 
made. 


Assuming  a natural  CH4  abundance  of  1.6  ppm[l6],  Table  7 
gives  the  CH4  absorption  coefficients  for  the  three  laser  lines  as 
extrapolated  from  the  experimental  data.  For  comparison  the  mea- 
surements reported  by  Spencer,  et  al.[24]  and  Deaton,  et  al.[26] 
are  also  presented  along  the  the  calculated  values  from  Chapter  II. 


TABLE  7 . 

METHANE  ABSORPTION  COEFFICIENTS  ASSUMING 
1 .6  ppm  CH^  IN  AIR  AT  SEA  LEVEL 


1 ine 

V 

(cm"^) 

k 

(km"^ ) 
calculated 

k 

( km"l ) 
OSU  exp. 

k 

( km"l ) 
Deaton  (a) 

k 

(km"^ ) 
Spencer  (b) 

2-1  P(6) 
2-1  P(7) 
2-1  P(8) 

2680.179 
2655.863 
2631 .068 

3.061x10"^ 

7.147x10"^ 

8.458x10"^ 

15.2x10"^ 
11 .3x10"^ 
8.59x10"^ 

11 .0x10"^ 
n .4x10"^ 
9.79x10"^ 

19.8x10"^ 

12.9x10"^ 

8.88x10"^ 

(a)  From  Deaton,  et.al.[25] 

(b)  From  Spencer,  et.al.[24] 
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CH^  CONCENTRATION  (ppm  X 10^) 


Fig,  65.  Measured  CH4  - Air  absorption  coefficient 
for  2-1  P(6)  line  at  2680,179  ctn'< . 


ABSORPTION  COEFFICIENT  (knrr>) 
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ABSORPTION  COEFFICIENT  (km-«) 


CH^  CONCENTRATION  (ppm  xIO^) 


Fig«  67,  Measured  CH4  - Air  absorption  coefficient 
for  2-1  P(8)  line  at  2631.068  cm'l . 


120 


The  r asured  absorption  coefficients  are  in  reasonable  agree- 
ment with  the  calculated  absorption  coefficients  only  on  the  2-1 
P(8)  line.  For  the  other  two  lines  the  agreement  is  njt  good  at  all. 
According  to  Benedict[29]  the  methane  data  on  the  AFCRL  line  data 
tape[3]  may  be  in  error.  Therefore  the  calculations  made  using  this 
data  would  also  be  in  error. 

The  CH4  measurements  reported  by  Deaton  were  made  using  the 
same  apparatus  and  procedures  as  for  the  N2O  measurements.  The^same 
calibration  factor  was  used  and  the  measurements  were  made  at  630 
torr  in  both  cases.  Considering  those  factors  the  agreement  of  the 
results  of  this  study  with  those  reported  by  Deaton  is  as  good  as 
could  be  expected. 

The  measurements  reported  by  Spencer  were  made  using  a 200  cm 
cell  and  pure  methane.  The  results  given  in  Table  7 were  obtained 
by  extrapolating  his  measurements  at  760  torr  to  1.6  ppm  assuming  a 
self-broadening  coefficient  of  1.  The  agreement  of  the  results  of 
this  study  with  those  reported  by  Spencer  is  reasonably  good  except 
for  the  2-1  P;6)  line.  On  this  line  some  error  could  be  introduced 
in  extrapolating  Spencer's  measurement  to  1.6  ppm. 

The  results  of  this  study  along  with  the  results  reported  by 
Deaton,  et.al . and  Spencer,  et.al.  confirm  that  the  methane  absorp- 
tion is  indeed  quite  small  compared  to  absorption  due  to  HDD  or  the 
water  vapor  continuum. 


C.  CO2  Absorption  Measurements 


The  calculations  presented  in  Chapter  II  indicate  CO?  absorp- 
tion coefficients  for  the  six  laser  lines  of  interest  in  the  3-2 
and  2-1  bands  of  about  lO""  or  lO'^  km'^  Myers[30]  has  investi- 
gated the  CO2  absorption  experimentally  using  a DF  laser  and  has 
found  that  the  absorption  coefficient  for  some  lines  is  four  orders 
of  magnitude  higher  than  the  calculations  predict.  This  absorption 
is  apparently  due  to  an  isotopic  or  weak  CO?  band  which  was  not 
included  in  the  AFCRL  compilation. 

Absorption  of  the  2-1  P(8)  DF  laser  line  by  pure  CO2  was 
measured  for  CO?  pressures  of  248,  503,  and  761  torr  in  a path  of 
0.7317  km.  The  optical  setup  and  detectors  were  the  same  as  for 
the  N?0  and  CH4  measurements.  Myers  found  that  this  line  had  higher 
absorption  than  any  of  the  other  six  lines  of  interest  in  this  study 
The  results  of  the  measurement  are  shown  in  Fig.  68.  The  results 
of  this  experiment  are  probably  less  reliable  than  the  results  of 
the  N?0  and  CH4  measurements  because  the  detector  calibration  was 
not  sjfficiently  well  known.  The  point  at  761  torr  is  probably  with 
in  1C%  with  the  other  two  points  being  somewhat  worse. 
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ABSORPTION  COEFFICIENT  (km 


Measur^u  CO2  absorption  v/OottiCiCnt  tor 
2-1  P(8)  line  at  2631,068  cnr' . 


The  measured  absorption  coefficient  was  1.6  km  at  761  torr . 
Assuming  a self-broadening  coefficient  of  1 for  CO2  this  corresponds 
to  5.3  X 10"^  km"'  at  330  ppm  CO2  and  1 atmo'  phere  total  pressure. 
Assuming  a self-broadening  coefficient  of  1.3  and  using  the  data  in 
Fig.  68,  a value  of  4.7  x 10"^  km*'  at  330  ppm  and  one  atmosphere  is 
obtained.  Table  8 lists  the  absorption  coefficients  obtained  by 
Myers  for  the  six  high  power  laser  lines  along  with  the  result  ob- 
tained in  this  study  for  the  2-1  P(8)  line. 


TABLE  8 

MEASURED  CO2  ABSORPTION  COEFFICIENTS  FOR  6 OF  LASER  LINFS 

I_'jn0  '^1  (km  ~ ) k (km  ) 

(cm“  ) Myers  (a)  this  study 


2-1  P(6) 

2680.179 

.19  X 10"^ 

2-1  P(7) 

2655.863 

1 .85  X 10"^ 

2-1  P(8) 

2631 .068 

9.26  X 10"^ 

3-2  P(6) 

2594.198 

6.85  X 10’^ 

3-2  P(7) 

2570.522 

5.24  X 10"^ 

3-2  P(S) 

2546.375 

4.21  x 10"^ 

(a)  From  MyersCaoJ 

The  difference  between  the  absorption  coefficient  for  the  2-1 
P(8)  line  measured  in  this  study  and  that  measured  by  Myers  is  not 
completely  accounted  for.  Myers  has  indicated  that  the  numbers  he 
has  quoted  are  not  as  accurate  as  would  be  indicated  by  the  number 
of  significant  figures.  In  any  case  the  CO2  absorptioi  coefficients 
are  small  compared  to  the  HDO  and  water  continuum  absorption  coeffi- 
cients. 


D . HDO  Absorption  Measurements, 

760  torr,  24°C 

Absorption  of  the  2-1  P(6),  P(7)  and  P(8)  and  the  3-2  P(6), 
P(7)  and  P(8)  OF  laser  lines  by  HDO-nitrogen  mixtures  was  measured. 
The  detectors,  optical  arrangement,  and  data  recording  procedures 
were  the  same  for  these  measurements  as  for  the  N2O  absorption 
measurements. 


123 


As  with  the  NoO  and  CH/j  absorption,  the  HDO  absorption  for  the 
six  laser  lines  of  interest  in  this  study  is  too  low  in  normal  atmo- 
spheric samples  to  be  measured  in  a 0.7317  km  path.  In  addition, 
the  calculations  presented  in  Chapter  II  indicate  that  the  HDO  ab- 
sorption is  the  same  order  of  magnitude  as  the  water-vapor  continuum 
absorption.  Therefore  it  was  necessary  to  enhance  the  HDO  concentra- 
tion relative  to  the  H2O  concentration  for  two  reasons.  First,  the 
HDO  enhancement  was  necessary  in  order  to  obtain  high  enough  absorp- 
tion to  accurately  measure  in  a 0.7317  km  path.  Second  it  was 
necessary  to  enhance  the  HDO  concentration  so  that  the  HDO  absorption 
was  much  greater  than  the  water  vapor  continuum  absorption.  The  small 
continuum  absorption  would  not  then  cause  any  appreciable  error  in 
the  measurement  of  HDO  absorption. 

In  these  experiments  it  was  necessary  to  devise  a careful  sam- 
ple preparation  procedure. 

In  the  steady  state,  HDO  does  not  exist  alone  but  rather  in 
equilibrium  with  H2O  and  D2O  according  to  the  following  equation: 


(130)  HoO  + DoO  ^ ^ 2HD0 

The  relative  concentrations  of  the  three  types  of  water  can  be  cal- 
culated from  the  equilibrium  constants,  K292°  = 3.543  and  K293 
3.506[24]  for  the  liquid  and  gas  phases  respectively. 

The  procedure  used  in  preparing  the  samples  for  the  work  pre- 
sented here  was  to  mix  a known  amount  of  D2O  with  a known  amount  of 
H2O  and  calculate  the  resulting  amount  of  HDO  using  the  equilibrium 
constant.  A formula  relating  the  final  amount  of  HDO  to  the  initial 
amounts  of  D2O  and  H2O  can  be  derived  as  follows: 

For  a mixture  of  H2O,  D2O,  and  HDO  where  the  molal  concentra- 
tions of  the  three  types  of  water  are  A,  B,  and  C respectively,  the 
equilibrium  concentrations  of  each  of  the  three  types  of  water  are 
related  by 


(131) 


AB 


= K 


where  the  appropriate  K for  the  liquid  or  gas  phase  is  used, 


Let  0^20  initial  volumes  of  H2O,  CQ2O  be  the  initial 

volume  of  D2O  and  Chdo  be  the  final  volume  of  HDO.  The  final  con- 
centrations of  H2O  and  D2O  are  then  (C^20  " ^Hgo^^^  (^D20  " 

Cunn/2)  respectively.  Substituting  these  definitions  for  A,  B,  and 
C 1 n Eq.  (131)  gives 
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(132) 


or 


(133) 


^H20  ^D20 


(C  + C )+ 

2 ^^H20  ^D20^  4 


= K 


Equation  (133)  may  be  written  as  follows: 


(134)  (4-K)  + 2K  (Cp^  n ■*■  n)  - 4K  Cu  ^ C 


HDO  '"D2O  ^H20^ 


'H2O  "D2O  = 0 


Equation  (134)  may  now  be  solved  for  C^qq  using  the  quadratic  formula: 


'1:3  5) 


'HDO 


'WkJ 


H20  ''D20 


Three  separate  mixtures  containing  different  HDO  concentrations 
were  used.  The  appropriate  mixture  was  chosen  for  each  laser  line  so 
that  the  absorption  could  be  measured  accurately  in  a 0.7317  km  path. 
The  samples  used  were;  for  the  2-1  P(6)  line,  .01%  D2O,  1.99%  HDO, 
balance  H2O;  for  the  2-1  P(7)  line,  .003%  D2O,  1.02%  HDO,  balance 
HoO;  and  for  the  2-1  P(8),  3-2  P(6),  3-2  P(7),  and  3-2  P(8)  lines, 

.T2%  D2O,  6.22%  HDO,  balance  H2O.  The  percentages  were  calculated 
using  the  gas  phase  equilibrium  constant  K = 3.506. 

There  is  a potential  problem  which  could  be  encountered  when 
the  enhanced  sample  is  introduced  into  the  White  cell.  This  is  caused 
by  the  fact  that  HDO  is  less  volatile  than  H2O  and  therefore  remains 
preferentially  in  the  condensed  phase,  It  has  been  reported[31]  that 
in  thermodynamic  equilibrium,  water  vapor  over  liquid  water  will  con- 
tain 8%  fewer  HDO  moleucles  than  the  liquid  water.  When  the  water 
vapor  is  being  introduced  into  the  White  cell,  the  liquid  and  gas 
are  certainly  not  in  thermodynamic  equilibrium,  so  the  HDO  concen- 
tration in  the  cell  could  be  quite  uncertain.  This  problem  was  over- 
come by  filling  a small  bottle  with  just  enough  of  the  specially 
prepared  water  to  fill  the  White  cell  with  the  desired  amount  of 
water  vapor  and  evaporating  the  sample  completely. 

Another  problem  which  could  affect  the  HDO  concentration  would 
be  preferential  adsorption  of  one  water  isotope  relative  to  the  other. 
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The  cell  has  been  tested  several  times  for  HoO  adsorption  and  none 
has  been  observed.  A mass  spectrometer  residual  gas  analyzer  was 
added  to  the  cell  so  that  the  HDO/H2O  ratio  could  be  monitored  as 
a function  of  time  in  order  to  evaluate  this  problem.  The  instru- 
ment was  not  accurate  enough  for  absolute  concentration  determi- 
nation. A water  vapor  sample  containing  an  enhanced  concentration 
of  HDO  was  admitted  to  the  cell  and  the  HDO/H2O  ratio  was  monitored 
using  the  mass  spectrometer  over  a period  of  several  hours.  Within 
the  accuracy  limits  of  the  instrument,  no  change  in  the  ratio  was 
observed.  A more  complete  discussion  of  this  experiment  is  given 
in  Appendix  D. 

Figures  69  through  74  show  the  results  of  the  HDO  absorption 

measurements  for  the  2-1  P(6),  P(7),  and  P(8)  lines  and  the  3-2  P(6), 

P(7),  and  P(8)  lines  respectively.  The  absorption  cell  path  length 
was  0.7317  km.  The  cell  originally  contained  15  torr  of  the  appro- 
priate special  water  sample  plus  nitrogen  to  a total  pressure  of 
760  torr.  The  lower  points  on  the  curves  were  obtained  by  partially 
pumping  the  cell  out  and  refilling  to  760  torr  with  nitrogen. 

For  each  line  a least  squares  fit  of  the  data  to  an  expression 
of  the  form  k = AP+BP*^  was  made,  where  k is  the  absorption  coeffici- 
ent and  p is  the  partial  pressure  of  enriched  wati^r  in  torr.  The  de- 
rived expression  for  a laser  line  is  presumed  to  be  a more  accurate 

characterization  of  the  absorption  coefficient  than  any  individual 
data  point.  The  expressions  for  the  measured  absorption  coefficients 
for  the  three  lines  are  presented  in  Table  9.  The  expressions  are 
valid  for  760  torr  total  pressure  with  p being  the  partial  pressure 
of  the  enriched  water. 

In  Table  10  the  expressions  from  Table  9 are  extrapolated  to 
an  assumed  abundance  of  HDO  to  H2O  of  0.03%.  This  was  done  by  multi- 
plying the  coefficients  A and  B by  ,03 /x  where  x is  the  percent  HDO 
concentration  of  the  mixture  used  to  make  the  measurements.  This 
procedure  assumes  that  line  broadening  caused  by  HDO  - HDO  collisions 
is  not  greatly  different  from  line  broadening  due  to  HDO  - H2O  col- 
lisions. These  expressions  have  been  evaluated  at  14.26  torr  for 
each  line  and  the  values  listed  in  Table  10  along  with  the  absorption 
coefficients  calculated  in  Chapter  II.  It  should  be  pointed  out  that 
the  expressions  for  absorption  coefficients  in  Table  10  are  only 
valid  at  24°C  and  760  torr  total  pressure.  The  temperature  depend- 
ence and  total  pressure  dependence  of  the  absorption  coefficients 
have  not  been  investigated. 

Since  the  samples  used  to  make  the  measurements  also  contained 
an  enhanced  D2O  concentration  as  well  as  enhanced  HDO,  it  is  possible 
that  some  error  might  be  introduced  due  to  D2O  absorption.  This 
possibility  was  investigated  by  measuring  the  absorption  of  a sample 
of  pure  D2O  vapor  broadened  with  dry  air.  The  D2O  vapor  pressure 
was  0.21  torr  which  is  12-500  times  the  D2O  partial  pressure  at  the 
highest  water  vapor  pressure  for  the  HDO  measurements  in  Figs.  69 
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TABLE  9 

MEASURED  HDD  ABSORPTION  COEFFICIENTS  FOR  SIX  OF  LASER  LINES 


concentration 


TABLE  10 


HDO  ABSORPTION  COEFFICIENTS  EXTRAPOLATED  TO  .03% 
RELATIVE  HDO  ABUNDANCE 


K (km  ')  .03%  HDO  K (km  ') 

760  torr  total  pressure  14.26  torr  fUO 


2-1  P(6) 
2-1  P(7) 

2- 1  P(8) 

3- 2  P(6) 
3-2  P(7) 
3-2  P(8) 


3.39x10  p + 2.05x10  ’'P 
-3  -5 

I-  t r oyi..T/-k  ^r\ 


5.24x10 

3.79x10 

9.54x10"^ 

7.35x10 

6.00x10"^ 

9.12x10' 

1 .74x10'^ 

7.18x10 

rV 

8.61x10'^ 

4.53x10 

rV 

2.46x10'^ 

1 .13x10 

ABSORPTION  COEFFICIENT 


2-1  P(8) 

P = 760  torr 


PARTIAL  PRESSURE  (torr  ) 


Fig.  71,  Measured  HDO-N2  absorption  coefficient  for 
2-1  P(8)  line  at  2631.068  cm-1 . 


H 

i 


ABSORPTION  COEFFICIENT  (fern” 


Fig,  72,  Measured  HDO-N2  absorption  coefficient  for 
3-2  P(6)  line  at  2594,198  cm'l . 


131 


6.2  PARTS  HDO 


PARTS  TOTAL  WATER 


PARTIAL  PRESSURE  I torr  ) 


Measured  HDO-N2  absorption  coefficient  for 
3-2  P(8)  line  at  2546,375  cm‘l . 


to  '4.  There  was  measurable  D2O  absorption,  however  if  the  measure- 
ments are  extrapolated  to  the  conditions  prevailing  in  the  HDO  ab- 
sorption measurements,  the  D2O  contribution  to  the  absorption  coef- 
ficients at  the  highest  water  vapor  pressure  is  .02  or  .03km~l  which 
is  less  than  the  unce^-tainty  in  the  experimental  data. 


E.  H2O  Absorption  Measurements 


Absorption  of  the  2-1  P(6),  P(7),  and  P(8)  and  3-2  P(6),  P(7), 
and  P(8)  DF  laser  lines  by  H20-air  and  H20-nitrogen  mixtures  was 
measured.  It  was  originally  thought  that  water  vapor  absorption  would 
be  too  small  to  measure  accurately  in  the  White  cell  at  room  tempera- 
ture. However  preliminary  spectrophone  measurements [28]  indicated 
higher  absorption  on  some  lines  than  had  bean  expected.  It  was  there- 
fore decided  that  with  a very  carefully  performed  experiment  it  might 
be  possible  to  measure  the  absorption  directly  using  the  White  cell. 

Improvements  in  the  White  cell  optics  which  were  discussed 
in  Chapter  II  made  it  possible  to  increase  the  path  length  to  1.341 
kilometers  with  only  a small  increase  in  insertion  loss.  Also  a new 
pulsed  laser  was  constructed  (described  in  Chapter  III)  which  was 
primarily  intended  to  eliminate  RF  noise  emission.  While  testing  the 
new  laser  it  was  also  found  that  the  average  power  was  improved  signi- 
ficantly over  the  earlier  model.  This  made  possible  the  use  of  thermo- 
pile detectors  which  eliminated  the  need  for  detector  linearity  cor- 
rection which  had  been  a source  of  error  with  the  measurements  using 
the  lead  selenide  detectors.  The  thermopiles  were  carefully  shielded 
in  styrofoam  to  reduce  the  effect  of  air  currents  and  short  term 
fluctuations  in  room  temperature. 

In  nrevious  measurements  the  ratio  of  input  power  to  output  povi/er 
for  the  evacuated  cell  and  the  cell  filled  with  the  desired  sample  had 
been  determined  cy  making  500  separate  measurements  of  the  ratio  over 
about  a one  minute  period  and  using  the  average  as  the  final  value. 

The  sample  transmittance  was  then  just  the  quotient  of  the  ratio 
obtained  with  the  sample  in  the  cell  and  the  ratio  obtained  with  the 
cell  evacuated.  The  absorption  coefficient  was  then  obtained  by 
dividing  the  negative  of  the  natural  logarithm  of  the  transmittance 
by  the  path  length. 

For  the  measurements  of  HDO,  N2O,  and  CH4  absorption  it  was 
possible  to  measure  the  absorption  coefficient  at  several  absorber 
concentrations  and  use  least-squares  curve-fitting  techniques  to 
reduce  the  uncertainty  of  the  measurements.  This  was  not  possible 
for  the  H2O  measurements  however  since  the  absorption  ac  the  highest 
water  vapor  pressure  was  still  quite  low. 

The  procedure  used  therefore  was  to  repeat  the  measurement 
of  input  to  output  ratio  described  above  several  times  in  a period 
of  a half-hour.  These  ratios  were  then  averaged  to  obtain  a final 
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value.  Also  the  ratio  of  input  power  to  output  power  for  the  evacu- 
ated cell  was  measured  both  before  and  after  the  measurement  of  the 
ratio  for  the  sample.  Transmittances  could  then  be  calculated  using 
each  empty  cell  ratio  separately  or  the  average  of  the  ratios.  The 
entire  procedure  was  repeated  two  or  three  times  for  each  line,  i.e., 
on  different  days  and  using  different  water  samples. 

The  optical  setup  shown  in  Fig.  44  was  used  for  these  measure- 
ments. Measurements  were  made  with  a path  length  of  1.341  km,  water 
vapor  pressure  of  14.3  torr  and  total  pressure  including  dry  air 
or  nitrogen  of  760  torr.  Both  dry  air  and  nitrogen  were  used  as 
broadening  gases,  and  within  error  limits  no  difference  in  absorption 
coefficient  was  observed.  Also  there  was  no  observable  difference 
between  the  transmittance  of  the  evacuated  cell  and  the  transmitt- 
ance of  the  cell  when  filled  with  either  dry  air  or  nitrogen. 

The  water  vapor  sample  was  introduced  into  the  cell  from  a 
bottle  containing  distilled  water  attached  to  the  sample-mixing  mani- 
fold. when  the  water  bottle  was  first  attached  to  the  manifold  it 
was  pumped  on  for  several  hours  to  remove  all  dissolved  gas.  The 
water  bottle  was  heated  about  5°C  above  room  temperature  in  order 
to  speed  the  filling  process.  About  hours  was  required  to  fill 
the  cell  with  14,3  torr  of  water  vapor.  The  water  vapor  pressure 
was  determined  approximately  using  the  mercury  mi cromanometer  as 
the  water  vapor  was  introduced.  When  the  desired  amount  of  water 
vapor  had  been  introduced,  the  cell  was  filled  to  760  torr  with 
either  nitrogen  or  dry  air  and  the  fans  at  each  end  of  the  cell  were 
turned  on  and  the  sample  was  allowed  to  mix  over  night.  The  actual 
water  vapor  pressure  was  determined  using  the  E.G.  & G.  Model  992 

Dew  Point  Hygrometer  at  the  time  the  transmittance  measurement  was 
made . 


in  Tables  11,  12,  13,  14,  15,  and  16 
for  the  2-(  P(6),  P(7),  and  P(8)  and  3-2  P(6),  P(7),  and  P(8)  DF 
laser  lines.  The  ratio  of  input  to  output  measured  with  the  sample 
in  the  cell  is  designated  R ; . The  ratio  measured  for  the  evacuated 
cell  before  the  sample  is  admitted  is  BK--,  and  the  ratio  measured 
when  the  cell  is  evacuated  after  measuring  the  sample  is  designated 
BKi2 


J.U  n-}P  average  values  for  each  line  are  compared  with 

the  HJO  measurements  presented  in  Table  10  calculated  H2O  values 
from  Table  4,  and  H2O  continuum  values  obtained  using  Burch's  data[4]. 

On  the  2-1  P(6)  line,  the  calculated  HoO  absorption  coefficient 
IS  very  small,  and  the  measured  HDO  absorption  coefficient  plus  the 
water  continuum  absorption  coefficient  determined  from  Burch's  mea- 
surements agree  quite  well  with  the  measured  H2O  absorption  coeffici- 
ent. Agreement  is  not  as  good  on  the  2-1  P(7)  line  although  the 
measured  HDO  absorption  coefficient  added  to  water  continuum  -absorp- 
tion coefficient  obtained  from  Burch's  measurements  is  certainly 
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within  the  experimental  scatter  shown  in  Table  12,  On  the  2-1  P(8) 
line  agreement  is  again  only  fair,  although  here  the  difference  might 
be  attributable  to  an  error  in  the  calculated  H2O  absorption  fcoeffi- 
cient.  For  the  3-2  P(6),  P(7),  and  P(8)  lines  the  measured  H2O 
absorption  coefficient  is  somewhat  greater  than  the  measured  HDO 
absorption  coefficient  added  to  the  water  continuum  absorption  co- 
efficient obtained  from  Burch's  measurements  (Column  4 in  Table  17) 
in  each  case.  For  all  three  lines  however  the  value  in  Column  4 of 
Table  17  is  within  the  experimental  scatter  shown  in  Tables  14,  15, 
and  16. 

One  conclusion  to  be  drawn  from  the  direct  measurements  of  H2O 
absorption  presented  here  is  that  the  water  continuum  absorption 
obtained  by  extrapolating  Burch's  high  temperature  measurements  is 
probably  correct.  A series  of  carefully  performed  spectrophone 
measurements  would  be  the  best  way  to  determine  the  water  continuum 
absorption  accurately.  This  was  beyond  the  scope  of  this  study  since 
a CW  DF  laser  was  not  available.  The  only  spectrophone  measurements 
of  the  water  continuum  performed  thus  far[28]  did  not  have  the  re- 
quired accuracy  due  to  experimental  difficulties. 

F.  Summary  of  Molecular  Absorption  for  the  Eight  Lines  Studied 

Table  18  gives  the  contribution  to  the  molecular  absorption  due 
to  each  absorber  at  each  of  the  8 laser  lines  investigated  for  the 
sea-level  mid-latitude  summer  model  whose  parameters  are  given  in 
Table  3.  The  values  listed  are  believed  to  be  the  most  accurate 
available  at  this  time.  For  N2O,  CH4,  and  HDO  the  values  given  are 
those  measured  in  this  study  for  those  lines  where  they  are  available. 
For  CH4,  the  absorption  coefficients  for  the  remaining  lines  are  those 
obtained  bv  Deaton,  et  al . [25]  or  the  calculated  values  from  Table  4. 
For  HDO,  the  absorption  coefficients  for  the  lines  not  measured  in 
this  study  are  taken  from  the  calculations  in  Table  4.  For  COo.  the 
absorption  coefficients  are  those  obtained  by  Myers  [30 3 except  for 
the  2 - 1 P(8)  line  which  was  measured  in  this  study.  The  contri- 
butions due  to  local  water  lines  are  obtained  from  the  calculations 
1n  Table  4,  The  water  vapor  continuum  absorption  and  the  pressure 
induced  nitrogen  contributions  are  those  obtained  by  Burch  [4]. 

For  each  line  the  total  absorption  coefficient  due  to  molecular 
absorption  is  presented.  From  the  total  molecular  absorption  coef- 
ficients the  transmittance  of  a ten  kilometer  sea  level  path  was 
calculated  for  each  line.  These  values  are  also  presented  in  Table  18. 
The  transmittance  over  a ten  kilometer  sea  level  path  varies  from 
0.31  for  the  2 - 1 P(7)  line  to  0.73  for  the  2 - 1 P(8)  line.  It 
should  be  pointed  out  that  the  actual  transmittance  through  a real 
atmosphere  would  probably  be  lower  than  shown  in  Table  18  because  of 
aerosol  scattering  and  absorption. 
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TABLE  n 

H2O  ABSORPTION  ON  THE  2-1  P( 
LINE  AT  2680.179  cm 

6j  DF  LASER 

Rl/BK-,1 

R1/BK12 

R2/BK21 

R2/BK22 


1.071/1.183 

1.071/1.182 

1.071/1,181 

1.071/1.192 


T 

K (Km"'') 

0.905 

0.074 

0.906 

0.074 

0.907 

0.073 

0.898 

0.080 

= 1.071/1  1825  0.906 


0.074 


.071/1.1865  0.903  0. )7fi 


TABLE  12 


H2O  ABSORPTION  ON  THE  2-1  P{7)  OF  LASER 
LINE  AT  2655.863  cm-'' 


Avg. 

K (Km-1)  K (Km-l) 


Rl/BKii 

Ri/BKi2 

R2/BK21 

R2/BK22 

R3/BK31 

R3/BK32 


,/BKn  + BKii 

'A  ' . 

BK21  + BK22 


BKq]  + BK32  \ 

«3/l — r—J 


196 

0.882 

0.094 

204 

0.876 

0.099 

204 

0.858 

0.114 

211 

0.853 

0.119 

211 

0.863 

0.110 

170 

0.893 

0.084 

= 1.055/1.200 

0.879 

0.096 

= 1.033/1.208 

0.855 

0.116 

= 1.045/1.191 

0.878 

0.097 

0.103 
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TABLE  13 

H2O  ABSORPTION  ON  THE  2-1  P(8)  OF  LASER 
LINE  AT  2631.068  cm-1 


K (Kni-l) 


Avg. 

K (Km"'') 


= 1.422/1.499 
= 1.422/1.484 
= 1.414/1.484 
= 1.414/1.498 


0.949 

0.958 

0.953 

0.94^ 


0.039 

0.032 

0.036 

0.043 


0.0376 


BK21  + BK22 


= 1.414/1.491 


0.953 


0.948 


0.036 


0.040 


0.038 


TABLE  14 


H2O  ABSORPTION  ON  THE  3-2  P(6)  OF  LASER 
LINE  AT  2594.198  cm"' 


Avg. 

T K (Km-1)  K (Km-1) 


R]/BKii  = 1.151/1.230 

0.936 

0.050 

R2/BK21  = 1.182/1.336 

0.885 

0.091 

R2/BK22  = 1.182/1.329 

0.889 

0.087 

R3/BK31  = 1.248/1.329 

0.939 

0.047 

R3/BK32  = 1.248/1.329 

0.939 

0.047 

R4/BK41  = 1.296/1.357 

0.955 

0.034 

R4/BK42  = 1.296/1.392 

0.931 

0.053 

Rl/BKii  = 1.151/1.230 

0.936 

0.050 

If  BK2]  + BK22^\ 

R2  A 2 J = 1.182/1.333 

0.887 

0.089 

^3  ~j  = 1-248/1.329 

0.939 

0.047 

+ BK42^ 

y 2 J = 1.296/1.375 

0.943 

0.044 

0.059 


0.057 
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TABLE  15 


H2O  ABSORPTION  ON  THE  3-2  P{7)  OF  LASER 
LINE  AT  2570.522  cm-1 


T K { Km" ^ ) 


Avg. 

K (Km-l) 


1.324/1.389 

0.953 

0.036 

1.324/1.397 

0.948 

0.040 

1.357/1.428 

0.950 

0.038 

1.357/1.422 

0.954 

0.035 

1.357/1.422 

0.954 

0.035 

1.357/1.403 

0.967 

0.025 

1.355/1.417 

0.956 

0.033 

1.355/1.420 

0.954 

0.035 

+ BK^2^ 

2 y = 1.324/1.393 

0.950 

0.038 

+ BK22\ 

2 1 = 1.357/1.425 

0.952 

0.036 

- 1.357/1.413 


0.035 


0.961  0.030 


0.035 


■]  = 1. 


355/1.419  0.955 


0.034 
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TABLE  16 

H2O  ABSORPTION  ON  THE  3-2 
LINE  AT  2546.375 

P(8)  OF  LASER 
cm“^ 

T 

K (Km 

Rl/BKn  = 

1.352/1.422 

0.951 

C 

t 

1 

R2/BK21  = 

1.394/1.469 

0.949 

C 

R2/BK22  “ 

1.394/1.429 

0.976 

C 

R3/BK31  = 

1.383/1.429 

0.968 

C 

R3/BK32  = 

1.383/1 .432 

0.966 

( 

Rl/BK^l  = 

1.352/1.422 

0.951 

! 

+ BKppN 

"2 —J  = 1.394/1.449 

0.962 

+ BK„  \ 

—j  = 1.383/1.431 

0.967 

TABLE  17 


H2O  MEASUREMENTS  COMPARED  WITH  CALCULATED  HoO  VALUES 
MEASURED  HDD  VALUES.  AND  THE  H2O  CONTINUUM 


Line 

1 

H?0 

calc, 

(km-l) 

2 

HDO 

exp 

(km-M 

3 

H2O  (a) 
cont. 
(km-1) 

4 

H2+3 

(km-1) 

5 

H2O 

measured 

(km-') 

2-1  P(S) 

2.:4x10“4 

5.24x10“^ 

2.12x10-2 

7.38x10-2 

7.5x10-2 

2-1  P(7) 

4.42x10“® 

9.54x10-2 

1.93x10“2 

1.15xl0-1 

1.03x1 0-1 

2-1  P(8) 

4.73x10-3 

6.00x10“3 

1.78x10“^ 

2.85x10-2 

3.8x10-2 

3-2  P(6) 

3.54x10-4 

1.74x10-2 

1.68x10-2 

3.46x10“^ 

5.7x10-2 

3-2  P(7) 

4.69x10“® 

8.6x10-3 

1.74x10-2 

2.60x10-2 

3.5x10-2 

3-2  P(8) 

4.44x10-5 

2.46x10-3 

1.86x10-2 

2.11x10-2 

3.1x10“2 

(a)  From  Burch  [4] 
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SUMMARY  OF  MOLECULAR  ABSORPTION  COEFFICIENTS  FOR  8 DF  LASER  LINES 


niiii 


The  purpose  of  this  study  was  to  determine,  as  accurately  as 
possible,  the  absorption  of  DF  laser  radiation  by  atmospheric  qases. 

Computer  programs  were  written  which  can  be  used  to  calculate 
the  absorption  of  radiation  by  atmospheric  gases  from  the  AFCRL 
Line  Compilation  [3].  Using  these  programs,  absorption  coefficients 
for  all  the  DF  laser  lines  and  all  the  atmospheric  absorbers  were 
calculated.  These  calculated  absorption  coefficients  were  used  to 
plan  experiments  to  measure  accurately  the  absorption  coefficients 
for  selected  DF  laser  lines.  Lines  selected  for  study  were  the  2-1 
P(6),  P(7),  and  P(8)  lines  and  the  3-2  P(6),  P(7),  and  P(8)  lines. 
For  these  lines  absorption  coefficients  V'ere  measured  for  N2^>  CH4, 
CO2,  HDD,  and  H2O.  For  a given  laser  line,  the  absorption  coeffic- 
ient was  measured  for  those  gases  which  were  significant  absorbers 
at  that  frequency,  and  not  necessarily  for  all  gases.  In  addition 
to  the  lines  listed  above,  the  NoO  absorption  coefficient  was  mea- 
sured for  the  2-1  P(10)  and  P(ll)  lines  since  they  were  the  only 
other  DF  laser  lines  for  which  N2O  absorption  was  significant. 

The  results  of  the  measurements  were  compared  with  the  cal- 
culations and  with  results  obtained  by  other  workers.  It  was 
found  that  the  calculated  values  for  the  N2O  absorption  coefficients 
agreed  very  well  with  the  measured  absorption  coefficients.  For 
the  other  gases  the  agreement  between  calculation  and  measurement 
was  not  as  good,  although  in  most  cases  the  calculations  and  mea- 
surements were  of  the  same  order  of  magnitude.  For  CO2  the  calcula- 
ted absorption  coefficients  are  about  10"7  km~\  Myers  [30]  found 
that  the  CO2  absorption  coefficients  were  actually  about  10~^  km~l . 
The  CO2  absorption  coefficient  was  measured  for  the  2-1  P(8)  line 
and  found  to  be  of  the  same  order  of  magnitude  as  that  measured  by 
Myers.  The  discrepancy  in  the  calculations  is  probably  caused  by 
a weak  isotope  band  which  is  not  included  in  the  AFCRL  tapes.  The 
H2O  absorption  measurements  confirmed  that  the  values  for  water  con- 
tinuum absorption  obtained  by  extrapolating  the  high  temperature 
measurements  of  Burch,  et.al.  [4]  are  close  to  the  correct  values. 

The  transmittance  of  a ten  kilometer  sea  level  path  for  the  Mid 
Latitude  Summer  model  [16]  was  calculated  from  the  most  accurate 
values  of  the  molecular  absorption  coefficients  for  each  of  DF  laser 


varied «V5’^  study  The  calculated  transmittance 
varied  from  31%  to  83%  for  a ten  kilometer  path.  The  ac+ual  trans- 

imttance  over  a ten  M loraeter  path  in  the  r^al  atmosphere  would  be 
somewhat  lower  since  attenuation  due  to  turbulence  and  aerosol^cat- 
tenng  and  absorption  has  not  been  accounted  for. 

measu^^IlL^  should  include  an  extension  of  these 

^ measurements  over  a range  of 

total  pressure  and  temperature.  There  is  also  a need  for  carefullv 
performed  spectrophone  measurements  to  more  accurately  determine  the 
water  contiuum  absorption  and  the  pressure  induced  nitrogen  absorption, 


APPENDIX  A 

UNITS  AND  CONVERSION  FACTORS  USED  IN  MOLECULAR 
ABSORPTION  CALCULATION  PROGRAMS 


The  information  i i this  appendix  has,  appeared  in 
earlier  reports  [32,33].  It  is  presented  here  for  complete- 
ness and  easy  reference. 


A,  Basic  Equation  for  Lorentz  Line 


(136)  «,n  T = - ku 

where  T = transmittance 

u = absorber  concentration 
k = extinction  coefficient 


and 


(137) 


S-  a. 

k = y 

V 4 r/  \2  2-, 

1 tt[(v-v.  ) + a.  ] 


where  S^.  = line  strength 
V = wavenumber 
v^.  = line  center  wavenumber 
a,j  = half-width  at  half  intensity. 


B.  Units  of  S and  u 


1.  It  is  necessary  to  specify  the  units  of  S and  u so 
that  the  product  ku  is  dimensionless, 

It  has  been  common  to  use  different  units  for  CO2  and 
other  gases  and  for  water  vapor.  Recently,  however, 
Cal  fee  and  otheis  have  begun  to  use  similar  units  for 
all  gases.  All  of  these  units  will  be  discussed. 
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Units  of  S 


HpO:  S is  in  cnrVgm  cm"2 

^ (pr-cm  and  gm  cm”2  are  the  same) 

CO2:  S is  in  cm*l/atm-cm. 

H2O  and  CO2  (uniform  system):  S is  in  cm'Vmolecules 

cm“2. 

3.  Units  of  Absorber  Concentration,  u . 

-2 

H2O:  gm  cm  or  pr-cm 

288.34  X 10  ^Ph  o.torr  Vr 

(138)  p = pt  = j pr-cm 

2 0^ 


for  reference 


0-9)  Ph  0 ' + .00367 

2 0 


1.05821 

-= gm/cm'’ 


CO2'.  atm-cm^ 


^^002  " PCO2  " 


It  is  necessary  to  specify  the  temperature  since  pressure 
depends  on  temperature  through  the  gas  law  p = NkT, 
i.e.,  Pj.g  is  CO2  pressure  in  atm.  at  T^. 

«,  - path  length  in  cm. 

2 

H2O  and  CO2  (uniform  system):  molecules/cm 
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4.  Conversion  factors  between  separate 
and  uniform  system[34]. 

19  2 

(140)  u[atm-cmgyp]  x 2.689  x 10  = u[inolecules/cni  ] 

22  2 

(141)  u[pr  cm^  q]  x 3.34  x 10  = u[molecules/cm  ] 

-1  -20  -1  -2 

(142)  S^[cm  /atm-cm^jp]  x 3.72  x 10  = S^[cm  /molecules  cm 

(143)  S^[cm'Vpi^  cm]  x 2.991  x 10'^^  = S^Ccm" Vmolecules  cm"^]. 

It  is  most  convenient  to  have  one  computer  program  for  all 
gases.  Thus,  the  uniform  system  is  to  be  recommended. 


C.  Conversion  of  (atm-cm)gjp  to  molecules 
cm~^[16] 

The  molecular  weight,  M,  of  any  gas  in  grams  occupies 
22.4136  liters  at  STP  [12]  so  one  cm3  of  any  gas  at  STP 
weighs 

M 

. j grams . 

2.24136x10^ 

Since  1 atm-cm<--rp  is  equivalent  to  one  cm  of  gas  at 
STP  per  cm2 

M -2 

(144)  1 atm-cm^-jp  = gm  cm  , 

2.24136x10^ 

Avogardro's  number  [35]  car  be  used  to  compite  the  number, 
n of  molecules  per  cm3  of  any  gas 

(145)  n = ■■  = 2.687016x10^°  molecules  cm'^. 

2.24136x10^ 

Hence  the  desired  result  is 


19  -2 

^146)  one  atm-cm^-j-p  = 2.6870x10  molecules  cm 
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D.  Conversion  of  Precipitable  Centimeters  of  Water 
to  mol  cm~^56] 


(147)  1 atm  -cmsTp  |^2q  = 


M 


2.24136x10 


18.016 

2.24136x10^ 


(148)  1 atm-cm^jp  |.j  q ~ 8.03798x10  ^ gm  cm"^. 

Using  the  constant  developed  in  paragraph  C»  above  we  have 
the  desired  result. 


(149) 


1 pr-cm  = 3.3429x10^^  mol  cm"^ 


E.  Partial  pressure  of  Water  Vapor  from 
Density  and  Tempera turel^~r. 


From  the  referenced  tables  we  have  the  relation 


(150) 

(5  • 

lO^a  4 p^/760 

= m m" 

where 

(151) 

a = 

weight  of  one  cubic  meter  of  dry  air  at  0°C 
and  760  torr  = 1.29278  kg. 

(152) 

6 = 

density  of  aqueous  vapor  relative  to  dry 
air  = .62168 

(153) 

Pa= 

partial  pressure  H^O  in  torr 

(154) 

a = 

coefficient  of  expansion  of  air  for  1°C  = 
,00366. 

Rearranging  we  liave  tlie  desired  result 


(155) 


p = .9456  (1  + .00366  T)  p torr 


where  T = temperature  in  deg  C 


= water  vapor  density  in  gm/m' 


p 


F,  Preci pi  table  centimeters  H2O  from  partial 
press u re,  temperature,  and  path  icr^qth 


(156) 


w = P' 


Using  tin;  relation  from  paragraph  C, 


(157/ 


1.05821  • 10‘“ 

1^“  00 36 7 {To|^~~T73. 1 3 ) 


where  ^ is  the  path  length  in  meters 

p^  is  water  vapor  pressure  in  torr. 

d 


Since  tlie  third  term  in  the  denominator  is  nearly  unity, 
we  can  write 


(158) 


2.8834x10"^  Pa 


y— pr-ciii 


6.  Conversion  to  atm-cm^jp  from  pressure, 
path  length,  and  telripeTature  of  sample  [16] 


atm- cm 


£Pa  273.13 


STP  760 


(159) 


atm-cm^.|.p  = 3.5938x10 


-1  ^ 


wliere 


absorber  pressure  in  torr 
path  length  in  cm 
temperature  in  deq  K. 


(160) 


H . ted  from  parti al  pressure. 

molecular  weight,  and  temperature 


This  is  of  course  the  Gas  Law 

p - 

( 

'T 


P = 1.6035x10  ^ ^ gin  cm  ^ 


where  p,  = pressure  in  torr 

a 

M = molecular  weight 
T = deg  K. 


1 1 Partial  pressure  of  water  vapor  from 
mixing  ratio  [14] 


(161) 


P,  = 


r P 


a .62168  + r 


torr 


where  r is  the  mixing  ratio  and  P is  the  total  pressure 
in  torr. 

From  the  data  in  the  AFCRL  models  [16],  one  should 
compute  mixing  ratio  as  follows: 


(162) 


Note  that  the  constant  .62168  in  Eq.  (7)  is  the  ratio 
of  the  molecular  weight  of  water  vapor  to  that  of  dry  air  at 
STP,  i.e.,  18. 0153/28. 97^5,  see  reference  35.  In  the 
Smithsonian  tables  [36]  che  value  used  is  0.62197  apparently 
based  on  older  molecular  weight  data. 


sure  from  density  for  AFCRL 


mul-lot.i  sumiiifr  model 


- X 760  X — 

0 1 loi  « 1 


“ 1.191  X lO'^ 

where  d = air  density  at  STP  = 28.9785(35), 

q = molecular  weight  of  ozone  = 47.9982  (35), 


-5  3 3 

and  6.0  x 10  g/m  is  the  density  of  ozone  and  1.191  x 10 

g/ttw  is  the  air  density,  both  at  sea  level,  of  the  AFCRL  mid 

latitude  surimer  model.  The  result  is  2.31  x 10-5  torr. 


APPENDIX  B 


ABSORPTION  CALCULATTjN  PROGRAM  LISTINGS 
A.  Single  Frequency  Calculation  Program 

A listing  of  the  program  is  given  on  pages  156  to  167.  The 
program  was  written  to  operate  on  the  tiectroScience  Laboratory 
time-sharing  system  and  contains  some  statements  which  are  incom- 
patible with  other  computer  systems. 

In  line  34  the  statement  OPTIONS  32K  allows  the  program  to 
use  over  60,000  words  of  program  storage  although  the  computer  only 
has  27,000  words  of  program  storage  available.  The  rest  of  the 
storage  is  on  disc  and  is  swapped  into  the  computer  memory  as  needed. 

In  line  49,  CALL  ESC  ($100)  causes  the  program  to  branch  to 
statement  100  immediately  if  the  escape  key  on  the  teletype  is  pressed. 

The  statement  CALL  COMDFL  ($100)  in  line  50  and  the  CALL 

C0MD( ) statements  in  lines  66,  73,  82,  96,  112,  119,  140,  162, 

209,  306,  364,  384,  389,  405  are  used  to  enable  the  user  to  branch 
easily  from  one  part  of  the  program  to  another.  When  the  'ALL 
C0MD( ) statement  is  executed  the  computer  types  < > on  the  tele- 

type and  waits  for  input  from  the  user.  The  user  responds  by  typing 
a character  string  of  1-6  characters  long  and  a carriage  return. 

The  computer  then  branches  to  the  first  executable  statement  follow- 
ing the  CALL  C0MD( ) statement  having  the  character  string  as  its 

argument.  For  example  if  the  statement  CALL  C0MD(5HPRESS)  in  line 
66  is  executed,  the  computer  will  type  < >.  If  the  user  types  PRESS, 
the  computer  will  branch  to  the  statement  in  line  70.  If  the  user 
types  SUB,  the  computer  will  branch  to  the  statement  in  line  86.  For 
another  computer  system  a similar  set  of  subroutines  could  be  written. 
They  would  have  to  be  written  in  assembly  language  and  would  require 
some  knowledge  about  how  arguments  are  passed  to  subroutines.  Alter- 
natively, a control  routine  could  be  written  in  Fortran  wliich  per- 
formed the  same  function  using  UF  statements  and  computed  GO  TO 
statements . 

The  CALL  ASSIGN  ( ) statements  in  lines  185,  218,  353,  and 

399  are  associated  with  the  ElectroScience  Laboratory  timesharing 
file  system  and  are  used  to  assign  logical  unit  numbers  to  physical 
devices  or  to  disc  files. 
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The  TPFILE  and  TPRKD  subroutines  called  ih  lines  244  and  245 
an:  system  subroutines  which  are  used  to  skip  files  and  records  on 
a magnetic  tape. 


For  the  single  frequency  calculation  program  there  is  a limit 
on  the  size  of  the  frequency  interval  over  which  calculations  may 
be  made  at  one  time.  This  limit  is  caused  by  the  size  of  the  array 
in  which  the  absorption  line  data  is  stored  (AD(4,5600)  in  line  4 ). 
For  a set  of  calculations  covering  the  frequency  range  VL  to  VH, 
the  program  requires  all  Ihe  data  for  absorption  lines  of  the  correct 
Substance  with  strength  greater  than  SLOW  in  the  frequency  range 
from  VL  - BOUND  to  VH  + BOUND.  For  the  ElectroScience  Laboratory 
computer  system  the  maximum  storage  allowed  for  the  program  and  the 
data  is  60,000  words.  TL’s  is  enough  room  for  the  calculation  pro- 
gram and  the  data  for  5600  absorption  lines. 


The  data  is  stored  in  one  4 by  5600  array  rather  than  four  560- 
word  arrays  in  order  to  increase  the  efficiency  of  the  program  on  a 
timesharing  computer. 


t 


1 

2 

3 

4 

5 

6 
y 
fi 
9 

10 

11 

12 

13 

14 

15 

16 
17 
16 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 
36 

39 

40 

41 

42 

43 

44 

45 

46 

47 
46 

49 

50 

51 

52 

53 


C 

C 

h# 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c***« 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


PROGRAH  NAME>  LEVEL 

COPYRIGHT  1975  THE  OHIO  STATE  UNIVERSITY 

CALCULATES  SINGLE  FREQUENCY»LEVEL  PATH 
ABSORPTION  coefficients, 

PROGRAM  authors*  C.  L-  TRUSTY  ANC  F,  3.  MILLS 

The  OHIO  state  university 
electroscience  laboratory 

1320  KINNEAR  road 
COLUMBUS,  OHIO  43212. 

USABLE  commands  ARE  PRESS,  BROAD,  NALFA,  ETa*  ASS, 
TEMP,  LINES,  EXc.  END.  DATA*  SUB,  STORE,  VAR, RESULT. 

CALCULATION  CAN  BE  HADE  FOR  A MAXIMUM  OF  175 
FRECUENCIES,  LAST  ENTRY  IN  DATA  FILE  SHOULD  BE  0.0 
IF  FILE  contains  FEWER  THAN  175  LINES.  FIRST 
LINE  MUST  BE  LOWEST  FREQUENCY  AND  LAST  LINE  HIGHEST. 

** 

LOGICAL  DEVICE  ASSIGNMENTS  FOR  THIS  PROGRAM 
0 AFCRL  LINE  DATA  TAPE  (INPUT) 

2 file  containing  CALCULATION  FREQ.  (INPUT) 

3 file  of  frequencies  AND  CALCULATED 
absorption  coefficients  (output) 

6 .OUT  file  (line  PRINTER  OUTPUT) 

7 MSG.3271A  (PROGRAM  WRITES  TIME  WHEN 
OFF-LINE  calculation  IS  COMPLETED) 

6 TELETYPE  (OPERATOR  COMMUNICATION, 

INPUT  OR  OUTPUT) 

JJ  OUTPUT  DEVICE  SELECTOR  6 OR  6) 

OPTIONS  32K 

LOGICAL  FIRST, RUOUND.  TELE 

Real  halfa 

DIMENSION  ICARD(5),ABS(5).lDATE(3) 

COMMON  JJ 

DIMENSION  Y(5,175).VL( 175).BX(7),CX(7).MA(7), ITIH(3) 
C8MM0N/Z/ AD  (4.5600,0  AY 

DATA  ITTY/3HTTY/,e0UHr/2O.0/,TEMPO/296.O/,S  .0W/.1E-26/ 
DATA  MA/16,44,49,44,?8.16,j7/,  IYF-S/ JhYES/ 

DATA  CX/. 62. .56. 0.5, 0.5, 0.5* 0.5, 0.5/ 

DATA  BX/1 .5,1 .0,1.5, I .0,1 .0*1. 5, I .0/ 

PRINT  LIST  OF  USABLE  C0MM4NDS  ON  TELETYPE. 

WRITE(8,2070) 

CALL  ESC(SIOC) 

CALL  COrlDFL(SlOO) 

FLAGS  USED  THRU  THE  PROGRAM* 

FIRST  — signifies  FIRST  TIME  THRU  THE  PROGRAM 
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54 

C 

RU6UND  - DATA  FILE  HAS  SEEN  REUBIIND 

55 

C 

TELE  FREOUENCIES  TB  BE  ENTERED 

ON 

teletype 

56 

C 

57 

JJ-6 

58 

FIRST-. TRUE. 

59 

RUBUND-. false. 

60 

TELE-. TRUE. 

61 

C 

62 

C 

WRITE  DATE  AND  TIME  BN  THE  TELETYPE 

• 

63 

C 

64 

CALL  CLOCK 

65 

GB  TB  150 

66 

100 

CALL  CBHDCSHPRESS) 

67 

C 

68 

c 

READ  THE  TBTAL  PRESSURE  IN  TBRR. 

69 

c 

70 

ISO 

WRITE(8, 20001 

71 

READ(6»-)P1 

72 

IF(FIRST)CB  T9  200 

73 

CALL  C0rtD(4HTEMP> 

74 

c 

75 

c 

READ  THE  TEMP  IN  DEG  F. 

76 

c 

77 

REWIND  0 

78 

RUBUND-. TRUE. 

79 

200 

WRITE* 8, 19801 

80 

REA! (8,-ItF 

81 

IF(FIRST)GB  To  300 

82 

CALL  C6rtP(3(!SU8) 

83 

c 

84 

c 

GET  SUBSTANCE  ID  NUMBER 

85 

c 

86 

300 

URITE(8.3020) 

87 

READ(8,>INSU8 

88 

C 

89 

C 

GET  ISBTBPE  INFORMATION 

90 

C 

91 

WRITE(6*3030) 

92 

read(8,-hsbtbp 

93 

IFIFIRSTIGO  '(&  400 

94 

RUBUND-. true. 

95 

REUIND  0 

96 

CALL  C3MD(3hA3S) 

97 

C 

98 

C 

READ  the  ASSBRBEK  PRESSURES  IN  TBRR, 

5 

MAXI PUN. 

99 

c 

100 

400 

URITE(8,2010) 

lOi 

C 

102 

C 

NTORRS  IS  THE  NUMBER  BF  ABS8RBER  PRESSURES  TO  BE  USED 

103 

C 

104 

READ(fl.-»NTBRRS 

105 

IF (NTeRRS.CT.5)NT0RR9-5 

106 

URITE‘8,2020) 

107 

c 

108 

c 

ABS(I)  are  the  ABSORBER  PRESSURES  IN 

tbrr 
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109 

110 
111 
112 

113 

114 

115 

116 
117 
lie 

119 

120 
121 
122 

123 

124 

125 

126 
127 
12S 

129 

130 

131 

132 

133 

134 

135 

136 

137 
13e 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 
163 


C 

READte.-l  (ABSM  1 , 1 -1  .NTORRS > 

IFIFIRSTIGO  To  500 
CALL  C0MD(5hPR0AD» 

C 

C READ  the  SELF-BROADENINC  COEFFICIENT. 

C 

500  URITEI6. 19901 
READie.-IQSOAli 
IFIFIRSTIGO  TO  tOO 
CALL  C0Hn(5HNALFA) 

C 

C read  the  NUHPER  of  H4LFUID1H3  WHERE  LINE  SHaRE 
C flODIFICATlOK  TAKES  PLACE. 

C 

600  JRITE(8,1970) 

READ(8,-)nALFA 
IFIFIRSTIGO  TO  700 
C 

C CALCULATE  CNORH,  THE  NORMALIZING  FACTOR  FOR  MODIFIED 

C line  shape,  see  mills  notebook  -4320  PC  104  FOR  DETAILS, 

c NOTE  IN  The  Following  lines  that  cnork  is  computed 

C at  two  places,  this  must  EE  DONE  IN  ORDER  FoR  THE 

C TIMESHARING  COMD  FUNCTION  TO  WORK  PROPERLY  FOR 

C BOTH  NALFA  AND  ETA. 

C 

IFIETA.E0.2.IG0  TO  650 

CN0RM-3.14159/I2.*( ATANINAlFAI+NALFA/I INALFA*NALFA+1.»* 

1 lETA-I .) ) ) ) 

GOTO  690 
650  CNORM-1. 

690  CALL  C0MDI3HETAI 
C 

C read  the  modifier  POWER  (ET*>. 

C 

700  WRITE(e.l960) 

READ(6,-IETA 

C 

C CHECK  TO  SEE  IF  ETA  IS  LESS  THAN  ONE. 

C 

IFlETA.GV.l.lGO  TO  710 
URITE(8,2o50) 

GO  TO  700 
C 

c calculate  cnorm,  the  normalizing  factor  For  modified 

C LINE  SHAPE.  SEE  MILLS  NOTEBOOK  -4320  PC  104  FOR  DETAILS. 

c 

7lO  IFIETA.E0.2.IGO  TO  730 

CN0RM*3.U15S/(2.*(ATANINALFA)+N ALFA/I  ( NALFA«NALF A4 I . I * 
IIETA-I.mi 
GOTO  740 
730  CNORM-1. 

740  IFIFIRSTIGO  TO  800 
CALL  COMDISHLINES) 

C 
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164 

C 

READ  CALCULATION  FREOUENCIES  INTO  VL,  175  MAXiNUH. 

165 

C 

166 

REWIND  0 

167 

RUeUND-.TRUE. 

166 

TELE  - .true- 

169 

close  2 

170 

C 

171 

C 

ARE  EREOUENCIES  TO  BE  ENTERED  6N  TELETYPE 

172 

C 

173 

600 

URITE(8. 10201 

174 

READ(6.1030)  I IN 

175 

IFIIIN.NE.ITTYITELE-.FALSE, 

176 

IFITELEICO  T6  155 

177 

C 

178 

C 

REQUEST  location  OF  FILE  WiTH  CALCULATION  FREQUENCIES 

179 

C 

LAST  LINE  OF  FILE  MUST  BE  0-0  IF  THERE  ARE  FEWER  THAN 

160 

C 

175  calculation  FREQUENCIES- 

161 

C 

162 

URITE(6»1C00) 

183 

READ(8, lOlOIFiLEL 

184 

READ(8-1010)USERL 

165 

CALL  ASSIGN(FILEL»USERL»2) 

186 

LUNlT-2 

167 

GO  TO  160 

188 

155 

LUNlT-8 

189 

160 

CONTINUE 

190 

IFILUMIT.EQ.? :UR ITE(6»2030> 

191 

C 

192 

193 

C 

c 

IVL  IS  THE  NUMBER  OF  CALCULATION  FREQUENCIES  USED. 

194 

IF(LUNIT.EQ.e)READ(6,-HVL 

195 

IFdUMT.EO  .6)  write  (8,20401 

196 

IF(LL'MT.E0.2)  IVL«175 

197 

DO  110  1-1, IVL 

198 

IF(LUN1T.EQ.P)WR1TE(8,1950»I 

199 

c 

20C 

c 

CALCULATION  FREQUENCY  FILE  MAY  HAVE  BAND  ID»S  FO®  USE 

201 

c 

REPORT, 3271T 

202 

c 

203 

READ(LUNIT-3000JVL( 1 1 

204 

3000 

F0RMAT(FI3.0) 

205 

IF( VL( n .EQ.O.OICO  to  130 

206 

110 

CONTINUE 

207 

130 

IVL-I-1 

206 

if(fipst)go  To  950 

209 

CALL  C0M0(4hDaTAJ 

210 

ruound-.true. 

211 

C 

212 

C 

REQUEST  LOCATION  OF  LIME  DATA 

213 

C 

214 

900 

URITE(8,950) 

215 

WPITE(8,I001) 

216 

READ(8,101C)FILED 

217 

CLOSE  0 

216 

CALL  *SSIGN(FILED,0.0,0I 
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219 

C 

220 

C 

GET  THE  FILE  AND  REC9PD  NUMBERS  eN  THE  TARE 

1 221 

C 

WHERE  THE  PERTINENT  ABS8RPtI0N  DATA  |S  LOCATED. 

• 222 

C 

223 

WRITEie. 30101 

224 

READia.-lNFlLE.NRKD 

225 

nfile«nfile-i 

226 

NRKD-NRKD-1 

227 

840 

IF(,N6T.FIRST»G6  T6  115 

226 

C 

229 

C 

SHOULD  calculation  BE  PERFORMED  OFF-LINE  ? 

230 

C 

231 

URITE(8. 30401 

232 

REAO(6.I030IIOUH 

233 

IFdOUM.EG.IVESICALL  DEASSN 

234 

FIRST*. false. 

235 

r 

236 

c 

J IS  THE  total  number  OF  SUBSTANCES  TO  BE  CONSIDERED. 

237 

c 

SET  NUMBER  equal  1 FOR  NOW 

236 

c 

239 

850 

J»1 

240 

C 

241 

C 

POSITION  data  tape. 

242 

C 

243 

REWIND  O 

244 

CALL  TPFlLE(O.NFlLE) 

245 

CALL  TPRKn(O.NRKD» 

246 

I ■ 0 

247 

C 

2(6 

c 

READ  DATA  FROM  AFCRL  LINE  DATA  TAPE, 

249 

c 

FOR  LINE  I.  AD(1..')  IS  THE  FREQUENCY,  AD(2,I>  IS 

250 

c 

the  LINE  strength,  Ari(3,I)  IS  THE  HALFUIDTH,  AND 

! 251 

c 

AD(4,II  IS  the  energy  OF  THt  LOWER  STATE  OF  THE 

2 

5 

152 

c 

transition. 

' 2 

154 

GO  TO  22 

i 2 

155 

20 

CLOSE  0 

f 2 

5e 

57 

GO  TO  23 

2 

I 2 

56 

So 

c 

r 

CHECK  FOR  SLOW 

2 

60 

21 

IF<  (AD(.',I  I.LT.SLOU)  ,0R.(  ADU,I  ) .LT  ,VL  ( I) -BOUND  ) ) l-l-l 

2 

61 

IF(AD<l,n.GT,VHlVL)AB0UND>GO  TO  27 

2 

62 

22 

I»I  + 1 

2 

63 

IF( I.GT.5600IC0  TO  27 

2 

2 

64 

65 

23 

r 

RE AD <0,3. END-20  I AD ( I , I ) , AD < 2, I) . AD ( 3. I ) . AD ( 4 , I ) , i STOP. N ABS 

2 

2 

66 

67 

•h. 

r 

CHECK  FOR  PROPER  SUBSTANCE  AND  ISOTOPE, 

2 

66 

IF(NAES,NE.NSU8»G0  to  23 

2 

69 

IF( ISCTOP.EO,0»GO  TO  26 

2 

70 

IF< ISOTOP, NE.ISTePICe  TO  23 

2 

71 

26 

IF(ADn,m?1.27,21 

2 

72 

27 

ICARDI JI-I-l 

2: 

73 

CLOSE  0 

1:74 

P-Pl 

275 

ICRD-ICARDU) 

276 

C 

277 

C 

TEMP  IN  OeG  K is 

278 

C 

279 

TEMP-5. /9.*(TF-32.) ♦273. 2 

280 

icrd»icard< J> 

281 

II-ICARD( J) 

282 

17-1 

283 

c 

284 

c 

CALCULATE  TEMPERATURC  C9RRECTieN  C9N3TANTS. 

285 

c 

286 

52 

CS1-(TEMP9-TEMP)/<TEMP9*TEMP*.6951 ) 

287 

53 

CS2-(TEMP9/TEMP)**BX(NSUB) 

288 

54 

CA-(TEMP9/TEMP)**CX(NSUB) 

289 

C 

290 

C 

CALCULATE  THE  C9PPLER  HALF-UIDTH  DIVIDED  BY  FREQUENCY 

29l 

C 

292 

DN2-AL9G(2.) 

293 

ALFaD-3.5812E-7*S0RT(TEMP/MA(NSUB) ) 

294 

C 

295 

C 

CHANCE  THE  STRENGTHS  AND  WIDTHS  T8  NEW  TEMPERATURE 

296 

C 

297 

D9  71  I-I.Il 

298 

AD ( 2. I ) -AD <2. I )-CS2»EXP(-Ad‘<»I >*C91 ) 

299 

C 

300 

C 

PRESSURE  C9RRECTI9N  T9  HALFRIDTH  (Ca)  IS  MADE 

301 

C 

IN  ABSV92  9UBR9UTINE. 

302 

C 

303 

AD(3. I )-A0(3> I> *CA 

304 

71 

C9NTINUE 

305 

G9  T9  116 

306 

115 

CALL  C9MD(3hEXE> 

307 

C 

308 

c 

N9TE«  IF  EXE  C9MMAN0  IS  GIVEN  AFTER  C9MPLETI9N 

309 

c 

9F  A CALCUI.ATI9N  WITH  N9  VARIABLE  CHANGES,  THE 

310 

c 

CALCULATION  WILL  BE  REPEATED  WITH  SAME  DATA  AND 

311 

c 

IF  9N-LINC  FNTER  CUTPUT  DEVICE  NUHBER  UILL  BE 

312 

c 

PRINTED  almost  IMMEDIATELY  LEADING  USER  T9  THINK 

313 

c 

that  Nf  NEW  CALCULATION  HAS  TAKEN  PLACE.  IF  NEW 

314 

c 

DATA  IS  WANTED  DUE  E.G.  T9  SUSPECTED  TAPE  READ 

315 

c 

PR9ELEM,  give  DATA  COMMAND  FIRST  WHICH  .ILL  CAUSE 

316 

c 

TAPE  T9  BE  REW9UNC  A‘iD  READ  AGAIN. 

317 

c 

318 

WRITE(8, 30401 

319 

READia.injOI lOUM 

320 

IF( IDUM.EC.IYeSICALL  DEASSN 

321 

IF(RW9UNDIG9  T9  850 

322 

IF( IDUM.EC. IYESIG9  T9  116 

323 

116 

C9NTINUE 

324 

D9  90  1-1,-IVL 

325 

D9  85  M-1.NT9RRS 

326 

c 

327 

c 

F9R  A 9NE  km  path  THE  ABSORBER  AMOUNT 

328 

c 

IN  M0LECULES7CM2  IS 

161 
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329 

C 

330 

U ■ ,733952E22*( ABS(H)/760,/TeMP)*1E5 

331 

C 

332 

C 

EFFECTIVE  PRESSURE  IN  ATM  iS, 

333 

C 

334 

P«(P1*(BR8AD-1.)«A8S(M) )/760. 

335 

V«VL(  I ) 

336 

CALL  A8SVe2<V,B8UND.l7. 1 1 .f » N ALFA. ETA* ALFAD ) 

337 

CAYL»CAY*U«CNORM 

338 

C 

339 

C 

CAYU  D8ES  N8T  HAVE  UNITS.  ]T  EQUALS  LH 

340 

C 

8F  TRANSMITTANCE. 

341 

C 

Y(H.»)  WILL  be  used  FftR  THE  ABS8RPTI8N  C8EFF 

342 

C 

IN  KM-1,  I.E.  Ln(T1/1.0 

343 

C 

t 

344 

Y(M.n-  + CAYu/1.000 

345 

85 

C8NTINUE 

346 

90 

C8NTINUE 

347 

C 

348 

C 

JRITE  MESSAGE  IN  MSC.3271A  TELLING  WHEN  THE  CALCULAT: 

349 

C 

IS  finished. 

350 

C 

351 

CL8SE  6 

352 

CL8SE  7 

353 

CALL  ASSICN(4hHSG  .5H3271A.7) 

354 

CALL  GETDaTI IDATE) 

355 

CALL  GETTlM(lTIM) 

356 

URITE(7,3050» ITIM. lUATE 

357 

CL8SE  7 

358 

C 

359 

C 

WHEN  reassigned  AFTER  8FF-LINE  CALCULATI6N  PPBGRAM 

360 

C 

WRITES  executive  SIGN  (<>>,  RESULTS  MAY 

361 

c 

THEN  BE  DBTaINEO  BY  USE  8F  STORE  COMMAND 

362 

c 

(abbreviated  VERSION!  8R  RESULT  COMMAND. 

363 

c 

364 

CALL  C8MD(6HFESULT,0) 

365 

c 

366 

c 

write  out  final  answers  if  calculation 

367 

c 

IS  ON-LINE  . 

368 

c 

369 

c 

3»0 

c 

SHOULD  calculation  RESULTS  BE  PRINTED 

371 

c 

ON  TELETYPE  OS  WRITTEN  ON  .OUT  FILE  ? 

372 

c 

373 

WRITE(8,1940) 

374 

8EAD(8,-» JJ 

375 

CALL  CLOCK 

376 

WRITE!  JJ.  i3  )B)((  NSUBI.CXI  NSUO)  .BOUND.SLOW.  ICARD(  J) 

377 

WRITE! JJ/206C ) CNOPM 

378 

DO  666  I«l,lVL 

379 

WRITE! JJ.120»VL(  I !,(  ABSIK  ) .Y(k. I ) . K» 1 j NTORRS > 

380 

120 

FORMAT(/»  AT  THE  FREDUENCY' F 1 0 .3// 

381 

1'  ABS  'eX*K»6X/'  T0RR'6X» I /KM» 

382 

221(/1X,0Rf6.2,1PEI1.3) ! 

( 
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i 

i 

4 


I 


383 

666 

384 

385 

C 

386 

C 

387 

C 

388 

389 

390 

C 

391 

C 

392 

C 

393 

C 

394 

395 

396 

397 

396 

399 

400 

401 

402 

10 

403 

404 

40S 

406 

407 

1 

4 08 

3 

409 

13 

410 

411 

412 

950 

413 

1000 

414 

lOCl 

415 

lOlO 

416 

1020 

417 

418 

1030 

419 

19''0 

420 

1950 

421 

I960 

422 

1970 

423 

1980 

424 

1990 

425 

20C0 

426 

2010 

427 

2020 

428 

2030 

429 

2040 

430 

2050 

43l 

2060 

432 

2070 

433 

434 

435 

2080 

436 

3010 

437 

3020 

CONTINUE 

CALL  COHDOhVaR) 

PRINT  LIST  OF  PReCRAM  PARAMETERS  ON  THE  TELETYPE 


URlTE(8.3060)Pl,ABsm,TF.NALFA,ETA,BROAD,NsUB,  ISOTOP 
CALL  C0MD(5HST0RE» 

NOU  WRITE  partial  DATA,  TF  DESIRED,  FOR  USE  IN  FURTHER 

data  reduction  procrams. 

WRITE(fi,4000) 

WRlTE(fl,1000» 

READ(8,1010)FILES 
READ(8,1010)USERS 
CLOSE  3 

CALL  ASS1GN(F|LES, USERS, 3» 

DO  10  I-1,1VL 
UHlTE(3.2080»VL(n,Y(l,I» 

CONTINUE 
CLOSE  3 

RUOUND-. False, 

CALL  C0MD(3HEND» 

CALL  EXIT 
F0RMAT(6E10.4J 

F0RmAT(2F1O.4,f5.2,F1O.2,38X,14, i3) 

FORMAT!'  BX«*F5.2/*  CX  -»F5*2,/ 

2'  BOUND  *'F7.2/,'  SLOW 

3£i,2,/*  NUMBER  OF  LINES  USED  •»,15/) 

FORMAT!'  where  IS  THE  LINE  DATA?' » 

FORMAT!'  ENTER  FILE  AND  USER  NAME  ON  TWO  LINES.*/) 
F0RMAT!'ENTER  .MTO  oh  .MTl  */) 

F0RMAT!A6) 

FORMAT!'  WILL  CALCULATION  FREQUENCIES  BE  READ  FROM* 

2*  A FILE  OR  THE  TTY  7*) 

FORMAT! A3) 

FORMAT!*  enter  OUTPUT  DEVICE  NUMBER,*) 

FORMAT! 13, 2H  -) 

FORMAT!'  enter  MODIFIER  POWER.*) 

FORMAT!*  enter  MODIFICATION  LOCATION  IN  HALfU I DTHS.* ) 
FORMAT!*  enter  THE  TFMP  IN  DEC  F.*) 

FORMAT!*  enter  THE  SEI.F-BRO  .DEN I NC  COEFFICIENT.*) 

FORMAT!*  enter  THE  TOTAL  PRESSURE  IN  T6Rk.*) 

FORMAT!'  >,0U  MANY  ABSORBER  PRESSURES  WILL  BE  USED?') 
FORMAT!'  enter  THEM  ON  ONE  LINE,  FREE  FORMAT.'/) 

FORMAT!'  HOW  MANY  LASER  LI.NES  WILL  BE  USED?') 

FORMAT!'  enter  THEM  ONE  PER  LINE  IN  INCREASING  ORDER.*/) 
FORMAT!*MfDIFlER  POWER  MUST  BE  STRICTLY  GREATER  THAN  11*) 
FORMAT!*  THE  NORMAL  1 7ATI0N  FACTOR  ■ ' ,F5.3> 

FORMAT!'  commands  ARE  '/ 

2'  PRESS,  ABS,  temp,  NALFA,  ETa,  LINES,  BROAD,' 

3'  DATA,  EXE,  END,  SUB,  STORE,  V AP, RE'ULT. *// ) 

FORMAT! 1X,F1C .3,1PE12.3) 

F0RMAT!'ENTER  tape  start  pointifile  number  and  RR-CORD  no.*) 
FORMAT!*EnTER  id  no.  of  absorber*/ 
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43S 

439 

3030 

440 

441 

3040 

442 

3050 

443 

3060 

444 

40C0 

445 

446 

c 

447 

c 

448 

c 

449 

c 

450 

451 

452 

453 

c 

454 

c 

455 

c 

456 

c 

457 

c 

458 

459 

460 

461 

462 

12 

463 

464 

14 

465 

15 

466 

467 

19 

468 

17 

469 

25 

470 

471 

472 

473 

474 

C 

475 

<• 

476 

C 

477 

c 

478 

479 

480 

c 

481 

c 

482 

c 

483 

c 

484 

485 

c 

486 

c 

487 

c 

488 

489 

490 

c 

491 

c 

492 

c 

2 » 1-H28  2*C02  3-6.1  4«N2p  5“C0  6-CH4  7«82»/» 

F8RMAT(*  ENTER  NUHRER  BE  PfSlRED  ISOTBPE.*/ 

1 * IF  ALL  ISPTBPES  arE  DESIRED  ENTE,'  0-*) 

FBRMATI'DB  YBU  UISH  T»  RUN  BEE-LINE?*) 

F8RMAT (*CALCULAT!BN  FINISHED  AT  *3A3*.  *3a3) 

F8RNAT(/F6.1»1PE11.2.0PF7.1»F7.1.F7.2»F6.1.i4,I«/) 

FBRNATI'UhErE  mill  data  be  STeREC?') 

END 

THE  F8LL8UING  SUBRBUTINE  CALCULATES  THE 
A0S8RPTIBN  C8EFEICIENT. 

SUBRBUTINE  APSVB2 I V.BBUN 0* I 7, I 1 ,P .N ALFA»ET  A, ALE AD ) 

REAL  NALFA 

CBMH8H/Z/AD(4,5600).CAY 

THE  TUP  EBLLBWING  LBBPS  DETERMINE  WHICH  ABSpRPTlBN 

LINES  WILL  BE  USED  TP  CALCULATE  THE  aBSBHPTIBN  CBEFFICIENT 

AT  THIS  FREQUENCY. 

15-1 

le-ii 

DB  14  I-l.Il 

IF(V-BeUND-AD(l»I))12.12.14 

I5-I 

G8  T8  15 
C8NTIHUE 
D8  19  K-I7.il 

IF( V*B8UNn-AC(l,K) ) 17.19.19 

C8NT INUE 

16  • K - 1 

CAYl-C.O 

CAY2-0.0 

CAY3-0.0 

DN2-AL0C(2.) 

DM-SGRT<DN2) 

THE  FPLLPUING  LBBP  SUMS  THE  CBKTR I BUT  1 8NS  FrBH 
ALL  THE  APSaRPTlBN  LINES. 

D8  46  I-I5.IB 
Y-A0S( V-ADl  1.  I ) ) 

THE  FPLLBUING  INSTRUCT'BH  CORRECTS  THE  LBREnTZ 
HALFUIDTH  FBR  PRESSURE. 

PT-AD(3.I)*P 

THE  F8LLPUINC  STATEMENT  CALCULATES  THE  DBPPLER  WIDTH, 

ALFD-AD( l.I )*ALFAD 
XV-DNl»Y/ALFD 

IF  absbrptibn  unE  is  far  away  FPBH  CALCULATIBN 
FREQUENCY,  use  LBRENTZ  SHAPE  REGARDLESS  BF  PRESSURE 
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493 

494 

495 

496 

497 
4 9a 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 

513 

514 

515 

516 

517 
516 

519 

520 

521 

522 

523 

524 

525 

526 

527 

528 

529 

530 

531 

532 

533 

534 

535 

536 

537 

538 

539 

540 

541 

542 

543 

544 

545 

546 

547 


C 

C 

c 

c 


c 

r 

c 


50 


C 

c 

36 


C 

C 

C 

42 


46 

C 

C 

c 

c 

c 

c 


c 

c 

c 

c 


10 


c 

c 

c 

c 

c 


3UT  IF  L8BENT2  HALFUinTH  IS  LESS  THAN  FIVE  TIHES 
THE  DOPPLEB  HALFUI’JTH.  USE  VBIGT  LINeSHAPE. 
OTHERUISE  USE  MOD'FIED  LBRENTZ  SHAPE. 


IFIXV.GE. 300.1  Ce  TO  53 
IF (PT/ ALFd.CE.S. IGO  T9  50 
YV-PT/ALFD*DN1 


CALCULATE  THE  V9IGT  SHAPE 

CAY3«CAY3  + AD(2.n*DM/(SaRT(3.l4l592i*ALFD)»V0IGT(XV,  VV) 

UO  TO  46 

auay-nalfa»pt 

IF(Y-AUAY136.36.42 

CALCULATE  THE  LORENTH  LINE  SHAPE  FOR  V HEAR  LINE  CENTER 

SUMl»An(2, I MPT/(Y**2»tPTl**2l 

CAY1«CAY1*SUH1 

GO  TO  46 


CALCULATE  MODIFIED  LINE  SHAPE  fop  V ) HALF A *A D 1 3. I) *P 

AUAY30-AUAY*AUAY 

CK2-ADI2. I )»PT/IAUAYSQ4FT*pT) 

SUH2>CK2*(AUAY/Y1**ETA 

CAY2-CAY2+3UM2 

CONTINUE 

summation  oviR  all  lines 

CAY  ■ 0.3163*(CAY1  ♦ CAY2+CAY3) 

absorption  COeFF,»  cay.  HAS  UNITS  OF  (M0L*CN-21-1 

RETURN 

END 


S'JBROUT/  'E  CLOCK 
COMMON  JJ 

DIMENSION  ITIMEI31.IDATE(3) 

CALL  GETDATdOATEl 

CALL  GETTIMI ITIMEI 

URITEUJ.IOI  I DATE.  I TIME 

FORMATC  date  »3A3»  TIME  '3A3/1 

RETURN 

END 

THE  FOLLOWING  SUBROUTINE  EVALUATES  THE  VOlCT 

profile  expression,  subroutine  WAS  provided 
BY  CHARLES  YOUNG. 


I 

i 


I 

\ 

i 

I 
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n n tm;S4-n  P,-. 131664, 0.0..«77081^E-1.0.0.-.S 
3 a?i412F*»  *0~ol*^36°i57iE-l»0.oi -.2  J64976E- 1,0.0,. ) J 1960  IE- 1*0.0,  •• 

3 JU5!loU3:2:r.!iJ^X«E-3.0.r.-.1173«rt-2.0.0.WJ.«20£^ 

. •«>  A 700B  77kE-4  ■ 0 «0  # • . 2f  1 #0  #0»  • E ..#U* 

5 •;:i!”S?««;5;S:j”«6«7.E;L-0.o;-.33...3.E-<.0.0..70,3303E-7, 
X-XIN 
Y-YIN 
X2-X*X 
Y2-Y«Y 

IF  (X-7.)  200,201,201 
200  IF  (Y-1.)  202,202,203 
203  HA<l)-0. 

CA<  l)-0. 

RBI  1 )«1  . 

CB( 1 )-0. 

RAI21-X 

CAI21-Y 

RBI2)-.S-X2+Y2 

CB(2)--2.*X«Y 

CB1-CEI2) 

UVl-0. 

DO  250  J"2,3l 

JMlNUS-J-1 

JPLUS-J+1 

FL8ATJ-JM1NUS 

RB1»2.«FU0ATJ*’RB(2) 

RAl»-FLeATJ*I2.*FL6ATJ-l.>/2. 

RA( Jvl 1cPu1*PAIJ)-C91*CA(J)*R*1*FA< J-*  > , 

CA( J*1 l-PHi«CA( J)*CP1*RA( J)*BA1*CAI JNINUS) 

RBI J+1I-«HI*FB(J)-CB!‘CB(J)^RA 



1 CBIJPLUi)«CBIJPLU3H 
IF  I ABS I UV-UVl 1-1 ,E -6 1251 ,250,250 

250  UVl-UV 

251  V81GT«UV/1 .772454 
RETURN 

202  IFIX-2. >301,301,302 

301  AINT  ■ 1. 

RAX"12.-*5.»X2 

kmax-rax-i 

K8-0 

DO  303  K*K8,Kf1AX 
AJ-MAX-K 

303  aINT*AINT«I-2,*X2)/I2.*AJ+1.1+1. 
g,-2.«X*AlNT 
C8  T8  304 

302  IF  IX-4.5)  305,306,306 
305  3I431'*0. 

R(44la0. 


J-42 

DO  307K-1.42 

9(J»«.4.X.8(J  + n-0(j*2)*A(j) 

307  J-J-1 

U»B(3)-B( 1 » 

SO  TO  304 
306  AINT-1.0 

HAX-2.+40./X 

AMAX-f^AX 

DO  306  K-l.HAX 

AINT«AINT.(2.*AMAX-1,»/(2.*X2»*1. 

306  AHAX-AMAX-1. 

U — AINT/X 

304  V»1.772454*EXP(-X2J 
H-,02 
JM-Y/H 

1F(JH)310.311.310 

311  H-Y 
310  Z-0. 

L-0 

DY( 1 )»0. 

312  DY(2)»H/2. 

DY(3»-DV(2» 

DY(4)»H 

318  AK ( 1 )«0. 

AM( 1 )»0. 

DO  313  J-1.4 
YY-Z+DY( JJ 
UU«U+.S»AK( J» 

VV«V*.5*AM( J» 

AK(J  + n»  2.*l  YY*UU*X*VV»*M 
AH( J+1  )«-2.*!l,*X*UU-YY*VV»*M 
IF(J-3I  313,314,313 

314  AK(4  )«2.*AK(4» 

AH(4)>A)1(4)aAM<4» 

313  CONTINUE 
Z"i.  'H 
L-L  + 1 

U-U+.1666667*(AK(2»*2.*AK<3>»AK(4»*AK(5) ) 

V«V*.16e6667«(  AH(2UAM(3>»AM(3)+AM(4»+A'115)» 

IFUM)  31F, 32(1, 315 

315  IF(L-JMI  v)18,3l7,320 
317  AJM-JM 

H-Y-AJM*H 
GO  TO  312 

320  V0IGT«V/1. 772454 
RETURN 
201  Fl-0. 

DO  330  J-1,2 

330  F1-F1+HH(  Jl/(Y2-MX-XX(  JM*lX-XX(J)  n+HH(J(/(  Y2+(X*XX(  Jl  )*(X  + XX(  J)  » 
1) 

V0ICT"Y*Fl/3. 1415927 

RETURN 

END 


B . ^prt^rg^^alculation  Program 

A listing  of  the  program  is  given  on  pages  169  to  182.  This 
program  was  also  written  to  operate  on  the  Ohio  State  University 
EloctroScience  Laboratory  timesharing  system  a"d  contains  some  of 
the  special  statements  described  in  Part  A for  the  single  frequency 
calculation  program.  Not  used  are  the  ESC,  COMDFL  and  COMD  sub- 
routines. 

The  data  storage  requirements  are  somewhat  relaxed  with  this 
program.  At  any  given  calculation  frequency  V,  the  computer  re- 
quires data  for  the  absorption  lines  having  a frequency  between 
V-BOUND  and  V+BOUND.  There  is  storage  in  the  program  for  data  from 
2500  absorption  lines.  Thus  if  BOUND  is  set  at  20  cm  , then  there 
may  be  no  more  than  2500  absorption  lines  in  any  40  cm-l  interval. 

The  absorption  line  density  per  40  cm"l  interval  on  the  AFCRL  tape 
is  less  than  2500  except  near  9.6  microns  where  tliere  are  many  ozone 
lines.  For  calculating  spectra  in  that  region  either  BOUND  must  be 
reduced  or  the  calculation  must  be  performed  on  a larger  computer. 

Some  modifications  have  been  made  to  the  subroutine  ABSVU3 
which  calculates  the  absorption  coefficient  at  each  frequency,  so 
that  it  operates  faster  than  the  subroutine  used  in  the  single  fre- 
quency calculation  program  or  the  subroutine  ABSCOE  written  by 
Deutschmann  and  Calfee  [16].  Each  time  the  subroutine  is  called, 
it  must  first  search  the  data  storage  array  AD(J,K)  to  determine 
what  range  of  K will  cover  the  absorption  lines  with  frequency  be- 
tween V-BOUND  and  V+BOUND.  In  ABSCOE  (Deutschmann  and  Calfee)  and 
in  the  subroutine  in  the  single  frequency  calculation  program,  this 
search  always  starts  at  K=1 . ABSV03  was  changed  so  the  search  starts 
^he  values  of  K determined  the  last  time  the  subroutine  was  called. 
This  saves  time  since  the  frequency  is  increased  by  only  a small 
amount  each  time  the  subroutine  is  called. 

Note  in  proof:  It  has  been  discovered  that  for  low-pressure, 

this  program  can  produce  "glitches"  in  the  final  plots.  The  dif- 
ficulty is  associated  with  the  transition  from  Lorentz  to  Voigt  line 
shapes  within  an  individual  plot.  This  in  turn  is  controlled  by  the 
program  logic  in  statements  533  to  542.  Future  versions  of  the 

program  will  include  a fix  for  this  problem.  Please  contact 
R.  K.  Long,  Ohio  State  University  ElectroSci ence  Laboratory  for 
further  information. 
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1 c 

2 C 

3 C 

4 C 

5 C 

6 C 

7 C 

8 C 

9 C 

10  C 

11  C 

12  C 

13  C 

14  C 

15  C 

16  C 

17  C 

16  C 

19  C 

20  C 

21 
22 

23 

24 

25 
Vi 
27 

26 
29 

•■’c 

31 

32 

33 

34 

3f  C 

36  C 

■37  C 

3«  C 

39  C 

40  C 

41 

42 

43 

44 

45 

46 

47  C 

46  C 

49  C 

50 

51 

52 

53 

54 


p:<6gb*m  name*  transplot 

COPYRIGHT  1975  THE  OHIO  STATE  UNIVERSITY 

THIS  PROGRAM  CALCULATES  THE  LOG  OF  THE  ABSORPTION 
COEFFICIENT  VERSUS  FREQUENCY  FROM  THE  AFCRL  LINE  DATA 
TAPE.  output  is  IN  A SERIES  OF  PLOTS  TFN  INCHES  UIDE. 

PROGRAM  AUTHORS*  G.  L.  TRUSTY  AND  F.  S.  MILLS 

the  OHIO  state  university 
electroscience  laboratory 

1320  KINNEAR  ROAC 
COLUMBUS.  OHIO  43212 


LASER  LINE  FREQUENCIFS  MUST  BE  ARRANGED  IN 
STRICTLY  ascending  ORDER.  l35  MAXIMUM,  IF' 

feufr  Than  125.  last  line  must  be  0.0 

OPTIONS  32K 
LOGICAL  ENDATa. first 

DIMENSION  IDATE(3).BX(7».CX<7».TITLEA(7).TITLEP(7) 
dimension  YllOOl ».VLIME(125).PP(7».B{7>.M(7).itIM{3) 
COMMON/2/AD(3.250C  ) . NABS  { 2SOO  ' . I(-,  1 5,  J6 
C0MM0N/A/PE(7)  ,U{7  ),pL.A0.a1.A2.VL.Vh.CAYi,',Av  FAD(7) 
COMMON/D/ENDATA.FIRST.NSUB{7I,numBR.CSI.CS2(7  ».CA(7).IS0T0P<7> 

DATA  IY2/JH  /.80UNn/2O.O/.TEMP0/296.0/.SL0U/ .1F.-26/ 
data  IYES/3HYES/ 

DATA  CX/.62.  .56.0. 3. P.5.0. 5»C. 5,0.5/ 

DATA  BX/1 .5,1. .0.1 .5. I. 0.1. 0.1. 5, 1.0/ 

DATA  M/18,44. 48. 44,28, 16.32/ 

CALL  F5RRI0) 

CALL  CLOCK 

NUMBER  IS  THE  TOTAL  NUMBER  OF  SUBSTANCES  TO  Be  C8NSIDERED.  II.F.  FOR 
C02  AND  H20  NUMBER  EQUALS  ?) 

NUMPLT  IS  The  number  OF  PLOTS  THAT  WILL  BE  MADE  UITH  EACH  PLOT 
CONTAINING  AN  INCREMENT  OF  1000.'»DELVl  WAVENUM3EPS 

URITE(8. 20001 
READ(8.-)NUMBR.NUMPLT 
WRJTEie, 30001 
READie,-)  VI.  DELVi 
V2«V14NUMPLT*LELVI 
DELVl-DELVl/lOOO. 

READ  THE  TOTAL  PRESSURE  IN  TORR 

URITE(S,40Q0) 

READie. -IpRESS 

PPLOT-PRESS 

PL-IOOO. 

WR1TE(8,6000) 
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A.' :: 


I 


I 


55 

56 

57 
56 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 
76 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 

101 

102 

103 

104 

105 

106 
107 
106 


C 

C READ  temperature  IN  DEC  f 
C 

READ(8,-)T 
983  TPL8T«T 
C 

C CONVERT  TEMPERATURE  T8  DEC  K 

C 

TEMP-273. 2*5. /9.*(T-32.) 

C 

C GET  ID  number,  isotope.  PARTIAL  PRESSURE.  AND  BROADENING 

C COEFFICIENT  FOR  EACH  ABSORBER. 

C 

URITEI8. 41001 
DO  400  I-t.NUMBR 
URITE(8.4900II 
READ(8.-|NSUB( 1 1 
J-NSUBI I ) 

URITEI8. 49501 
READ(8.-)I80T0P(JI 
URITEI8. 50001 
C 

C READ  ABSORBER  PRESSURE  IN  TORR 

C 

READ(8i.-)PP(  J» 

URITEI8. 55001 
C 

C READ  SELF  BROADENING  COEFFICIENT 

C 

READ(8.-)B( J| 

400  CONTINUE 
C 

C REQUEST  continuum  INFORMATION. 

C 

URITE(8. 61001 
READ(8.8500)  1 IN 
IFdlN.NE.IYESI  GO  TS  MO 
URITE(8.620C) 

READ(8.-)a0 

UR1TE(8.6300) 

READ<e.-)Al 
WK1TE(8.64C0I 
READ<8.-)A2 
URITEI8. 65001 
READ<8.-»VL.VH 
GO  TO  150 
140  AO-0. 

Al-0. 

A2-0. 

C 

C GET  TITLE  INFORMATION  FOR  pLOTS 

C 

ISO  URITEI8. 42001 

READI8.1011)TITLEA 


•? 


1 

! 


ifriTt-'  'liiiiiMiiiKTtf'i 


READ<8*10ll)TrTLCP 


109 

no 

111 

112 
113 
1.  i 

115 

116 
117 
lie 

119 

120 
121 
122 

123 

124 

125 

126 
127 
126 

129 

130 

131 

132 

133 

134 

135 

136 

137 
13a 

139 

140 

141 

142 

143 

144 

145 

146 

147 
146 

149 

150 

151 

152 

153 

154 

155 

156 

157 
15B 

159 

160 
161 
162 
163 


C 

C REQUEST  LBCATIBN  BP  ABSBRPTION  LINE  DATA 
C 

WRITE  18,950) 

URlTEie.llOO) 

READie.lOlOlFiLED 
CALL  ASSIGNIFILED. 0.0.0) 

UR1TE1A.9400) 

READ16  -InFILE.NRKD 

NFILE-NFILE-I 

NRKD-NRKD-1 

c 

C REQUEST  Tape  unit  NUH8EP  FBR  writing  PLBT  'TaPE. 
c 

WRITEia.1200) 

WRITE18.1100) 

READie.lOlOlFiLEP 
CALL  AS3IGN(FILEP.0.0,2) 

WRITE (8. 1300) 

READ(8.-)NPFILE 

NPFILE-NPFILE-1 

C 

C REQUEST  LBCATIBN  BF  LASER  LINE  DATA. 

C 

WRITE18.9500)  1 

WPITE(8.1000) 

READ(8,1010)FilEL 
READ18.  lOlOMiSERL 
C 

C 3HBULD  PRBGRAM  BE  RUN  BFF-LINE? 

C 

UPITE(8,8700) 

READ(8.8500)IIN 

IBAK6N-0 

IF(IIN.EQ.IYES)IBAKCN«1 

C 

C CBMPUT:  ABSBRBER  C8NCENTRAT18N(H0L*C«-2)  FBR  EACH  aBSBRBER 

c 

□B  310  I«I.NUnBR 

J-NSUB(I)  i 

W(J)«9.65726E20*PP( J)*PL/TE«P 
310  CBNTINUE 
C 

C CBMPUTE  EFFECTIVE  PRESSURE  FBR  EACH  ABSBRBER. 

C 

□B  300  I-l.NUHBR 
J-NSUBl  I ) 

PEI J)-(PRESS*(B(J)-1.  )*PPI J» )/760.  j 

300  CBNTINUE 

C i 

C CBMPUTE  CSl  ! 

C I 

CS1»<TEHP8-TEmP)/(TEMPB»TEmP«,6951)  ! 
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COMPUTE  CS2  and  CA  FOR  EACH  SUBSTANCE 


164 

165 

166 
167 
166 
169 
1>0 
1/ 

172 

173 

174 

175 

176 

177 
176 
179 
160 
181 
182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 
217 
216 


C 

c 

DO  410  I«l>7 

CS2( I l-(TEMPO/TEhP)*«SX(  I > 

CA( I l-{TEMpe/TEMP)»»CX(  I ) 

410  CONTINUE 
C 

C COMPUTE  DOPPLER  HALF-UIDTH  FOR  EACH  SUBSTANCE 

C 

DO  420  I»l>7 

ALFADm-3.5612E-'7«SQRT{TEHP/MU  > > 

420  continue 
C 

C POSITION  TAPES  PROPERLY 

C 

REWIND  0 

CALL  TPFILE(0,NFILE> 

CALL  TPKKD(0»NRKD> 

URITE(8«8800> 

REWIND  2 

CALL  TPFILE(2,NPFILE) 

CLOSE  2 
WRITE(8.8900) 

ENDATA-. false. 

FIRST  ■•TRUE- 
IARI-1 
C 

CCCCC  RFAD  THE  LASER  LINES  INTO  vLINE  FROM  FILEL>USeRL 
C 

CALL  CLSE(4> 

CALL  ASSIGN|FILEL»USERL.4) 

00  110  JF-1,125 
READ(4,8600IVLINEUF) 

IF( VLINEIJFI .EO. 0.0150  TO  l70 
110  CONTINUE 
170  CONTINUE 

IF( IBAKGN.EO.l )CALL  DEASSN 
CALL  ASNPLTI512.2.FILEP) 

DO  700  N»1»NUMPLT 
DO  500  I-l.lOOl 
Y( I )»0. 

500  CONTINUE 

DO  600  I-l.tOOl 
V-V1»( I-l >*nELVl 
C 

C CALL  data  to  make  SURE  THERE  IS  ENOUGH  DATA  TO  DO  THE 
c calculation. 

C 

CALL  DATAIV, round. VI, V2.IUAKGN, SLOW) 

C 

C CALL  A0SVO3  To  GALCUL'TE  THE  ABSf^RPTlON  COEFFICIENT. 

C 

CALL  A0SVO3IV, BOUND) 

Y(  I )«Y(  I I'ACAYW 
600  CONTINUE 
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219 

220 
221 

610 

222 

C 

223 

C 

224 

c 

225 

2'16 

227 

220 

620 

229 

230 

231 

C 

232 

C 

233 

c 

234 

235 

236 

237 

630 

230 

239 

240 

c 

241 

C 

242 

C 

243 

244 

245 

640 

246 

247 

C 

240 

C 

249 

C 

250 

251 

252 

650 

253 

C 

254 

C 

255 

C 

256 

C 

257 

C 

250 

C 

259 

260 
261 
262 

C 

263 

264 

265 

266 
267 
260 

59 

269 

270 

210 

271 

272 

273 

DO  610  1-1*1001 
Y( I )-ALOClO(Y( I » » 

CONTINUE 

FIND  BIGGEST  YU) 

YMAX--1  .E35 
DO  62C  1-1*1001 
IF(  Y(  1 ) .CT.YNAX)YMAX-VU  ) 

CONTINUE 
I YMAX-YMAX 

IF( YHAX.GT.O) lYMAX-I YMAX41 

FIND  smallest  YU) 

YHIN-1  .E35 
Of  630  1-1*1001 
IF(Y(  l).LT.YmN)YMlN-Y(I) 

CONTINUE 

lYMIN-YMIN 

1F(YHIN.LT.0)IYM1N-IYMIN-1 

SUBTRACT  lYMlN  FROM  ALL  YU> 

DO  640  1-1*1001 
Y( 1 >-Y( I >-IYMIN 

continue 

IYMAX-IYMaX-IYMIN 

NORMALIZE  Y FOR  PLOTTING 

DO  650  1-1*1001 
Y(  1 )-6.5*Y( I 1/IYMAX 

continue 

THE  NEXT  LOOP  PLOTS  THE  CALCULATED  DATA 


FIRST  A diversion  TO  GET  THE  PEN  WRITING 

CALL  P(  ■)T(0.0*-1.*3) 

CALL  PLeT(4.C*-I.*2) 

CALL  PL0T(0.0,-I.*2) 

DO  59  IN-1*1001 
X-0.01*t 1N“1) 

YP-YI  IN) 

IF( IN.EO.l )CALL  PL0T(X*YP*3) 

CALL  PL0T(X.YP,2» 

continue 

CALL  PLOT  lO.n. 0.0*3) 

IFIVLINEUARI  ).CT.V1+DELV1«1000.  ) CO  TO  220 

IF ( VL  INE< I API ) ,EO.O.n)GO  TO  220 
IF(  VL  INEU  ARI  ) .LT.Vl  )G0  TP  215 
X - (VLINEUARD  - VI  )/DEL\'l/100, 


L I 

i 


274 

CALL  PL6T(X.0.0.2) 

275 

CALL  PL6HX,6,5.2) 

276 

CALL  PLOT{X,0.0.3) 

277 

215 

IAR1-IAK141 

276 

ce  TO  210 

279 

220 

CALL  PL6T1  10.0.0..-'  < 

280 

CALL  PLaT<0..6.5,3) 

261 

230 

CALL  PL3Tno.0,6.5,2) 

262 

C 

283 

C 

plot  title  INF6RHATI6N 

284 

C 

265 

CALL  SVneoL (0.0. 6.8,0. 1 . IShTEMP  (DEG  PM  .Q.Q.ISI 

266 

CALL  NUnOEP( 1.50,6.6.0.1 .TpLOr. 0.0,1 1 

267 

CALL  SVMBsL(C.0,7.0,0.1.l5MAh3UNT  (T8(TPM  .0.0.151 

266 

call  symbol ( 1 .50.7.0.0.1.T1 rLE».r.0.42> 

269 

CALL  SYMBeL(0.0.7.2,0.|,16MA8SPP?ER3l  .0.0.161 

290 

CALL  SYMB6L( 1.50.7.2,0.1, title A. 0.0.421 

291 

CALL  SYMBOL (7. 1 ,7.0,0. 1 ,23hT6TAL  PRESSURE  (TORRM  .0.0*231 

292 

CALL  CHAMCEIPPLOT.PUMI 

293 

CALL  SYH8oL(9.4.7.0,0.1.DUH.O.O,5> 

294 

CALL  CETDaT ( IDatEi 

295 

CALL  SYMB8L(7.1,7.2,0.1. IDaTE,0.0,9) 

296 

C 

297 

ccccc 

296 

DYS-IYMAX/6.5 

299 

SPC-1  .0 

300 

XNT«100.»DELV1 

301 

IP(XNT.LT..7)SPC-2.0 

302 

CALL  XXIS(0.0.O.O.10HWAVENUMBEP,“10,i0.0.0.O,V1,'XNT.1C.0, 

303 

CALL  AXIS<O.O.O.O.IY?.-1,“10.0,C.0.0.0.1 .O.SPC.l 1 

304 

CALL  FACT6RI0.5I 

305 

CALL  AXIS(0.0.0.0.IY2.-1.-20..0.C.0.0,1.0,0.4, 1 ) 

306 

CALL  FACTeR(l.O) 

307 

YFIN«10.*.IYHIN 

306 

CALL  L5AXS(lC,,0..|Y2.  -1  . -«> .5. 90 .0. YM I N. n YS  1 

309 

CALL  LGAXS(0.0,O.C.1BHAPS.  C6EFF.  ( KM- 1 > . 1 8. 6 .5,90 .0. YM I N, 1 

310 

CALL  PI-OT(  12.0,0.0,-31 

311 

CALL  PLOTlO. 0,0. 0.999) 

312 

CALL  ASNPlT(5i2.2,FILEP) 

313 

V1»V1+I000.*DELV1 

314 

c 

315 

c 

WRITE  message  in  MSG.3271A  TELLING  DATE  AND  TIME  EACH 

316 

c 

CALCULATION  U*S  COMPLETED. 

317 

c 

318 

CLOSE  7 

319 

CALL  ASSIGNI4HMSC  ,5H3271A,7) 

320 

CALL  GETTIMIITIMI 

321 

WRITE  1 7,9000  In, IT IM.IOATE 

322 

CLOSE  7 

323 

700 

CONTINUE 

324 

CALL  EXIT 

325 

C 

326 

C 

FORMATS 

327 

c 

328 

950 

FORMAT!*  where  IS  THE  ABSORPTION  LINE  DATaT*  1 

174 


I 

I 


4 


329 

9400 

F8RMAT ( * 

FNTER  tape  start  pOINTI  file  NUMBER  AND  RECORD  NUMBER 

330 

9500 

FORMAK' 

where  is  the  lASER  LINE  DATA?*) 

33l 

1000 

FORMAT { » 

ENTER  FILE  AMD  USE«  NAME  ON  TWO  LINES.*/) 

332 

1010 

F0RMAT(a6) 

333 

1011 

F0RMAT(7A6) 

334 

1 100 

FORnAK' 

enter  .mto  or  ,mti*) 

335 

1200 

F3RMAT(* 

WHICH  TAPE  UNIT  WILL  PLOT  DATA  Be  WRITTEN  ON?*) 

336 

1300 

FORMAK* 

ENTER  number  OF  FIRST  FILE  TO  BE  WRITTEN.*) 

337 

C 

NUMBR  OF 

absorbers  • NUMBR 

336 

c 

NUMBER  OF  PLOTS  ■ NUMPLT 

339 

2000 

FORMATC 

FNTFR  THE  NIIMPFR  oE  ABSORBERS  AND  THE  NUMBER  OF  PLOTS 

340 

3000 

FORMAT ( * 

FNTER  THE  beginning  WAVFNUMBER*/ 

341 

1*  AND  ThF.  UaVENUMBER/PLOT.*  ) 

342 

4000 

FORMAT ( * 

enter  the  total  pressure  in  TORR.*) 

343 

4 100 

FORMAT(* 

enter  id  number,  partial  pressure.  'iND  broadening*/ 

344 

1 

coefficient  for  each  absorber. */) 

345 

4 200 

FORMAT ( * 

enter  absorber  description  information  on  two  lines*/ 

346 

1 

OF  UP  TO  42  characters  E*CH.  THE  FIRST  LINE  WILL*/ 

347 

2 

appear  immediately  TO  THE  RIGHT  OF  ABSORBFRSI  .*/ 

346 

3 

THE  SECOND  LINE  WILL  APPEAR  TO  THE  RIGHT  OF*/ 

349 

4 < 

AMOUNT ! TORR » : AND  IMMEDIATELT  BELOW  THE  FIRST  LINE.*/ 

350 

4900 

FORMAT ( * 

enter  id  number  of  ABSOPSER  no.* 11/ 

35i 

1 

1-H20  2-C02  3»03  4«N20  5»C0  6-CH4  7»02*/> 

352 

4950 

FORMAT (' 

FNTER  NUMBER  OF  DESIRED  ISOTOPE.*/ 

353 

1 

IF  all  ISOTOPES  AoE  DESIRED  ENTER  Q.*) 

354 

5000 

F0RmaT(* 

enter  the  ABSORBER  AMOUNT  1 IJ  TORR.*) 

355 

5500 

FORMAT ( * 

ENTER  THE  .SflF-BR»ADENING  COEFFICIENT.*) 

356 

5600 

FORMAT ( * 

enter  the  path  length  in  meters.*) 

357 

6000 

FORMAT ( * 

FNTER  THE  TEMPERATURE  !►'  DEGREES  F.*) 

356 

6100 

FORMAT (* 

CO  YOU  WISH  A CONTINUUM?*) 

359 

6200 

FORMAT(* 

ENTER  The  coefficients  fop  the  polynomial*/ 

360 

1 

K ! V 1 •At)4Al*v*/2*V**2.  W)ERE  K IS  In  KM-1  AND*/ 

361 

2 

V IS  UAVENUM3EIT.*/ 

362 

3 » 

aO-  *) 

363 

6300 

FORMAT(* 

Al«  *) 

364 

6400 

FORMAT!* 

A2«  *) 

365 

6500 

FORMAT!* 

enter  frequency  LIMITS  FO  CONTINUUM.  VL  AND  VH*/) 

366 

3500 

F0KMAT!A 

3) 

367 

8600 

FORMAT!F 

12.0) 

366 

8700 

FORMAT!* 

DO  TOU  WISH  TO  RUN  THIS  PROGRAM  IN  BACKGROUND?*) 

369 

8800 

FORMAT ! ' 

The  Data  Tape  is  positioned.*) 

370 

8900 

FORMAT  !' 

The  Plot  tape  is  positioned.*) 

371 

9000 

FORMAT  ! ' 

PLOT  '13*  finished  at  *3a3*.  *3a3) 

372 

END 

373 

C 

374 

C 

THE  FOLLOuING  SUBROUTINE  RFAOS  DATA  FROM  TAPE  AND 

375 

C 

MODIFIFS 

IT  FOP  USE  IN  THE  program.  IT  ALSO  CHECKS  BEFORE 

376 

C 

EACH  CALOL'LATION  is  HADE  TO  MAKE  SURF  THERE  IS  SUFFICIENT 

377 

C 

DATA  TO 

PERFORM  THE  CALCULAflON. 

378 

C 

379 

subroutine  nATA(V.B0UND.Vl.V2.IBAKGN.3L9U) 

380 

LOGICAL 

endata. first 

381 

C0MM0N/Z/aD!3.250U).NABS!250C). IR . 15.16 

382 

C0MM0N/D/ENDATA.FIRST.NSUP!7».NU“'BR.cS1.CS2!7».CA!7».iS0T0P!7) 

\ 
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|riyHg 


363 

364 

365 
36b 
367 

366 
36V 

390 

391 

392 

393 

394 

395 

396 

397 
396 

399 

400 

401 

402 

403 

404 

405 

406 

407 
406 

409 

410 

411 

412 

413 

414 

415 

416 

417 
416 

419 

420 

421 

422 

423 

424 

425 

426 

427 

428 

429 

430 

431 

432 

433 

434 

435 

436 

437 


c 

C IF  THIS  is  the  first  TIME  THIS  St'3R6UTlNE  HAS  BEEN  CALLED 
C G6  IMMEDIATELY  T6  INPUT. 

C 

IFIFIRSTIGO  T6  500 
C 

C IF  THERE  IS  ENOUGH  DATA  IN  ADU.K)  AND  NAPS(K)  TO  PERFORI’ 

c the  calculation  at  this  frequency,  or  if  he  end  of  The 

C AFCRL  DATA  TAPE  HAS  BEEN  REACHED,  RETURN. 

C 

IFI  (ADI  I.IH  I .GT.IV^BOUNDn.OR.ENDATA)  RETURN 
C 

0 IF  MORE  data  IS  NEEDED  TO  PERFORM  THE  CALCULATION  AT  THIS 

C FREOUFNCY,  FIRST  MOVF  THE  dATa  IF  AD(J.K)  AND  MARS(K) 

C UHICH  IS  still  needed  from  MICH  K TO  LOU  K. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
DO  100  I-IH.l.-l 

1F(  At(  1,  1 1 .LT,  (V-BOUNOnCO  TO  200 
100  CONTINUE 
200  I«1A1 

K«1 

DO  400  J-I.IH 
DO  300  M-1,3 
AD(M.KI*AD<M,JI 
300  CONTINUE 

NABS(KI-NA8SUI 

K«K*1 

400  CONTINUE 
GO  TO  600 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
500  K«1 

FIRST". false. 
vlou«vi-bound 

VHIGH»V2*5OUN0 
GO  TO  600 
C 

C NOU  READ  THE  DATA  INTO  AD<J*K)  AND  NABS(K)  UNTIL  ADIJ.K) 

C IS  FULL  OR  enough  DATA  HAS  BEEN  READ  TO  COMPLETE  ALL 

C THE  DESIRED  CALCULATIONS  OR  THE  END  OF  THE  AFCRL  DATA 

C TAPE  has  been  reached. 

C 

C FOR  AN  ABSORPTION  LINF  K,  aD(1,K)  Is  TH^  FREQUENCY, 

C AD(2,KI  IS  the  line  strength,  AD(3,K)  IS  THE  K-iLF-UlDTH, 

C EPF  IS  energy  of  the  LOUER  STATE  OF  THE  TRANSITION, 

C ISTOF  IS  THE  ISOTOPE  I DEN r 1 F I c AT  1 ON,  AND  NABS(K) 

C IS  THE  substance  IDENTIFICATION  NUMBER. 

C 

590  CLOSE  0 

600  READ(0.1000,END»590) AL'(1,K),AD(2,K).aD(3.K»,EpP,I3T0P,NABS(K) 

C 

C DOES  THIS  line  HAVE  A LOUER  FREPUENCY  THAN  JS  NEEDED 
C FOR  THE  CALCULATION,  OR  IS  THE  LINE  STRENGTH  TOO  LOU? 

C 

IF<(AD(l,Kl.LT,(V-BeUNDn.0B.(AD(2,I).LT.SL0UnG0  TO  600 
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43«  - 

439  C 

440  C 

441 

442  C 

443  c 

444  C 

445  C 

446 

447  C 

448  C 

449  C 

450  690 

451 

452  700 

*;^3  710 

454 

455  790 

456 

457  800 

458  C 

459  C 

460  C 

461 

462 

463  C 

464  C 

465  C 

466  810 

467  C 

468  C 

469  C 

470 

471 

472 

473 

474 

475  880 

476  890 

477  900 

478 

479 

480  C 

481  C 

482  C 

483  C 

484 

485 

486  loco 

487  2000 

488 


IS  THIS  THE  END  OF  TIE  DATA  TAPE? 

IF(AD(1.K).E0. 12075. 863)<:9  TO  79C 

IIF(ADM.K)  .CT.VHICHIGO  to  880 

IS  THIS  substance  to  be  considered  in  the  calculation? 

D8  700  I-l.NUMBR 

IF ( NSUP ( I ) .eC.NABSIK ) )G0  TO  800 

CONTINUE 

IF( .NOT.ENDaTaICO  to  600 

GO  TO  890 
ENDATA».TRUE. 

GO  TO  690 
J-Na8S(K) 

IS  THIS  ISOTOPE  TO  BE  CONSIDERED  IN  yHE  CALCULATION? 

IF( ISOTOPI J1 .EO.OIGO  to  810 
IF( ISOTOPI J) .NE.ISTOPIGO  TO  7lO 

CORRECT  THE  LINE  STRENGTH  FOR  TEPPERATURE. 

AD(2.K)«AD(2.K)«CS2( J)*EXPf-EPP»CSl ) 

CORRECT  THE  HALF-WIDTH  FOR  TEMPERATURE. 

ADI3.K l■AD(3.K)*CA^ J) 

IF(ENDATA)GO  TO  900 
K»K*1 

IFIK.LE. 2500)00  TO  600 
GO  TO  890 
ENDATA-.TRUE. 

K-K-1 

IH«K 

I5«l 

I6«t 

Ir  program  is  on-line  TYPE  MESSAGE  ON  THE  TELETYPE  THAT 

data  mas  been  read  in.  ■ 

IF(  IBAKGN.EQ.O) WRITE  18, 2000) 

RETURN 

FeRMAT(2E10.4,F5.3.F10.3,38X.I4,l3) 
formati'one  set  of  data  has  Been  Read.*) 

END 


489 

490 

491 
49? 


cccccccccccrccrrcccccccccr-cccccccccccrccccccccccrcrccrcrccccccrcccc 

C THE  FOLLOWING  SUBROUTINE  CALCULATES  THE  ABSORPTION 

C COEFFICIENT. 


i 


493 

494 

495 

496 

C 

497 

C 

498 

C 

499 

C 

500 

C 

501 

502 

503 

C 

504 

505 

10 

506 

20 

507 

506 

30 

509 

50 

510 

40 

511 

512 

513 

514 

515 

55 

516 

C 

517 

C 

518 

C 

519 

520 

521 

522 

C 

523 

C 

524 

C 

525 

526 

527 

C 

528 

C 

529 

C 

530 

C 

531 

532 

c 

533 

c 

534 

c 

535 

c 

536 

c 

537 

c 

538 

c 

539 

540 

541 

542 

c 

543 

c 

544 

c 

545 

546 

c 

SUBROUTINE  ABSV03I V»BOUND> 

COMMON /Z/ ADI 3»2500)»NABS I 25OO) » IM» 15. 16 

CeMMON/A/PE(7l,W|7)»Pi..AO»Al»A2»VL.VH.CAYW»Al.FAD(7) 

THE  FOLLOUING  loops  nETFRMlNE  WHICH  ABSORPTION  LINES 

WILL  BE  USED  TO  CALCULATE  THE  ABSORPTION  COEFFICIENT  AT  THIS 

FREQUENCY. 

DO  20  I-I5.IH 

IF(V-B0UND-AD(1,I) I10.10.20 
I5-I 

CO  TO  30 
CONTINUE 
DO  50  K-I6.IH 

IFI  V+BOUND-ADn.KnAO.50.50 
CONTINUE 
I6-K-1 
CAYl-0.0 
CAY2-0.0 
CAY-0.0 
DN2-AL0G(2.) 

DN1-SQPT(DN2> 

THE  FOLLOUING  LOOP  SUMS  THE  CONTRIBUTIONS  FROM  ALL  THE 
ABSORPTION  LINES. 

DO  60  I-I5.I6 
J-NABSI I ) 

Y-ABSIV-ADU.m 

CORRECT  the  DOPPLER  HALFUIDTH  FOP  FREQUENCY. 

ALFD-ADI 1. I »*ALFADI J) 

xv-dm-y/alfo 

THE  FOLLOUING  STATEMENT  CORRECTS  THE  LORENTZ  HALF-UlDTH 
FOR  PRESSURE. 

ALFL-ADI3.I )*PE(J» 

IF  THE  ABSOPTTION  LINE  IS  FAR  AUAY  FPOM  CnLcULATION 
FREOUEMCY,  USE  LOPENT?  SHAPE  REGARDLESS  OF  pRESSURF.1 
BUT  IF  LORENTZ  HALF-WIDTH  iS  LESS  THAN  FIVE  TIMES  THE 
DOPPLER  UIDTH.  USE  THE  VOIcT  LINE  SHAPE.  OTHERUISE 
USE  THE  LORENTZ  SHAPE- 

IF(XV.CP.3()0.)  GO  TO  59 
IFIALFL/ALFD.GE.5.IG0  TO  5’ 

YV-ALFL/ALFD«0N1 

CALCULATE  THE  VOIGT  SHAPE 

SUMS- AD ( 2. I )*UIJ>«DNI/(S0RT(3.1‘I1592)»ALFD)*V0IGT(XV.YV) 
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547 

540 

549 

59 

550 

C 

551 

C 

552 

c 

553 

60 

554 

555 

556 

557 

C 

558 

C 

559 

C 

560 

C 

561 

C 

562 

70 

563 

564 

565 

80 

566 

c 

567 

C 

568 

C 

569 

570 

c 

57l 

c 

572 

c 

573 

c 

574 

c 

575 

c 

576 

c 

577 

578 

579 

580 

581 

ccccc 

582 

c 

583 

c 

584 

c 

585 

586 

c 

587 

c 

588 

589 

590 

591 

592 

10 

593 

594 

595 

ccccc 

596 

c 

597 

c 

598 

c 

599 

600 
601 

c 

CAY-CAV+o.3ie3*SUi13 
60  TO  eo 
IF(Y-3. »60,60.70 

CALCULATE  THE  LORENTZ  LINE  SHAPE  FOR  V NFAR  LINE  CENTER 

SUMl-AD(2.n*AD(3.i)/(Y.Y*AD(3.I)«AD(3*I)*pE(J)*PE(J)) 
SUMl'-SUHi  *pE(  J)«W(  j)*o.3t03 
CAY1«CAY14SUM1 
GO  TO  00 

CALCULATE  THE  LORENTZ  LIME  SHAPE  FOR  V FAR  rRnll  LINE 
CENTER.  IGNORING  THE  HALF-uIDTH  SQUARED  TERM  IN  THE 

denominator  For  increased  calculation  speed. 

SUM2*AD(2.I»*AD(3,I)/(Y*Y) 

SUM2*SUM2*PE(J»*u(J»*0.3103 

CAY2«CAY2-»SUf12 

CONTINUE 

SUMMATION  OVER  ALL  LINES. 

CAYU«CAY14CaY2+CaY 

the  following  ALLOWS  A CONTINUUM  UHIPM  IS  A FUNCTION  OP 
FREQUENCY  TO  BE  ADDED  TO  ThE  LOCaL  UNC  ARSoRmirjp^TRUN 
THE  polynomial  WHICH  EXPRESSES  THE  CONYINUU**  IS  XIVI,  WHERE* 

i;  zv. 

IF((V.LT.VL».0R.(V.CT.VH) ) return 
CAYU«CAYW+(AC4AI»V4A2*V*V»*PL*1 .E-3 
RETURN 
END 


THE  FOLLOWING  SUBROUTINE  PRINTS  DATE 
SUBROUTINE  CLOCK 


AND  time  on  the  TELETYPE. 


DIMENSION  ITIME(3»,IDATE(3) 

CALL  GETDaTI IDATEI 
CALL  GETTIMI ITIMEI 
URITE(8.10)IDaTE.ITIME 

FORMAT!*  date  *3a3»  TIME  *3A3/) 

RETURN 

END 

subroutine  chance  for  F10.4  To  AS  CONVERSION 
FOR  NUMBERS  IN  THE  RANGE  9999.9999  TO  .0001 

subroutine  CHaNGEIFN'JH,  A5) 

INTEGER  FlVAl  (5),A10(10».SCR(<*> 


602 

603 

604 

605 

606 

607 

608 

609 

610 
611 
612 

613 

614 
6!5 
616 

617 

618 

619 

620 
621 
622 

623 

624 

625 

626 

627 

628 

629 

630 

631 

632 

633 

634 

635 

636 

637 

638 

639 

640 

641 

642 

643 

644 

645 

646 

647 

648 

649 

650 

651 

652 

653 

654 

655 

656 


C FIND  DECIMAL  POINT  LOCATION.  LOC  UlLL  BE  NUMBER  OF 

C DIGITS  TO  left  OF  DEC  IF  POSITIVE. 

C 

LOC- ALOGlCi(FNuH)*  1.00001 
LFTL0C-5-L0C 
1F(lFTLOC.GT.5)LFTLOC-5 
ENCODE  I lO.lOOO.ScRlFNUM 
DECODE!  1C, 20CCi. SCR)  AlO 
C 

C PUT  THE  5 CESlRED  CHARACTERS  INTO  FIVAl 

C 

DO  100  1-1-5 
J-LFTLOC*! 

FIVaI  I 1 )-AlO( J) 

100  CONTINUE 

C 

C 

C CHANGE  FROM  5A1  TO  A5 

C 

ENCODE (6>300C.SCR)FIVA1 
DEC0DE(6.40QG.SCR)A5 
1000  FORMAT(F10.4) 

2000  FORMATdOAl) 

3000  FORMATISAl) 

4000  F0RMAT(A5) 

RETURN 

END 

C 

C THE  FOLLOUING  SUBROUTINE  EVALUATES  THE  VOIGT  PROFILE 

C EXPRESSION. 

C subroutine  WAS  PROVIDED  BY  CHARLES  YOUNG. 

C 

FUNCTION  VOIGT  (XIN.YIN) 

DIMENSION  RA(32),CAI32).RB|32),CPI32),9I44),AKI5),AMI5)»DYI4) 
1>HH(2)-XX|2I-A(42) 

data  HM  /. 8049141. .8l3l2e3E-Cl/, XX  / .5246 4 76, 1 . 650680/, A 

1 /O.O, .2,C.,-.184.0.0,.l55S4,0.0,-.l2l664,r).0».fl7708l6E-l,0.0*-.5 

2 851412E-1 ,C.  ,.3621573E-1»0.0,-.20P4O76e-1»0.0,.1 I 1960lE-l,0.C,- 

3 ,5623190E-2,C.C,  .?648763F-2,0.r,-.i  1732678-2 » 0 . 0 ,. 4 ?,99520E-3,C .0 

4 ,-. 1 53363 1 E-3.0.C,.7?28775e-4,C.C,-. 256555 1E-4,0.0».8662C74E-5»0. 

5 0»-.278763eE-5,C.C,.fl566674E-6,0.0,-.25l8434E-r»C.O,.7093602E-7/ 
X-XIN 

Y-YIN 

X2-X-X 

Y2-Y*Y 

IF  (X-7.)  200,201,201 
200  IF  (Y-1.)  202,202,203 
203  RA(l)-0. 

CA(l)-0. 

RB(1)-1. 

CB(1 )-0. 

RA(2)-X 

CA(2)-Y 

RB(2;-.5-X2aY2 

CB(2)--2.*X-Y 
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CBI-CB)2) 

UVl-0. 

D9  25C  J«2»3l 

JMINU9-J-1 

JPLUS-J+1 

FLOATJ-JhINUS 

RB1»2.*FL8ATJ*RB<2) 

RAl"-FLPAT.r»(2,«Fl9ATJ-l  .)/2. 

RA(  J + 1 )“Pel*RA(J)-CB1*CA< J)-M?a1«RA(J-1) 

CA(  J*1  J)+CD1*RAU)+RA1«CA(  JmINU3> 

RB ( J*1  )-PBl •F9( J)-CP1 *C9( JI+Ra  1*RB(J-1) 

CB( J+1 )"R?1 ‘Cat J)+CB1 *RS( J)+Ra1 *CB( J-1 ) 

UV-(CA( JPLUS>*RB( JPLIIS)-PAt JPLUS)*CB( JPLUS) )/(RB( JPLUS>**2* 

> CB(  JPLUS)«CP(JPLl'S>) 

IF(ABS<UV-UV1 >-1.E-6)25I,250,250 

250  UV1»UV 

251  V9IGT‘^UV/1. 772454 
RETURN 

202  IF( X-2.  )301#301»302 

301  aINT  - 1. 

MAX-12. +5. «X2 

KMAX-MAX-1 

K9-0 

D8  303  K-K9.KMAX 
AJ-MAX-K 

303  AINT-AINT«(-2.*X2)/(2.*AJ*1.1*1. 

U«-2-*XaAINT 

G9  T9  304 

302  IF  (X-4.5)  305,306,306 

305  9(43)«0. 

‘3(44)«0. 

J«42 

D9  307K-1,42 

B(J)-.4«X«9(J'U)-BlJ.t2>+AIJ> 

307  J-J-1 
u-8(3)-B( 1 » 

S9  T9  304 

306  AINT-1.0 
MAX-2. ♦40. /X 
AMAX-MAX 

D9  308  K»1,MAX 

AINT-AJN'*(2.*AMA'.-I.)/(2.«X2)  + 1. 

308  AMAX-AMAX-l. 

J--AINT/X 

304  V-1 .772434*EXP(-X21 

H-.02 

JM-Y/H 

IF(JM)310, 311,310 

311  H-Y 
310  E-0. 

L-0 

0Y( 1 )-0. 

312  DY(2)-H/2. 

UY(3)-DY(2» 

0Y(4)«H 


3la  AK(1I>0. 

AHfn-O. 

00  313  J-1.4 
YY-Z*mj) 

UU"U*.5*AK{J» 

VV«V*.5*AM{Jt 

AKlJ  + n-  2.*<VY*UU*X«VV»*H 
Ah< J*1  )»-2.*{i.  + x*UU-YY*VV)*H 
IFU-3)  313,314,713 

314  AK(4I-2.*AK(4| 

AM{4)«aM{4I*AM(4! 

313  CONTINUE 

L-C  + l 

U"U*.1566667*(AKC2)+2.*AK(3)-»AK{4)*ak  (5)  ) 

V«V*. 1665667*  <AM<2 ) + *M<1)*aM(3)*AN{ 4 ) ♦AM151 ) 

IFIJH)  315,320.315 

315  IF<L-JN)  318,317,320 
317  AJM-JM 

M»Y-AJH*H 
GO  TO  312 

320  VOlGT-V/1 . 772454 
RETURN 
201  F1«0. 

no  330  J«l,2 

330  F1“F1+HH<J)/<Y2*(X-XX'.  Jl)*(X-XX<J)))*HH(J)/{Y2*{x*XX{3>)*{X*XX{J)) 
1) 

V0IGT«Y*Fl/3. 1415927 

RETURN 

END 
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APPENDIX  C 

OPTICS  DESIGN  PROGRAM  LISTING 


1 C r.  s,  HILLS  NOVEMJER  2,197? 
i C 

3 C FHlOl 

•*  C 

5 c la^er  Spot  size  program 

6 C 

7 C This  program  uses  the  FGUaTIONS  on  which  the  COLLINS 

8 C CHART  IS  bA'iEO.  RLFERENCE:  S,  a.  COLLINS,  jR., 

9 C AFPL.  OPT.  3,  1263(1964) 

10  C 

11  c This  program  calculates  the  spot  size  for  a plane- 

12  C spherical  resonator,  up  to  )0  focusing  mirrors  or  lenses  mat 

13  C PE  PLACEO  ANTWHEHE.  THE  OUTPUT  OK  THE  LASER  IS 

14  C CONSIOFHEU  To  PE  from  the  Si'ithiCAL  mIRRoR. 

15  c distances  are  specified  with  respect  to  the  plane  mirror. 

16  c the  user  specifies  the  Radius  of  curvature  o''  the  output 

17  c mirror,  its  Index  of  refraction,  the  mirror  separation, 

18  C ANO  THE  LOCATION  ANO  FOCAL  LENGTH  OF  THE  FOCUSING 

19  C elements, 

20  C when  ALL  parameters  ARE  SPECIFIED,  THE  PROGRAM  WILL 

21  C TTPE  'ENTER  LOCATION' 

22  C ANO  THE  USER  TYPES  THE  LOCATION  OF  A POINT  OF  INTEREST 

23  C WiTh  RESPECT  TO  THE  PLANE  MIRROR  OF  THE  L»SER.  ThE  PROGRAM 

24  C types  the  spot  SIZE  ANO  THE  OISTANCE  TO  THE  FOCUS 

25  C OF  the  beam.  THE  PROGRAM  THEN  TYPES  'ENTER  LOCATION'  AGAIN 

26  C ANO  the  user  MAY  PROCEEO, 

27  C 

28  dimension  BI(12),ZM1(12(,ZM2(12),FM(10) 

29  CALL  ESC(SIO) 

30  ZM1(1)=0, 

31  Zm2(1)=0. 

32  10  URITE(8,20) 

33  20  FoRMA'M/'LASeR  SPOT  SIZE  PROGRAM'/// 

34  1'ENTeR  output  mirror  RADIUS,  INDEXt  MIRRoR  SEPARATION'/) 

35  REA0(8,-)  B1,RN,Zm1(2) 

36  URITE(8,3S) 

3T  35  FORMAT) 'ENTER  WAVELENGTH'/) 

38  REA0(B,-)WAVE 

39  WRITE(B,25) 

40  25  FqRMAT( 'enter  NUMBER  OF  ELEMENTS'/) 

41  R£A0(S>-)N 

42  00  33  1=1, N 

43  WRITE(8.30)I 

44  30  FoRMAT( 'ENTER  POSITION,  FOCAL. LENGTH  OF  ELEMENT  NO,  ' 

45  1,I2/) 

46  REA0(8,-)ZM1(I+2),FM(D 

47  33  continue 

48  C 

49  c PI  IS  OUTPUT  Mirror  radius  of  curvature 

50  c Rn  is  index  of  refraction  of  output  mirror  material 

51  C Zm1(2)  is  mirror  SEPERATioN 

52  C ZmI  is  position  OF  FOCUSING  mIRROR  OR  LENS 

53  c Wave  is  the  wavelength  of  thf  laser 

54  c FM  is  focal  length  of  focusing  mirror  or  lens 

55  C 


I 


5G 

C 

57 

C 

58 

59 

60 

C 

61 

C 

62 

6* 

C 

6R 

C 

65 

C 

66 

67 

C 

68 

C 

69 

C 

70 

71 

72 

C 

7S 

C 

7H 

75 

c 

76 

c 

77 

78 

79 

c 

80 

c 

81 

82 

c 

89 

c 

8H 

85 

c 

86 

c 

87 

88 

c 

89 

c 

90 

91 

c 

92 

c 

99 

9H 

c 

95 

c 

96 

97 

c 

98 

90 

99 

c 

100 

c 

101 

c 

102 

60 

109 

109 

70 

105 

C 

106 

C 

107 

108 

109 

110 

C 

THt  intersection  POINT  OF  IHF  Hspl  LINE  aNO  THE 
constant  OISTAMCE  circle  DESCRIBEO  by  ZxTl  Is  GIVEN 

PY  the  coordinates 

VIsSORTI  I l./2./Z«H2>  )*»2I 

Yisi./ni 

the  a CIRCLE  WHICH  also  INTFoSECTS  this  point  is  given  bt 
Blll|slXl**2TYI»*?»/2./Xt 

the  output  mirror  acts  like  a negative  lens  of  focal 
length  F1=-8I/(RN-1. > 

F13-B1/(NN-1.> 

the  lens  effect  of  the  mirror  causes  Y1  to  be  DEcREASEO 
BT  l./Fl.  XI  remains  the  SAME.  JUST  OUTSIOE  THE  MIRROR 
X2SX1 

Y2rTl-l./Fl 

THE  NEW  B CIRCLE  IS  GIVEN  BT 
HI(2)s(X2**2*Y2**2>/'2,/x2 

the  new  z reference  is  given  by 

2m2(2)=T2/(X2»*2+Y2**2) 

no  NO  I=1«N 

at  the  FOCUSINi;  element,  alpha  IS  GIVEN  pY 

AlPHAsAtAN(2.*(ZM1II  + 2»-tM1| I*i)+2M2I !♦!) )*BI I I + U > 

at  the  focusing  EuEMENT  the  Y-VaLUE  is  given  by 

TM=BHIT1)*SIN(2.*ALPHA 

ANO  THE  X-VALUE  IS 
xm*c;;:*htbhiti»*cosi2.*alphA) 

the  focusing  element  oecreasfs  YM  by  X./FM 

YM*TM-X./FM(II 

The  new  b circle  is  given  by 

BI(IT2»=<XM**2+TM**2»/2,/XM 

the  new  z reference  is  given  by 

Zrt2(IT2»MTM/IXM**2+YM**?) 

continue 

THE  calculation  PORTION  STARTS  HERE 
WRITE(e.TO) 

FoRMAK/’ENTER  location*) 

REAOIB.-IZ 

check  to  SL'-.  where  in  THE  OPTICAL  SYSTEM  Z IS 

MsN*l 

00  7S  1 = 1. M 

IF(Z.GT.ZM1(ItI))G0  To  75 
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APPENDIX  D 

MASS  SPECTROMETER  MEASUREMENTS  OF  HDO  CONCENTRATION 


A small  mass  spectrometer  residual  gas  analyzer  was  used  to 
monitor  the  concentration  of  HDO  relative  to  H2O  in  the  absorption 
cell.  First  the  mass  spectrometer  was  baked  out  using  a heat  lamp 
and  heating  tape  to  remove  as  much  residual  water  as  possible. 
After  baking,  the  pressure  in  the  mass  spectrometer  system  was 
10^  torr.  Next  a sample  of  water  containing  6.22%  HDO  was 
admitted  to  the  absorption  cell  by  completely  evaporating  a 
small  amount  of  the  enriched  water  contained  in  a bottle,  and 
nitrogen  was  added  to  a total  pressure  of  760  torr.  The  sample 
was  allowed  to  mix  for  one  hour,  and  then  a mass  spectrum  of  the 
sample  was  recorded.  The  portion  of  this  spectrum  between  16  and 
19  is  the  bottom  curve  in  Fig.  75.  The  sample  was  admitted  to 
the  spectrometer  through  a small  leak  valve  and  continuously 
pumped  out  through  the  diffusion  pump.  The  spectrum  was  recorded 
at  a pressure  of  approximately  10“5  torr.  In  Fig.  75  mass  number 
16  corresponds  to  oxygen,  17  to  OH,  18  to  H2O  and  OD,  and  19  to 
HDO. 


Next  the  absorption  cell  was  closed  off  from  the  mass  spectrom 
eter,  and  the  sample  was  allowed  to  remain  in  the  cell  and  mix  over 
night.  During  this  time  the  mass  spectrometer  was  again  baked  out 
in  order  to  remove  the  residual  water  from  the  first  sample.  After 
baking,  the  mass  spectrometer  was  again  at  a pressure  of  10"®  torr. 
The  next  day  another  mass  spectrum  was  recorded  the  same  way  as  the 
first.  The  portion  of  this  spectrum  from  mass  number  16  to  19  is 
the  top  curve  in  Fig.  75. 

The  amplitude  of  the  peaks  at  mass  number  18  and  19,  cor- 
responding to  H2O  and  HDO,  is  greater  in  the  second  spectrum  than 
in  the  first.  However,  in  both  cases  the  ratio  of  the  amplitude 
of  the  peak  at  18  to  the  amplitude  of  the  peak  at  19  is  about  10 
or  12.  Therefore  the  ratio  of  HDO  concentration  to  H2O  concentrati 
did  not  change  although  20  hours  elapsed  between  the  time  the 
first  spectrum  was  recorded  and  the  timv->  the  second  spectrum  was 
recorded.  The  conclusion  then  is  that  the  absorption  cell  walls 
do  not  absorb  one  water  isotope  more  than  the  other. 


It  was  not  possible  to  determine  the  absolute  concentration 
of  HDO  compared  to  H2O.  It  was  only  possible  to  determine  that 
the  concentration  did  not  change  while  the  sample  was  in  the 
cell.  The  absolute  concentration  must  still  be  determined  from 
the  amounts  of  D2O  and  H2O  used  to  mix  the  water  sample,  and  the 
known  equilibrium  constants. 
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